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Abstract 



In this thesis we study the dynamics of higher- dimensional gravity 
in a universe emerging from a brane collision. We develop a set of 
powerful analytic methods which, we believe, render braneworld cos- 
mological perturbation theory solvable. Our particular concern is to 
determine the extent to which the four- dimensional effective theory 
accurately captures the higher- dimensional dynamics about the cos- 
mic singularity. 

We begin with a simple derivation of the low-energy effective action 
for braneworlds, highlighting the role of conformal invariance, before 
showing how the effective action for a positive- and negative-tension 
brane pair may be improved using the AdS/CFT correspondence. 

We then solve for the cosmological perturbations in a five-dimensional 
background consisting of two separating or colliding boundary branes, 
as an expansion in the collision speed V divided by the speed of 
light c. Our solution permits a detailed check of the validity of four- 
dimensional effective theory in the vicinity of the event correspond- 
ing to the big crunch/big bang singularity. We show that the four- 
dimensional description fails at the first nontrivial order in (V/c) 2 . 
At this order, there is nontrivial mixing of the two relevant four- 
dimensional perturbation modes (the growing and decaying modes) 
as the boundary branes move from the narrowly-separated limit de- 
scribed by Kaluza- Klein theory to the well-separated limit where grav- 
ity is confined to the positive-tension brane. 

We highlight the implications of this result for cosmology, in particular 
for the propagation of a scale-invariant spectrum of density pertur- 
bations across the bounce in a big crunch/big bang universe. The 
generation of curvature perturbations on the brane is also examined 
from a five- dimensional perspective. 

Finally, as an application of our methods, we develop a new colliding- 
brane solution of the Hofava-Witten model of heterotic M-theory. 



Contents 



1 Introduction 1 

2 The early universe 7 

2.1 The horizon problem 7 

2.1.1 Inflation 8 

2.1.2 A collapsing universe 11 

2.2 Classical cosmological perturbation theory 13 

2.3 Perturbations in inflation 16 

2.3.1 Massless scalar field in de Sitter 16 

2.3.2 Including metric perturbations 19 

2.4 Perturbations in a collapsing universe 22 

2.4.1 The story of C and $ 23 

2.4.2 Power spectra 25 

2.4.3 The relation between ( and <3> 27 

2.5 Frontiers 28 

2.5.1 Matching at the bounce 28 

2.5.2 Breaking the degeneracy: gravitational waves 29 

2.6 Summary 30 

3 Braneworld gravity 32 

3.1 The Randall- Sundrum model 32 

3.1.1 The Israel matching conditions 33 

3.2 Static solutions 34 

3.3 Cosmological solutions 35 

3.3.1 Generalising Birkhoff's theorem 36 

vii 



CONTENTS 



3.3.2 The modified Friedmann equations 39 

3.4 A big crunch/big bang spacetime 41 

3.4.1 Trajectories of empty branes 43 

3.4.2 Brane-static coordinates 43 

3.5 Four- dimensional effective theory 46 

3.5.1 Exact versus inexact truncations 46 

3.5.2 The moduli space approximation 47 

4 Conformal symmetry of braneworld effective actions 52 

4.1 Derivation of the effective action 52 

4.1.1 Conformal gauges 55 

4.1.2 Generalisation to other models 56 

4.1.3 Cosmological solutions 57 

4.2 AdS/CFT 59 

5 Solution of a braneworld big crunch/big bang cosmology 62 

5.1 Introduction 62 

5.2 Three solution methods 66 

5.3 Series expansion in time 70 

5.4 Expansion in Dirichlet and Neumann polynomials 72 

5.4.1 Background 72 

5.4.2 Perturbations 74 

5.5 Expansion about the scaling solution 76 

5.5.1 Scaling solution for the background 77 

5.5.2 Evolution of the brane scale factors 81 

5.5.3 Analytic continuation of the bulk geometry 83 

5.5.4 Higher-order corrections 88 

5.5.5 Treatment of the perturbations 91 

5.6 Comparison with the four-dimensional effective theory 95 

5.6.1 Background 97 

5.6.2 Perturbations 100 

5.7 Mixing of growing and decaying modes 103 

5.8 Summary 105 

viii 



CONTENTS 



6 Generating brane curvature 107 

6.1 Introduction 107 

6.1.1 £ in the four-dimensional effective theory 108 

6.1.2 C on the branes 108 

6.2 Inducing brane curvature from a five- dimensional perspective . . . 109 

7 Colliding branes in heterotic M-theory 113 

7.1 Heterotic M-theory 113 

7.2 A cosmological solution with colliding branes 115 

7.2.1 Equations of motion 115 

7.2.2 Initial conditions 116 

7.2.3 A conserved quantity 117 

7.2.4 The scaling solution 118 

7.3 Discussion 120 

8 Conclusions 122 
A Israel matching conditions 124 
B Five-dimensional longitudinal gauge 128 

C Detailed results 130 

C.I Polynomial expansion for background 130 

C.2 Polynomial expansion for perturbations 132 

C.2.1 All wavelengths 132 

C.2. 2 Long wavelengths 135 

C.3 Perturbations from expansion about the scaling solution 136 

C.4 Bulk geodesies 139 

References 149 



ix 



Chapter 1 
Introduction 



There is no excellent beauty that hath not 
some strangeness in the proportion. 

Francis Bacon 

One of the most striking implications of string theory and M-theory is that 
there are extra spatial dimensions whose size and shape determine the particle 
spectrum and couplings of the low energy world. If the extra dimensions are 
compact and of fixed size R, their existence results in a tower of Kaluza-Klein 
massive modes whose mass scale is set by R~ l . Unfortunately, this prediction is 
hard to test if the only energy scales accessible to experiment are much lower than 
R . At low energies, the massive modes decouple from the low energy effective 
theory and are, for all practical purposes, invisible. Therefore, we have no means 
of checking whether the four-dimensional effective theory observed to fit parti- 
cle physics experiments is actually the outcome of a simpler higher-dimensional 
theory. 

The one situation where the extra dimensions seem bound to reveal them- 
selves is in cosmology. At the big bang, the four-dimensional effective theory 
(Einstein gravity or its string-theoretic generalisation) breaks down, indicating 
that it must be replaced by an improved description. Already, there are sugges- 
tions of improved behaviour in higher-dimensional string theory and M-theory. If 
matter is localised on two branes bounding a higher- dimensional bulk, the matter 
density remains finite at a brane collision even though this moment is, from the 
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perspective of the four-dimensional effective theory, the big bang singularity [1- 
3]. Likewise, the equations of motion for fundamental strings are actually regular 
at the collision in string theory, in the relevant background solutions [4, 5]. 

We will adopt this scenario of colliding branes as our model of the big bang. 
Our focus, however, will not be the initial singularity itself, but rather, the dy- 
namics of higher- dimensional gravity as the universe emerges from a brane col- 
lision. The model we study - the Randall-Sundrum model [6] - is the simplest 
possible model of braneworld gravity, consisting of two empty Z2-branes (or orb- 
ifold planes) of opposite tension, separated by a five- dimensional bulk with a 
negative cosmological constant. We develop a solution method for the bulk ge- 
ometry, for both the background and cosmological perturbations, in the form of 
a perturbative expansion in (V/c) 2 , where V is the speed of the brane collision 
and c is the speed of light [7]. Our solution allows us to track the evolution of 
the background and cosmological perturbations from very early times right out 
to very late times, providing a benchmark against which the predictions of the 
four- dimensional effective theory can be tested. 

We will find that the four- dimensional effective theory is accurate in two 
limits: that of early times, for which the brane separation is significantly less 
than the anti-de Sitter (AdS) radius L; and that of late times, for which the 
brane separation is significantly greater than L. In the former limit, the brane 
tensions and the warping of the bulk become negligible, and a simple Kaluza-Klein 
description consisting of four-dimensional gravity and a scalar field applies (the 
gauge field zero mode having been eliminated by the Z 2 projections). In the latter 
limit, however, in which the branes are both widely separated and slowly- moving, 
the physics is qualitatively very different. Rather than being uniform across the 
extra dimension, the low energy gravitational zero modes are now localised around 
the positive-tension brane, as shown by Randall and Sundrum [8]. Nevertheless, 
the four-dimensional effective theory describing this limit is identical, consisting 
of Einstein gravity and a scalar field, the radion, parameterising the separation 
of the branes. 

Surprisingly, however, our five- dimensional solution reveals that in the tran- 
sition between these two limits - from Kaluza-Klein to Randall-Sundrum gravity 
- the four- dimensional effective theory fails at first nontrivial order in (V/c) 2 . In 
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effect, the separation of the branes at finite velocity serves to excite massive bulk 
modes, which curb the accuracy of the four-dimensional effective theory until 
their decay in the late-time asymptotic region. This process generates a striking 
signature impossible to forge within any local four-dimensional effective theory; 
namely, the mixing of four-dimensional cosmological perturbation modes between 
early and late times. 

This mode-mixing is conveniently described in four- dimensional longitudinal 
gauge, in which the sole physical degree of freedom associated with adiabatic 1 
scalar perturbations is encoded in the Newtonian potential $4. For sufficiently 
long wavelengths such that {kt^l <C 1, the four-dimensional effective theory pre- 
dicts that 



where £4 is four-dimensional conformal time, and A4 and B4 are constants pa- 
rameterising the amplitudes of the two perturbation modes. The first mode, with 
amplitude A4, represents a curvature perturbation on constant energy density or 
comoving spatial slices, while the second mode, with amplitude B4, corresponds 
to a local variation in the time elapsed since the big bang. (In an expanding 
universe, the curvature perturbation comes to dominate over the time-delay per- 
turbation at late times, and hence is often referred to as the 'growing' mode. In 
a collapsing universe, however, the two roles are reversed and the growing mode 
corresponds instead to the time-delay perturbation). 

Now, if the dynamics were truly governed by a local four-dimensional effective 
theory, the perturbation amplitudes A4 and B4 would be constants of the motion. 
From our five-dimensional solution, however, we can compute the actual asymp- 
totic behaviour of the four- dimensional effective Newtonian potential, at early 
and at late times, by evaluating the five- dimensional metric perturbations on the 
positive-tension brane 2 . This allows us to identify the four- dimensional effective 
mode amplitudes A4 and B4 in terms of the underlying five- dimensional mode 

i.e. perturbations which do not locally alter the matter history. 
2 Explicitly, the four-dimensional effective metric is related to the metric on the positive- 
tension brane via a conformal transformation, and so the anti-conformal part of the four- 
dimensional effective metric perturbation - namely, $4 - is equal to the anti-conformal part of 
the induced metric perturbation on the brane. 



. -D4 

$4 = A 4 - - T , 
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amplitudes, which are truly constant. We find that, while the four-dimensional 
effective theory prediction (1.1) does indeed hold in the limit of both early and 
late times, the four-dimensional effective mode amplitudes A4 and B4 are mixed 
in the transition from early to late times. (For example, if the system starts 
out purely in the time-delay mode at small £4, then one ends up in a mixture 
of both the time-delay and the curvature perturbation modes as £4 — > 00). This 
mixing first occurs at order (V/c) 2 , reflecting the fact that the four-dimensional 
effective description holds good at leading (zeroth) order. Equivalently, parame- 
terising the mixing by a matrix relating the dimensionless mode amplitudes A4 
and B^iL 2 /V 2 ) at early times to their counterparts at late times, we find that 
this matrix differs from the identity at order (V/c) 2 . 

Although this mixing is very small in the limit of a highly non-relativistic 
brane collision, it is nonetheless of great significance when considering the origin 
of the primordial density perturbations. Inflation provides one possible for ex- 
planation for the origin of these primordial fluctuations by postulating an early 
phase of quasi-de Sitter expansion, in which quantum fluctuations of the infla- 
ton field are exponentially stretched and amplified into a classical, scale-invariant 
pattern of large-scale curvature perturbations. 

In a scenario such as the present, however, in which the universe undergoes 
a big crunch to big bang transition described by the collision of two branes, an 
alternative mechanism is feasible. In the ekpyrotic model [2], and its cyclic exten- 
sion [1, 9], the observed scale-invariant perturbations first arise as ripples on the 
two branes, produced as they attract one another over cosmological time scales. 
These ripples are later imprinted on the hot radiation generated as the branes 
collide. So far, this process has only been described from a four-dimensional ef- 
fective point of view, in which the perturbations are generated during a pre-big 
bang phase in which the four- dimensional Einstein-frame scale factor is contract- 
ing [10-13]. 

A key difficulty faced by the ekpyrotic and cyclic scenarios is a scale-invariant 
spectrum of perturbations is generated only in the time-delay mode parameterised 
by B4, which is the growing mode in a contracting universe. Naively, this mode 
is then orthogonal to the growing mode in the present expanding phase of the 
universe, which is the curvature perturbation mode parameterised by A4. (In an 
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expanding universe the time-delay perturbation mode B^ decays rapidly to zero). 
In order to seed the formation of structure in the present universe, therefore, it 
is essential that some component of the growing mode time-delay perturbation 
in the collapsing phase be transmitted, with nonzero amplitude, to the growing 
mode curvature perturbation post-bang 1 . 

Clearly, our discovery that the four-dimensional perturbation modes mix in 
the full five-dimensional setup sheds new light on this problem. The most simple- 
minded resolution would be for the ekpyrotic mechanism to generate a scale- 
invariant spectrum in B^ long before the collision (when the branes are relatively 
far apart and the four- dimensional effective theory is still valid), with a piece of 
this subsequently being mixed into the curvature perturbation A± shortly before 
the collision. Yet the failure of the four- dimensional effective theory at order 
(V/c) 2 , coupled with the likely five- dimensional nature of any matching rule used 
to propagate perturbations across the singularity, suggests that instead a reformu- 
lation of the problem in terms of purely five-dimensional quantities is necessary. 

To this end, we define the curvature perturbation on the brane directly in 
five dimensions, showing how it differs from the curvature perturbation in the 
four-dimensional effective theory at order (V/c) 2 . We then proceed to analyse the 
generation of brane curvature through the action of an additional five- dimensional 
bulk stress (over and on top of the negative bulk cosmological constant) serving 
to pull the branes together. Under the assumption that this additional bulk stress 
is small, a quantitative estimate of the final brane curvature perturbation prior 
to the collision can be obtained, although further progress is dependent upon the 
introduction of a specific five- dimensional model for the additional bulk stress. 

Finally, it is worth bearing in mind that there is good hope of experimentally 
distinguishing inflation and the ekpyrotic mechanism over the coming decades 
through their very different predictions for the long-wavelength spectrum of pri- 
mordial gravitational waves [30]. 

In addition to the cosmological ramifications of the work in this thesis, an- 
other line of enquiry being actively pursued is the application of our methods to 
more sophisticated braneworlds stemming from fundamental theory, such as the 

1 This paradox is a recurring theme in the literature, see e.g. [14-29]. 
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Hofava-Witten model of heterotic M-theory. Our efforts to date form the basis 
of the concluding chapter of the present work. 

#** 

The plan of this thesis is as follows: we begin in Chapter § 2 with an in- 
troduction to the early universe, encompassing cosmological perturbation theory 
and the generation of the primordial density perturbations, via both inflation 
and the ekpyrotic mechanism. In Chapter § 3, we review the Randall-Sundrum 
model and braneworld cosmology, before proceeding to introduce the low energy 
four- dimensional effective theory. The subsequent chapters present original re- 
search material: in Chapter § 4, we investigate the role of conformal invariance 
in the braneworld construction, showing how the form of the effective action up 
to quadratic order in derivatives is fully constrained by this symmetry [31]. We 
also consider how to improve the effective action for a pair of branes of opposite 
tension through application of the AdS/CFT correspondence. Chapter §5 forms 
the heart of the present work. In this chapter, we present our solution methods 
for the bulk geometry, and apply them to calculate the behaviour of both the 
background and perturbations in a big crunch/big bang cosmology [7]. We high- 
light the failure of the four-dimensional effective theory, and compute explicitly 
the mixing of four-dimensional effective perturbation modes. In Chapter §6, we 
revisit the generation of curvature perturbations from a five-dimensional perspec- 
tive, considering the action of an additional five-dimensional bulk stress. Finally, 
in Chapter § 7, we present work in progress seeking a cosmological solution of 
the Hofava-Witten model with colliding branes, in which the five-dimensional 
geometry about the collision is that of a compactified Milne spacetime, and the 
Calabi-Yau volume at the collision is finite and nonzero. 



6 



Chapter 2 

The early universe 



Now entertain conjecture of a time 

When creeping murmur and the poring dark 

Fills the wide vessel of the universe. 



Henry V, Act IV. 



2.1 The horizon problem 

Observations of the cosmic microwave background (CMB) allow a precise deter- 
mination of the nature of the primordial density perturbations. To date, these 
observations require the primordial density perturbations to be small-amplitude, 
adiabatic, Gaussian random fluctuations with a nearly scale-invariant spectrum 
[32]. Even though this is almost the simplest conceivable possibility, to generate 
density perturbations of this nature demands new physics beyond the Standard 
Model. 

The essence of the problem is simple causality: a universe originating in a big 
bang has a particle horizon associated with the fact that light has only travelled 
a finite distance in the finite time elapsed since the big bang. For a universe with 
scale factor a given in terms of the proper time t by a = (t/t ) p , the size of this 
horizon in comoving coordinates is 




(2.1) 
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where the comoving Hubble radius = (da/dt) -1 = p~ l tf ) t 1 ~ p and tj is some 
initial time. For standard matter, p < 1 (e.g. p = 1/2 for radiation, p = 2/3 
for matter), and so the size of the horizon is always increasing. At sufficiently 
late times, the integral is dominated by its upper limit; light travels the greatest 
distance at late times, and the horizon grows as d ~ IK -1 ~ t l ~ p . 

Observations of the CMB, however, indicate that the Universe was quasi- 
homogeneous at the time of last scattering on scales much larger than the size 
of the causal horizon at that time. Although the angular scale subtended by the 
horizon at last scattering is only ~ 1°, the temperature of the microwave sky is 
uniform to one part in a hundred thousand over much larger angular scales. This 
puzzle of explaining why the universe is nearly homogeneous over regions a priori 
causally independent is known as the horizon problem. 

2.1.1 Inflation 

One possible resolution of the horizon problem is that the early universe un- 
derwent a period of inflation, defined as an epoch in which the scale factor is 
accelerating, i.e., a > 0, where dots denote differentiation with respect to proper 
time. The comoving Hubble radius !K _1 is therefore shrinking during inflation, 
since "K = a > 0. (Note, however, that the converse, shrinking comoving Hubble 
radius implies inflation, is not true, as we will see in the next section). 

To see how inflation solves the horizon problem, let us adopt a simple model 
in which the scale factor expands exponentially, a = exp(H(t — to))- The size 
of the comoving horizon, equal to the comoving distance travelled by a light ray 
during inflation, is given by 

/** c\t' 

<W*) = / ™=0Cr 1 -Xr 1 , (2.2) 

where the comoving Hubble radius JC -1 = (aH)~ l . Thus, due to the exponential 
expansion of the scale factor, after a sufficient number of e-foldings (defined as 
iV = ln(a/aj)), the integral is dominated by the lower limit; light rays travel 
the greatest distance at early times and we find a near-constant horizon size 
(iinfl. ~ IJQ" 1 . By choosing the time of onset of inflation tj to be sufficiently small, 
and the proper Hubble parameter H to be sufficiently large, we can always ensure 
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that the horizon size at the end of inflation is sufficiently large to solve the horizon 
problem. 

As a rough estimate of the number of e-folds of inflation required, we stipulate 
that the comoving distance travelled by light during the radiation era must equal 
the comoving distance travelled by light during inflation, i.e., the causal horizon 
at the end of inflation equals the size of our past light-cone at that time, since we 
live approximately at the end of the radiation era. This ensures that the portion 
of the surface of reheating visible to us (and hence any later surface, such as the 
surface of last scattering) is causally connected. 

From (2.1) with p = 1/2, we see that the comoving distance travelled by 
light during the radiation era is approximately equal to the comoving Hubble 
radius at matter-radiation equality, "K^ow From (2.2), this must then equal the 
comoving Hubble radius at the onset of inflation, 'K^ 1 . The shrinking of the 
comoving Hubble radius during inflation is therefore matched by the growth of 
the comoving Hubble radius during the radiation era. Since during inflation 
"K" 1 ~ a -1 , whereas during radiation domination JC _1 ~ a, 



-1 ~~ or- 1 



(2.3) 



where a rn is the scale factor evaluated at reheating, etc. Finally, since the 
scale factor is inversely proportional to the temperature during radiation dom- 
ination, we deduce that the number of e-folds of inflation required is roughly 
N ~ ln(T rn /T now ) « ln(10 12 GeV/1 meV) ~ 55, assuming that reheating occurred 
approximately at the GUT scale 1 . 

Having dealt with the horizon problem, we can now ask about the nature of 
the matter required to sustain a period of inflation. From the Friedmann equation 
for a homogeneous and isotropic universe, 

a 1 , , 
- = --(„ + 3P) (2.4) 

(where we have set 8irG = 1), we see that inflation (a > 0) requires matter with 
density p and pressure P such that p + 3P < 0. This corresponds to an exotic 



1 Note that if inflation does terminate at the GUT scale, observational upper bounds on the 
density of monopoles constrain the GUT scale to be less than 10 12 GcV [33]. 
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equation of state with uj = Pj p < —1/3. Ruling out a cosmological constant with 
uj = — 1 on the grounds that inflation must at some point terminate, the next 
simplest candidate is a massless scalar field <p, minimally coupled to gravity, and 
with a potential V(tp). The stress tensor then takes the form of that for a perfect 
fluid, with p = tp 2 /2 + V and P = <p 2 /2 — V. Inflation therefore requires that the 
potential energy is greater than the scalar field kinetic energy, up 2 < V. 
From the action 

S=\J ^[R-{dipf-2V{ip)], (2.5) 
the equations of motion are 

3H 2 = ~V 2 + V, (2.6) 
= <p + 3H<p + V )ip , (2.7) 
H = -\^\ (2-8) 

where the Hubble parameter H = a /a, and the last equation follows from the 
first two. A particularly useful limit in which the dynamics are simplified is the 
slow-roll regime, in which the scalar field kinetic energy is negligible, tp 2 <C V, 
and the system is over-damped, ip Hip. This holds provided the slow-roll 
parameters 

e = i(liy«i£ (29) 

V = + (2.10) 

' V Hip 2 H 2 ' v ; 

are much smaller than unity (the approximate relations following in the limit 
when this is true) . Under conditions of slow roll, the number of e-foldings obtained 
during inflation is simply 

N= / —H(ip) w / ^-dtp. (2.11) 

Finally, let us introduce a simple model that will prove useful, obtained by 
setting the equation of state parameter uj = P/p to a constant. In this case 
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the dynamics are exactly solvable (independently of considerations of slow roll), 
yielding the scaling solution 




(2.12) 



where to and Vq are given in terms of the constant parameter phy t G = p — 1 and 
Vq = p (1 — 3p)/to- The parameter p is in turn related to the equation of state 
parameter to by p = 2/3(1 + to). We recover slow-roll inflation in the limit where 
p>l(w~-l), since the slow-roll parameters are e — l/p and rj — 2/p . In this 
case the potential takes the form of a positive-valued exponential. 



2.1.2 A collapsing universe 

An alternative resolution of the horizon problem is that the big bang is not the 
beginning of the universe. Instead, the present expanding phase of the universe 
would be preceded by a collapsing phase, with a concomitant transition from big 
crunch to big bang. 

A concrete example is provided by the simple model with constant to discussed 
above. From (2.12) with p > 1, we obtain an expanding universe for times 
< t < oo. If instead we take < p < 1, we obtain a contracting universe with 
to < and the time coordinate taking values in the range — oo < t < 0. The 
potential now corresponds to a negative- valued exponential. Since u > —1/3, 
the universe is not inflating and a < 0. Nonetheless, the comoving Hubble radius 
is still shrinking during the collapse: re-expressed in terms in of conformal time 
t = — [t/(p — l)] 1_p , the scale factor and scalar field in (2.12) are 



^= V^lnM, (2.13) 
1 — p 



where — oo < r < 0. (Note that for < p < 1, the range — oo < r < is 
mapped to — oo < t < 0, whereas for p > 1, the same range is instead mapped 
to < t < oo). For general p, the magnitude of the comoving Hubble radius is 
therefore 

|:k|- 1 = |(p~ 1 -iH, ( 2 - 14 ) 
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where we have taken the absolute value since in a collapsing universe 3~C < 0. With 
this definition, a shrinking comoving Hubble radius corresponds, in a collapsing 
universe, to the condition !K = a < 0. In contrast, in an expanding universe, a 
shrinking Hubble radius implies "K = a > 0, and hence inflation. 

For a universe undergoing a bounce, the resolution of the horizon problem is 
trivial. The comoving distance travelled by light rays during the collapsing phase 
is as given in (2.1), except that now, owing to the shrinking of the comoving 
Hubble radius during the collapse, the integral is dominated by its lower limit: 
provided the collapsing phase began at some sufficiently negative initial time ti, 
the comoving distance travelled scales as By increasing the duration of 

the collapsing phase, the horizon size at the big bang can be made arbitrarily 
large. A lower bound is given by the size of our past light cone at the big bang, 
yielding the estimate | HHT^ | 1 ~ IK^ow' ^- e -> ^ ne shrinkage of the comoving Hubble 
radius during the collapse is matched by its growth in the subsequent expanding 
phase. 

During both inflation and a collapsing era then, the comoving Hubble radius 
shrinks inside fixed comoving scales. The two scenarios differ, however, when 
interpreted in physical coordinates. While the proper distance corresponding 
to a fixed comoving length scales as a = (t/to) p , the proper Hubble radius is 
H^ 1 = t/p. In an inflating universe (p > 1) therefore, the physical wavelength 
corresponding to a fixed comoving wavevector is stretched outside the Hubble 
radius as t — > oo. In a collapsing universe (0 < p < 1), however, the opposite 
happens: as t approaches zero from negative values, the shrinking of the Hubble 
radius outpaces the shrinking of physical wavelengths. 

Finally, for completeness, we note that a stability analysis [12] shows that the 
background solution (2.12) is a stable attractor under small perturbations, both 
in the expanding case where p > 1 (u < — 1/3), and also in the contracting case 
provided that p<l/3(o;>l). This conclusion remains valid when we include 
more general forms of matter. Intuitively, in an inflating universe with uj ~ — 1, 
this is because the scalar field potential energy density is nearly constant, whereas 
the curvature (ex a~ 2 ), matter (oc a~ 3 ), radiation (oc a~ A ) and other forms of 
energy density all decrease as the scale factor grows. Likewise, in a collapsing 
universe with u 3> 1, the scalar field kinetic energy density increases as a~~ 3 ( 1+Ui ^ as 
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the scale factor shrinks, whereas curvature, matter, radiation and other forms of 
energy density increase at a slower rate. Consequently, in each case the dynamics 
are largely insensitive to the initial conditions: after a few e-folds of expansion 
or contraction, the solution converges rapidly to the attractor. This resolves 
the so-called flatness and homogeneity problems, explaining why the present-day 
universe is so uniform and geometrically flat, without recourse to fine-tuning of 
the initial conditions. 



2.2 Classical cosmological perturbation theory 

We now turn our attention to the behaviour of classical perturbations in linearised 
general relativity. Considering only the scalar sector relevant to the density per- 
turbations, the perturbed line element for a flat FRW universe can be expressed 
in the general form 

ds 2 = a 2 (r) (-(1 + 20) dr 2 + 2B :i drdx i + ((1 - 2^)% + 2E :ij ) daW) , (2.15) 

where the spacetime functions 0, ip, B and E are linear metric perturbations. 
Under a gauge transformation £ M = <9 l £), the perturbations transform as 

<P -> (2.16) 

^ - ^ + ^°, (2.17) 

B -> £ + (2.18) 

E -> E-Z, (2.19) 

where the primes indicate differentiation with respect to the conformal time r. 
We may then construct the gauge-invariant Bardeen potentials [34] 

$ = + <K(B - E') + (B - E')', (2.20) 
# = ^ - <K(B - E 1 ). (2.21) 

Another very useful gauge-invariant quantity is the curvature perturbation on 
comoving hypersurfaces, (, defined in any gauge by 

C = ^-^p? (2-22) 
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where the potential q satisfies q ti = 8T®. (Thus, in comoving gauges characterised 
by 8T® = q ti = 0, k 2 ( is proportional to the curvature of spatial slices, given by 
i?(3) = -4k 2 ^p). 

We can easily relate this definition to an equivalent one in terms of the Bardeen 
potentials as follows. Working in longitudinal gauge for convenience (Bl = El = 
0, and hence </>l = $ and ipL = where the subscript L will be used to denote 
quantities evaluated in this gauge), the perturbed G° Einstein equation yields the 
momentum constraint equation 

(<5G°) L = ~ [V + X*] ti = W)x = QL,i (2-23) 

(where 87rG = 1 still). Making use of the background Friedmann equation in the 
form %' — "K 2 = —a 2 (p + P)/2, we therefore find (in any gauge now) 

C = ^- n ^^ ] = g-| + (2.24) 

where we have reverted to proper time in the last equality. 

The remaining components of the perturbed Einstein tensor, in longitudinal 
gauge, are: 

(5G° ) L = — [3IK 2 $ + 3M' - V 2 ^] , (2.25) 

(5G\) L = — [2f + 4M'- V 2 (*-$) +2M' + 2(2%' + "K 2 )^] 8) 

+\s ik (*-$) Jk . (2.26) 
er 

Specialising to the case of a single scalar field, the absence of anisotropic stresses 
(8Tj = for % 7^ j) forces the two Bardeen potentials to coincide, thus $ = 

Evaluating the perturbed stress tensor to linear order and making use of the 
background equations of motion in conformal time, 

= 3J{ 2 -^' 2 -aV, (2.27) 
= ^ + 2%' + aV„ (2.28) 
= %' - % 2 + ^' 2 , (2.29) 
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we find the linearised Einstein equations in longitudinal gauge read 



V 2 $ - 3M$' - (IK' + 2:K 2 )$ 



1 1 



(2.30) 



(2.31) 



$" + 3IK$' + (:K' + 2:K 2 )$ 



(2.32) 



where 6(fL denotes the perturbation of the scalar field in longitudinal gauge. 

We now subtract (2.30) from (2.32). Using the background equation of motion 
(2.28), we can further subtract (4J{ + 2(p'q/<p' ) times (2.31) to find the following 
second-order differential equation for the Newtonian potential: 



(Note the absence of ensures that this equation is now gauge invariant). 
Combining (2.33) with (2.24) (recalling $ = W), we then have 



This tells us that on large scales, where the gradient term is negligible, ( is 
conserved (except possibly if <p' = at some point during the cosmological evolu- 
tion, e.g. during reheating when the inflaton field undergoes small oscillations). 
In fact, even if <p' Q does vanish at a particular moment in time, it turns out that 
( is still conserved outside the Hubble radius. For quite general matter, it can 
be shown (see e.g. [35]) using the perturbed continuity equation that ( is con- 
served on super-Hubble scales, provided only that the perturbations are adiabatic 
{i.e., every point in spacetime goes through the same matter history. Note this 
automatically holds for single field models). 

This accounts for the utility of ( in the calculation of inflationary perturba- 
tions: once a mode exits the Hubble radius, ( for that mode is conserved. Despite 
our ignorance of the processes occurring during reheating, the initial conditions 
for a mode re-entering the Hubble radius in the radiation era are fully specified 
by the conserved value of £ calculated as the mode left the Hubble radius during 
inflation. 




(2.33) 




(2.34) 
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2.3 Perturbations in inflation 

Inflation stretches initial vacuum quantum fluctuations into macroscopic cosmo- 
logical perturbations, thereby providing the initial conditions for the subsequent 
classical radiation-dominated era. Remarkably, the nearly scale-invariant spec- 
trum of perturbations predicted by inflation exactly accounts for present-day 
observations of the CMB. The computation of the primordial fluctuations arising 
in inflationary models was first discussed in [36-40]; the present section com- 
prises a review of these classic calculations. (Other useful review articles include 
[35, 41-43]). 

2.3.1 Massless scalar field in de Sitter 

Before launching into a full calculation of the quantisation of a scalar field in- 
cluding gravitational backreaction, let us first sharpen our claws on a simple toy 
model; namely, a massless scalar field on a fixed de Sitter background. The cos- 
mological scale factor is given by a oc exp(Ht), which in conformal time reads 
a = —1/Ht. (Note the conformal time takes values in the range — oo < r < 0). 
The action for a massless scalar field in this geometry is 

S=-\j d 4 oV^(^) 2 = \J drd 3 xaV 2 - (V^) 2 ], (2.35) 

where differentiation with respect to conformal time is denoted by a prime. 
Changing variables to x — at P an d integrating by parts, we obtain 

S= l ~ |drd 3 x[ X /2 -(Vx) 2 + ^X 2 ], (2.36) 

where the kinetic term is now canonically normalised. The fact that the scalar 
field lives in de Sitter spacetime rather than Minkowski results in a time- dependent 
effective mass 

m 2 cS = -- = -- 2 . (2.37) 
The quantum field x may then be expanded in a basis of plane waves as 

X(r,x) = JL J d'k[^x k (T)e^ + alxl(r)e-^% (2.38) 



16 



2.3 Perturbations in inflation 



where the creation and annihilation operators a' and a satisfy the usual bosonic 
commutation relations 

[%, %'] = 4'] = °' fa> «y = - ( 2 - 39 ) 

The function Xfc( r ) is a complex time-dependent function satisfying the classical 
equation of motion 

xl + (k 2 ~ ~)Xk = 0, (2.40) 
a 

representing a simple harmonic oscillator with time- dependent mass. In the case 
of de Sitter spacetime, the general solution is given by 

X k = Ae~^(l-^)+Be^(l + ^). (2.41) 

Canonical quantisation consists in imposing the following commutation relations 
on constant-r hypersurfaces: 

[ X {r,x), x(t,x')} = [7t x {r,x), n x (r,x')} = (2.42) 

and 

[X(r,x), n x (r,x')}=ih8(x-x'), (2.43) 

where the canonical momentum n x = 5S/5x = x' ■ Substituting the mode expan- 
sion (2.38) into the commutator (2.43), and making use of (2.39), we obtain the 
condition 

XkXt ~ Xlx'k = ift, (2-44) 

which fixes the normalisation of the Wronskian. 

Specifying a particular choice for Xfc( r ) corresponds to the choice of a partic- 
ular physical vacuum |0), defined by a^|0) = 0. A different choice for Xk{ T ) leads 
to a different decomposition into creation and annihilation operators, and thus 
to a different vacuum. In our case, the most natural physical prescription for the 
vacuum is to take the particular solution that corresponds to the usual Minkowski 
vacuum Xk ~ exp(— ikr) in the limit k\r\ ^> 1. This is because the comoving 
Hubble radius |t| _1 is shrinking during inflation: provided one travels sufficiently 
far back in time, then any given mode can always be found within the Hubble 
radius (i.e. we can always choose \t\ sufficiently large such that k\r\ ^> 1). Once 
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the wavelength of the mode lies within the Hubble radius, the curvature of space- 
time has a negligible effect (i.e. the effective mass term in (2.40) is negligible) 
and the mode behaves as if it were in Minkowski spacetime. Taking account of 
the normalisation condition (2.44), we thus arrive at the Bunch-Davies vacuum 

» = VS^ 1 " f>- (2 ' 45) 

The power spectrum is defined via the Fourier transform of the correlation 
function, 

(O|0(a?i)0(i? 2 )|O) = J d 3 fc e ^-*)^S, (2.46) 

where 

-2 ? -3m _ i,„ i2 \Xk\ 



2^*"% = Irf = (2.47) 
With this definition, the variance of <p in real space is 



(v(x) 2 ) = 9 v d\nk, (2.48) 



<p 

and so a constant IP^ corresponds to a scale- invariant spectrum, i.e., one in which 
the field variance receives an equal contribution from each logarithmic interval of 
wavenumbers. 

Thus we find that modes well outside the Hubble radius, i.e., those for which 
k\r\ <C 1, possess the scale-invariant spectrum 

* (|Q • ( 2 - 49 ) 

(Note also that in the opposite limit k\r\ ^> 1, in which the mode is well inside 
the Hubble radius, we recover the usual result for fluctuations in the Minkowski 
vacuum, = h(k/2ira) 2 ). 

Looking back, we can trace the origin of this scale invariance to the — 2/r 2 
effective mass term in the classical equation of motion (2.40). Thanks to this 
term, once a given mode \k exits the Hubble radius, it ceases to oscillate (i.e., 
the mode 'freezes') and instead scales as r _1 in the limit as r — > 0. Since this 
amplification sets in at a time r ~ k -1 , the late-time mode amplitude picks up an 
extra factor of k~ l on top of the k~ 1 ^ 2 factor from the initial Minkowski vacuum, 
leading to scale invariance, i.e., \xk\ 2 k~ 3 . 
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2.3.2 Including metric perturbations 

We now turn to the full calculation of the spectrum of inflationary perturbations 
including gravitational backreaction. 

Recalling our discussion of classical cosmological perturbation theory, there 
are a total of five scalar modes; four from the metric ((f), ip, B and E), and the 
perturbation of the inflaton itself, Sip. Gauge invariance allows us to remove two of 
the five functions by fixing the time and spatial reparameterisations x M — > 
where £ M = <9*£), while the Einstein constraint equations remove a further two. 
We are left with a sole physical degree of freedom. 

In order to derive the perturbed action up to second order in the metric pertur- 
bations, we will work in the ADM formalism, following [44]. Foliating spacetime 
into spacelike 3-surfaces of constant time, the metric can be decomposed in terms 
of the lapse N(t,x t ), the shift N l (t,x l ) and the spatial 3-metric hij, as 

ds 2 = -iV 2 dt 2 + hijidx* + N i dt)(dx j + N j dt). (2.50) 

Spatial indices are lowered and raised using the spatial 3-metric. The 4-curvature 
can be decomposed in terms of the spatial curvature and the extrinsic cur- 
vature tensor 

Kij = ^ (kj - ViNj - VjNi) (2.51) 
according to the Gauss-Codacci relation [45] 

R = R {3) + KijK ij - K 2 , (2.52) 

where K = K\ and we have omitted a total derivative term. The relevant action 
(2.5) is thus 

S = i J dtd^xVh NR {3) - 2NV + N- 1 (E ij E i: > - E 2 ) 

+N-\<p - iVV,*) 2 - Nh ij <p ti cpj ] , (2.53) 

where, for convenience, we have defined Eij = NKy. 

The ADM formalism is constructed so that (p and play the role of dynam- 
ical variables, while iV and iV* are simply Lagrange multipliers. Our strategy, 
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after fixing the gauge, will be to solve for these Lagrange multipliers and then 
backsubstitute them into the action. 

Considering only scalar modes, a convenient gauge is 

5<p = 0, hij = a 2 (r) [(1 - 2C)<y , (2.54) 

where a(r) is the background scale factor but ( is a first-order perturbation. 
In this gauge the scalar field is unperturbed, hence 5T® = and the gauge is 
comoving. (Note also that the gauge is fully specified by the above). 

The equations of motion for N l and N are the momentum and Hamiltonian 
constraints respectively, which in this gauge take the form 

V i [N- 1 (E i j -S i j E)] = 0, (2.55) 
ij(3) -2V- N- 2 (E tj E ij - E 2 ) - N~ 2 (pl = 0. (2.56) 

At linear order, the solution of these equations is 

where d 2 is a shorthand notation for <5 lJ didj. 

Substituting the above back into the action (2.53), and expanding out to 
second order, we obtain a quadratic action for £. Note that for this purpose 
it is unnecessary to compute N or N l to second order. This is because these 
second order terms can only appear multiplying the zeroth order Hamiltonian 
and momentum constraint equations dL/dN and dL/dN l , which vanish for the 
background. Direct replacement in the action gives, up to second order, 

S = - [ dtd 3 xae-<(l-^-)[4d 2 (-2(d i () 2 -2a 2 Ve- 2 <} 
2 J H 

We-^(l - |)^[-6(C - H) 2 + cp 2 }, (2.58) 

where we have neglected a total derivative and (di() 2 = 5 13 di(dj(. After inte- 
grating by parts, and using the background equations of motion (2.6) to (2.8), 
we obtain 

S = \J dtd 3 x|| [a 3 C 2 - a(d t C) 2 ] = -\J dW=sf| GTfyC&C), (2-59) 
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where is the background FRW metric. 

To connect with our earlier discussion of a scalar field in de Sitter space, let us 
pass to conformal time and canonically normalise the kinetic term by introducing 
the variable u = z(, where z = a(p' Q /'H. We find 

S = \ f drd 3 x \v' 2 - (diuf + — (2.60) 
Z J z 

the action for a scalar field in Minkowski spacetime with a time-dependent mass 
z" jz. (Recall that in de Sitter spacetime the effective mass was instead a" /a). 
Each Fourier mode obeys the equation of motion 

< + (k 2 - -)v k = 0, (2.61) 



or equivalently, 

C£ + 2(-)Cfc + fc 2 G = 0. (2-62) 

z 

(Since ( is gauge- invariant, and z depends only on background quantities, this 
relation in fact holds in any gauge). 

Under conditions of slow-roll, the evolution of </?o and H is much slower than 
that of the scale factor a, and so to leading order in slow-roll one finds 



z" a" 



(2.63) 

z a 

Consequently, all the results of the previous section now apply to our variable v in 
the slow-roll approximation. The correctly normalised expression corresponding 
to the Bunch-Davies vacuum is approximately 

»^{f k e ~ lkT ( 1 -^- (2 ' 64) 
In the super-Hubble limit, k\r\ 1, this becomes 

^ „ * « (2.65) 

V 2k kr V 2k k K J 

where we have used a ~ —1/Ht and the asterisk indicates that, to get the most 
accurate approximation, we should evaluate the value of H at the moment the 
mode crosses the Hubble radius. 
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Hence, on scales much larger than the Hubble radius, the power spectrum for 
C is given by 

k 3 k| 2 ft (Hi 



where we have reverted once more to physical time, and the quantities labelled 
with asterisks are to be evaluated when the mode exits the Hubble radius. This 
is the celebrated result for the spectrum of scalar cosmological perturbations 
generated from vacuum fluctuations during slow-roll inflation. 

The spectrum is not strictly scale-invariant, however, since the values of the 
Hubble parameter and the scalar field velocity evolve slowly over time. This 
introduces a small additional momentum dependence conveniently parameterised 
by the spectral index n s , where 



(2.67) 

using the slow-roll parameters as defined in (2.9) and (2.10). A spectral index 
n s = 1 therefore corresponds to an exactly scale-invariant spectrum, whereas n s < 
1 and n s > 1 correspond to a red and a blue spectrum respectively. ( Equivalent ly, 
we have k 3 \( k \ 2 oc /c™ s_1 ). 

A useful heuristic to understand the parametric dependence of the power 
spectrum (2.66) is as follows. Consider a mode just about to cross the Hubble 
radius and freeze out. On a constant-time slice the quantum fluctuations in the 
inflaton field will be of order 5(f ~ ii?*, approximating the background as de Sitter 
space and using (2.49). Comoving surfaces, characterised by 5ip = 0, are therefore 
offset from surfaces of constant time by a time delay 5t ~ Stp/tp ~ H*/<p. The 
exponential expansion of the background spacetime as a ~ exp(if„,t) then leads 
to a comoving curvature perturbation £ ~ 5a /a ~ H*5t ~ H%/<fi, consistent with 
(2.66). 



2.4 Perturbations in a collapsing universe 

We have seen that models formulated in terms of gravity and a single scalar 
field possess only one physical scalar degree of freedom at linear order. In the 
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case of inflation, we chose to parameterise this degree of freedom in a gauge- 
invariant fashion through use of the comoving curvature perturbation (. The 
utility of this variable lies in its conservation on super-Hubble scales, allowing 
the amplitudes of modes re-entering the Hubble radius during the radiation era 
to be calculated independently of the detailed microphysics of re-heating. In the 
case of a collapsing universe, however, the situation is more subtle. The final 
spectrum of perturbations in the expanding phase is, in general, sensitive to the 
prescription with which the perturbations are matched across the bounce. In 
particular, it does not necessarily follow that ( is conserved across the bounce. 
In this section we will therefore adopt a more generic approach, computing the 
behaviour of both ( and the Newtonian potential $ during the collapse. 

2.4.1 The story of ( and <3> 

Let us start by collecting together a number of our earlier results concerning £ 
and $. These variables satisfy the second order linear differential equations given 
in (2.33) and (2.62), and are related to each other via (2.24) (recalling that the 
absence of anisotropic stress for scalar field matter implies that $ = \l/) and 
(2.34). To express these relations more compactly, it is useful to introduce the 
surrogate variables u and z/, defined as 



where the background quantity z = a(f' /!K as before. In particular, u and v 
have the same /c-dependence as $ and £, and hence the same spectral properties. 




v = zC, 



(2.68) 



Defining 9 = l/z, (2.33) and (2.62) are 













(2.69) 



z 



and the relations (2.24) and (2.34) become 





(2.72) 
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Choosing a boundary condition for u and v corresponds to specifying a vacuum 
state for the fluctuations. As usual, we will take this to be the Minkowski vacuum 
of a co moving observer in the far past, at a time when all comoving scales were 
far inside the Hubble radius. Thus, in the limit as r — > — oo, 



■3/2 -ikr 



u -> iVh(2k)- 
v -> Vh(2k)- l/2 e- lkT . 



[2.73) 
[2.74) 



(It is easy to check from (2.71) that these two boundary conditions are equivalent). 

Assuming the equation of state is constant while all wavelengths of cosmo- 
logical interest leave the Hubble radius, from our scaling solution (2.12) with 
uj = 2/3p — 1 , we find 



V 



z" 
z 



(p- l) 2 r 2 ' 
p(2p-l) 

( p _l)2 r 2- 



(2.75) 
(2.76) 



The solutions of (2.69) and (2.70) are then easily obtained in terms of Hankel 
functions: 



u(x) 
v(x) 



V2 [ A ^H^(x) + A^H^( 

bWh£\x) + b®h 



1/2 



X) 



(2.77) 
(2.78) 



where x = k\r\ is a dimensionless time variable, A^ 1 ' 2 ) and B^ 1 ' 2 ) are constants, 



and the order of the Hankel functions if. 



(1,2), 



X) IS 



a = [(6"/6)t 2 + 1/4] 1/2 
/5 = [(z'V^^ + l/^ 1 / 2 



1+p 
1 — p 
l-3p 



1 — p 



(2.79) 
(2.80) 



In the far past (x — > oo), when comoving scales are well inside the Hubble 
radius, the asymptotic behaviour of the Hankel functions is 



x) 




2 

— exp 



S7T 7T 
±1 \ X 

2 4 



(2.8i; 
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(where the + sign corresponds to H, 
(2.74) then imply 



(x)). The boundary conditions (2.73) and 



it- 



is 



\ 2 (hirx/4k) 1/2 H { p\x), 



(2.82) 
(2.83) 



where Ai and A2 are the k- independent complex phase factors 



exp[i(2a + 3)vr/4], 
exp[i(2/? + 1)tt/4]. 



(2.84) 
(2.85) 



Pausing to catch our breath, we notice a curious duality: the equation of 
motion for u ((2.69), along with (2.75)) is invariant under p —>■ 1/p [11]. More- 
over, the boundary condition, (2.73), is independent of p (as is natural, since 
the boundary condition is imposed in the far past when comoving scales are well 
inside the Hubble radius). Consequently, our expressions (2.79) for a, and (2.82) 
for u, are invariant under p — > 1/p. The spectrum of fluctuations of the New- 
tonian potential $ in an expanding universe with p > 1 is therefore identical to 
that in a collapsing universe with p — 1/p < 1. (The same cannot be said for £, 
however, as (2.76) is not invariant under p — > 1/p). 

2.4.2 Power spectra 

Returning now to the calculation of the power spectra for £ and $, at late times 
when comoving scales are well outside Hubble radius (x — > 0), the asymptotic 
behaviour of the Hankel functions is 



where s > and T(s) is the Euler gamma function. On scales much larger than 
the Hubble radius, the power spectra for ( is therefore 




(2.86) 




(2.87) 
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where asterisked quantities are to be evaluated at Hubble radius crossing (k = 
aH), and the /^-independent numerical factor A(p) is 



We immediately see that the power spectrum for ( is only scale-invariant when 
(3 = 3/2, cancelling the x-dependence on the right-hand side of (2.87). Equiva- 
lently, the spectral index (given by — 1 = 3 — 2/3) must be unity for scale 
invariance. From (2.80), this requires p — > oo (u rs —1), corresponding to an in- 
flating universe in the slow-roll limit (recall that the slow-roll parameter e = 1/p ). 
In this limit, the numerical factor A(p) tends to unity and we recover our earlier 
result (2.66). 

On the other hand, the power spectra on super-Hubble scales for the Newto- 
nian potential $ is given by 



The spectral index for $ is thus n$ — 1 = 1 — 2a, with scale invariance requiring 
a = 1/2. From (2.79), this is satisfied in only two cases: firstly, that of an expand- 
ing universe undergoing slow- roll inflation (p — > oo); and secondly, the case of a 
very slowly collapsing universe with p -C 1 {ui ^> 1). (In this latter scenario, the 
scalar field is rapidly rolling down a steep negative exponential potential). From 
our earlier remarks about duality, this result is exactly as expected: the power 
spectrum (2.89) is invariant under p — > 1/p, since a is invariant. Nonetheless, 
it is remarkable that a nearly scale-invariant spectrum can be obtained, without 
inflation, on a background which is very nearly static Minkowski. 

To sum up, while in the case of slow-roll inflation both $ and ( acquire a 
scale-invariant spectrum well outside the Hubble radius, in the case of a slowly 
contracting universe, $ is scale-invariant but ( is blue (n^ — 1 ~ 2, indicating 
greater power at short wavelengths). 

Our analysis has been restricted to the case in which the equation of state is 
constant. More generally, however, when the equation of state is a slowly varying 
function of time, one can introduce the 'fast-roll' parameters e and f] defined by 



Hp) = r 2 (/5) 



4?(1 -pf 
2irp 2 



(2.88) 






7/ 



1 




(2.90) 
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Analogously to the case of inflation, one can then express the spectral index in 
terms of these fast-roll parameters. A short calculation [12] yields the result 

n s -l = -4(e + 77). (2.91) 
2.4.3 The relation between £ and $ 

Since a generic scalar fluctuation can be described using either ( or $ with equal 
validity, the fact that one variable possesses a scale-invariant spectrum, while the 
other does not, is at first sight puzzling. To understand this behaviour, it is useful 
to probe further the relationship between ( and This is given by 

'-faff). ^ 

which follows directly from the definition of £ given in (2.24), and the equation 
of motion H 2 / H = —p. 

Expanding once more the Hankel function H^(x) as x — > 0, but this time 
retaining the subleading term as well, we find 

= -ir( a ) (£f - l -e— T{-8) (|) S + 0(x 2 - x 2+ °). (2.93) 

Up to fc-independent numerical coefficients, the leading and subleading behaviour 
of $ as t — > is therefore 

$ „ k~ 3 / 2 ~P r 1 -* + k~ 1/2+p f , (2.94) 

where we have used (2.82) and expanded all exponents up to linear order in p 
(implying r ~ t 1_p ). The leading term in this expression, the 'growing' mode, 
is responsible for the scale invariance of $. Since a/H ~ t 1+p , however, we 
immediately see from (2.92) that this term makes no contribution to £■ Instead, 
the leading contribution to ( comes from the subleading (or 'decaying') mode in 
$, which scales as k~ l / 2+p . The spectrum for £ is therefore boosted by a factor 
of k 2+2p relative to scale invariance, yielding a blue spectrum with — 1 w 2. 

For completeness, from (2.83), the growing and decaying modes for (, in the 
limit as t — > 0, are 

£ ~ /c" 1/2+p tg + /c 1/2_p f 1 -* t 2p , (2.95) 
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where again we have expanded the exponents to linear order in p and we are 
neglecting fc-independent numerical coefficients. (Incidentally, it should be borne 
in mind that in a collapsing universe the terms 'growing' and 'decaying' are gauge- 
dependent. The behaviour $ ~ t~ 1_p is the growing mode in Newtonian gauge, 
yet the same physical fluctuation 'decays' as ( ~ t 1_3p in the gauge specified by 
(2.54)). 

2.5 Frontiers 

2.5.1 Matching at the bounce 

Having seen how perturbations are generated in a collapsing universe, we must 
now ask how these perturbations relate to observable quantities in the present-day 
universe. 

In the expanding phase of the universe t is positive and tending towards +oo, 
with p — 1/2 during radiation dominance, followed by p = 2/3 during matter 
dominance. From the analysis leading to (2.94) and (2.95) (but without expand- 
ing the exponents), we see that the dominant contribution to both $ and ( at 
late times comes from the time-independent piece, and that both time-dependent 
pieces decay to zero. The $ ~ ( ~ const, mode is thus the 'growing' mode 
in an expanding universe, and provides the dominant contribution to physical 
quantities as modes re-enter the Hubble radius. Compatibility with observation 
therefore requires this mode to have a scale-invariant spectrum. 

From (2.94), however, we see that the time-independent term in <3> does not 
have a scale-invariant spectrum in the collapsing phase. Nevertheless, it can be 
argued that the matching rules across the bounce induce a generic mixing between 
the two solutions, allowing the scale-invariant spectrum of the $ ~ t~ 1 ~ p mode 
in the collapsing phase to be inherited by the $ ~ const, mode in the expanding 
phase [2, 12]. In other words, the mixing at the bounce must match a piece of 
the scale-invariant growing mode in the collapsing phase to the growing mode 
in the expanding phase, even though from a naive perspective these modes are 
orthogonal. 
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A proposal realising this in the full five-dimensional braneworld setup was 
made in [46], based on the analytic continuation of variables defined on time 
slices synchronous with the collision. (For subsequent criticism of this proposal, 
however, see [47]). Since the nature of the matching at the bounce continues to 
be a subject of active research we will not pursue it further here. Instead, as 
explained in the Introduction, our chief concern will be to understand how the 
problem of propagating a scale-invariant spectrum of density perturbations across 
the bounce is modified, when restored to its true higher-dimensional setting. 

2.5.2 Breaking the degeneracy: gravitational waves 

As we have seen, the power spectrum of the Newtonian potential in a collapsing 
universe is identical to that of its expanding dual. Fortunately this degeneracy is 
broken by tensor perturbations, providing a crucial signature for future observa- 
tions. 

Returning to our calculations in Section § 2.3.2 based on the ADM formalism, 
to compute the action for tensor perturbations we can set 

6<p = 0, h ij = a 2 (r)[8 ij + lij }, (2.96) 

where the tensor perturbation 7^ satisfies djjji = ja — and is automatically 
gauge- invariant. Inserting this into the action (2.53) and expanding to quadratic 
order, we find 

S = ^ J dtd 3 x [a^jijjij - adk'Yijdk'Yij]. (2.97) 
Expressing 7^ in a basis of plane waves with definite polarisation tensors as 

/r\ 3 k 

where en = h % e^ = and ef ■(&) e^(fc) = 25 ss i, we see that each polarisation mode 
obeys essentially the same equation of motion as a massless scalar field: 

7" + 2%i + k 2 j = 0. (2.99) 
This can be recast in canonical form by setting x — a J, yielding 

x" + (k 2 - —)x = 0. (2.100) 
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Again, the standard vacuum choice is the Minkowski vacuum of a comoving 
observer in the far past, corresponding to the boundary condition 

X -» (h/2k) 1/2 e~ ikT (2.101) 

as r — > — oo. To solve for \i one need only observe that the equation of motion 
for x is identical to that for v given in (2.70) (since z oc a for constant u, hence 
z"jz = a" /a). The boundary conditions (2.101) and (2.74) are also identical. 
The solution for \ therefore follows from (2.83); 

x = v = \ 2 {h7rx/4k) 1/2 H^{x), (2.102) 

where x = k\r\, and (3 and A 2 are as defined in (2.80) and (2.85) respectively. 
Defining the tensor spectral index ht such that fc 3 |x| 2 oc fc nT_1 , on super- Hubble 
scales the mode freezes and we have, from (2.86), 



n T - 1 = 3 - 2p = 3 



3p 



1 — p 



(2.103) 



This expression is not invariant under p — > 1/p: an expanding universe with 
a given p produces a tensor spectrum which is much redder than the tensor 
spectrum produced in a contracting universe with p — 1/p as, from the above, 
nr < tit — 2. 



2.6 Summary 

In this chapter we have reviewed two independent scenarios for generating the 
scale-invariant spectrum of primordial density fluctuations required to fit obser- 
vation. In the first scenario, inflation, the early universe undergoes a brief period 
of accelerated expansion in which physical wavelengths are stretched outside the 
quasi-static Hubble radius. The alternative is a slowly collapsing universe followed 
by a big crunch to big bang transition, as proposed in the cyclic and ekpyrotic 
models. In this scenario scale-invariant density perturbations are generated in 
the growing mode of the Newtonian potential during the collapsing phase, as the 
Hubble radius shrinks rapidly inside the near-constant physical length scales of 
the perturbations. 
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Thus far, all our considerations have been couched in the framework of four- 
dimensional effective theory. In the following chapter we leave the prosaic world 
of four dimensions, and recast our view of the universe in five dimensions. 
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Chapter 3 
Braneworld gravity 



Now it is time to explain what was before obscurely said: 
there was an error in imagining that all the four elements 
might be generated by and into one another... 
There was yet a fifth combination which God used 
in the delineation of the universe. 

Plato, Timaeus. 

One of the most striking ideas to emerge from string theory is that the universe 
we inhabit may be a brane embedded in, or bounding, a higher-dimensional 
spacetime [48]. The key to this picture is the notion that gravity, being the 
dynamics of spacetime itself, is free to roam in the full higher-dimensional bulk 
spacetime, whereas the usual standard model forces are confined to the brane. In 
this section we review the simplest phenomenological model of such a scenario, 
the Randall-Sundrum model, focusing on its cosmological properties and on its 
four-dimensional effective description. For a broader survey of the braneworld 
literature, see in particular [49-51]. 

3.1 The Randall-Sundrum model 

The first Randall-Sundrum model [6] comprises a pair of four-dimensional positive- 
and negative-tension branes embedded in a five-dimensional bulk with a negative 
cosmological constant. For simplicity, a Z2 symmetry is imposed about each 
brane, so that the dimension normal to the branes is compactified on an S ll /Z 2 
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orbifold. This scenario provides the simplest possible setup incorporating branes 
with a nontrivial warp factor in the bulk. 
The complete action takes the form: 

S = 2f 6 V=g(±mlR-A)-Yl jf {*n\K* + ^ + Matter) , (3-1) 

where the first group of terms represent the action of the bulk (777,5 is the five- 
dimensional Planck mass, R is the five- dimensional bulk Ricci scalar, and A < 
is the bulk cosmological constant), and the second set of terms encode the action 
of two branes of tension a ± 1 each with its own induced metric g^ b (where Latin 
indices run from to 4), and (optional) matter content specified by £ matter- 

To deal with the assumed Z 2 symmetry about each brane we have adopted 
an extended zone scheme in which each brane appears only once, but there are 
two identical copies of the bulk (hence the factor of two multiplying the bulk 
action). Each brane presents two surfaces to the bulk, and therefore contributes 
two copies of the Gibbons-Hawking surface term proportional to K^, the trace of 
the extrinsic curvature K^ b (itself defined as the projection onto the brane of the 
gradient of the outward-pointing unit normal, i.e., K^ b = g ±c a V cnf) . The 
purpose of these surface terms is to cancel out the second set of surface terms 
arising from the integration by parts of the Einstein-Hilbert action (required to 
remove second derivatives of the metric, furnishing an action quadratic in first 
derivatives only). 

3.1.1 The Israel matching conditions 

In practice, however, it is usually simplest to adopt a "cut and paste" approach 
when dealing with branes: given a solution of the bulk equations of motion, cor- 
responding to a solution of five-dimensional gravity with a negative cosmological 
constant, we can create a Z2-brane simply by gluing together two copies of this 
geometry. The total stress-energy T ab induced on the brane is then given in 
terms of the jump in extrinsic curvature across the brane via the Israel matching 
conditions [52]: 

T ah = -m\ {[K ab ]-h ab [K\), (3.2) 
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where h ab denotes the induced metric on the brane. (A simple derivation of 
this result is provided in Appendix A, where we study the Einstein equations in 
the presence of distributional sources). To evaluate the jump [...], we introduce 
a normal vector to the brane and calculate the extrinsic curvature term in the 
bracket on both sides of the brane using the same normal. We then subtract the 
terms calculated on one side (the side from which the normal points away) from 
the terms calculated on the other side (the side to which the normal is pointing). 
It is easy to see that the final result does not depend on the choice of normal: 
reversing the normal changes the order of the subtraction, but also the sign of 
the extrinsic curvatures. 

Alternatively, we could have evaluated the extrinsic curvatures on each side 
using two different normals (pointing away from the brane in each case). Their 
sum is then equal to the jump defined above. In the case of a Z2-symmetric brane, 
the jump therefore reduces to twice the extrinsic curvature at the brane, eval- 
uated using the outward-pointing unit normal. This allows the Israel matching 
conditions to be re-written as 

mlK ab = -~(T ab --h ab T). (3.3) 

If the stress-energy on the brane is fixed in advance, then these equations con- 
strain the embedding of the brane into the bulk geometry. (We will see an explicit 
example of this when we come to consider brane cosmology: the Israel matching 
conditions will determine the trajectory of the brane through the bulk, and hence 
the form of the induced Friedmann equation on the brane). 

3.2 Static solutions 

An obvious solution for the bulk geometry is to take a slice of five-dimensional 
anti-de Sitter (AdS) spacetime, 

ds 2 = dY 2 + e 2Y/L r]^dx^dx v , (3.4) 

where L is the AdS radius (defined as L 2 = — 677tJ/A), and the indices fi and v run 
from to 3. Let us consider bounding the slice of AdS with flat branes located 
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at constant Y = Y ± , with Y + > Y~ . The coordinates x M then parameterise 
the four dimensions tangent to the branes, and Y parameterises the direction 
normal to the branes. Determining the stress-energy on the branes via the Israel 
matching conditions (3.2) as discussed above, a simple calculation shows that 
T^ v = =F (6ml /L) g^ u . To have static branes located at fixed Y ± therefore, we see 
that the branes must be empty apart from a constant tension cr ± = ±6m|/L (so 
thatT± =-<7±<7±). 

With the tensions tuned to these values there is then a continuous one- 
parameter family of static solutions, parameterised by the distance between the 
branes Y + — Y~ . (The dependence on the overall position in the warp factor can 
be removed by a re-scaling of the coordinates tangential to the brane). Heuristi- 
cally, the stability of the system derives from a balance between, on the one hand, 
the cosmological constant and gravitational potential energy of the bulk, and on 
the other hand, the energy stored in the tensions on the branes. The warp factor 
dictates that an increase in the volume of the bulk (rendering the bulk energy 
more negative) is always accompanied by a corresponding increase in the area 
(and hence energy) of the positive-tension brane, and/or a decrease in the area 
(hence increase in the energy) of the negative-tension brane. 

A useful generalisation of the static solution above is found by taking the bulk 
line element to be 

ds 2 = dY 2 + e 2Y/L g^(x)dx"dx v , (3.5) 

where g^x) is any Ricci-flat four- dimensional metric. (In particular, one could 
choose the Schwarzschild solution, leading to the braneworld black string [53]). 
Again, with the brane tensions appropriately tuned, a solution exists for branes 
located at arbitrary constant Y = Y ± . Later, we will use this solution as the 
starting point in our derivation of the four- dimensional effective theory via the 
moduli space approximation. 

3.3 Cosmological solutions 

Another obvious choice for the bulk metric is the AdS-Schwarzschild black hole. 
This satisfies the five-dimensional Einstein equations with a negative cosmolog- 
ical constant, and moreover, has the interesting property that the horizon is a 
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three-dimensional manifold of constant curvature, which may be either positive, 
negative, or zero [54]. (This is just one example of the considerably richer phe- 
nomenology of black holes in higher dimensions). Our suspicions are immediately 
alerted to the possibility of embedding a three-brane of constant spatial curva- 
ture into this bulk geometry, thereby obtaining a cosmological solution of the 
Randall-Sundrum model. 

More precisely, the line element for AdS-Schwarzschild takes the form 

ds 2 = -f{r)dt 2 + r\r)dr 2 + r 2 dQ 2 3 , (3.6) 

where the constant-curvature three-geometry, 

d^ = (i_A; X 2 )^dx 2 + X 2 d^, (3.7) 

is either open, closed or flat according to whether k = — 1, +1 or respectively, 
and where dQ 2 denotes the usual line element on an S 2 . The function f(r) — k — 
fi/r 2 + r 2 /L 2 , where L is related to the bulk cosmological constant A = —Qml/L 2 
as before, and the black hole mass parameter /i is an arbitrary constant. If fi 
is positive we have a black hole in AdS, but we can also choose fi to be zero, 
yielding pure AdS, or even negative, describing a naked singularity in AdS. (If 
/i = 0, the different choices for k then correspond respectively to closed, open or 
flat slicings of AdS). 

In the next section bar one, we will show how to embed a brane in this 
bulk geometry to obtain an induced metric of FRW form. The corresponding 
Friedmann equations then follow from the Israel matching conditions, assuming 
perfect fluid matter on the branes. First, however, we will reassure ourselves that 
AdS-Schwarzschild is indeed the most general possible choice for the bulk metric, 
granted only cosmological symmetry on the branes. 

3.3.1 Generalising Birkhoff's theorem 

We begin by recalling Birkhoff's theorem in four dimensions: that the unique 
spherically symmetric vacuum solution of Einstein's equations is the static 
Schwarzschild geometry. The assumption of spherical symmetry completely de- 
termines the form of the metric, up to two arbitrary functions of time and the 
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radial coordinate. Application of the Einstein equations then fixes these two func- 
tions in terms of the radial coordinate alone, yielding the Schwarzschild solution. 

Physically, Birkhoff's theorem rules out the emission of gravitational waves 
by a pulsating spherically symmetric body. Furthermore, in a configuration con- 
sisting of a number of infinitely thin, concentric spherical shells of non-vanishing 
stress-energy, with everywhere else being held in vacuum, the gravitational field 
at any point is unaltered by the re-positioning of shells both interior and exterior 
to this point (provided only that interior shells do not cross to the exterior and 
vice versa). Gravity reduced to (1 + 1) dimensions by means of a symmetry is 
therefore completely local. 

In the present, five-dimensional context (where we also have a bulk cosmolog- 
ical constant), a similar theorem applies: the assumption of cosmological sym- 
metry on the branes, amounting to three-dimensional spatial homogeneity and 
isotropy, constrains the form of the metric sufficiently that, upon application of 
the Einstein equations, we obtain a unique solution. With the right choice of 
time-slicing, this solution is moreover static. 

In analogy with our discussion of concentric shells, the physical import of 
this generalised Birkhoff theorem is that the bulk geometry seen by a brane is 
completely unaffected by the presence, and even the motion, of any other branes 
placed in the same bulk spacetime. (Except, of course, when collisions occur). 

To derive the theorem explicitly, following [55] (but using different sign con- 
ventions), we start by writing the bulk line element, without loss of generality, 
as 

ds 2 = e 2 °B~ 2 /\-dr 2 + dy 2 ) + 5 2 / 3 d^, (3.8) 

where a(r,y) and B(r,y) are arbitrary functions, and dQ 2 is as defined in (3.7). 
This represents the most general bulk metric consistent with three-dimensional 
spatial homogeneity and isotropy, after we have made use of our freedom to write 
the (r, y) part of the metric in a conformally flat form. We have further chosen 
the exponents with the aim of simplifying the bulk Einstein equations as much 
as possible. The latter read G a b = — (k/m\)g ab = (6/L 2 )g ab , or equivalently, 
R a b = —{^/L 2 )g a b- For the metric (3.8), we then obtain the following system of 
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partial differential equations: 



B yy — B TT = [12L- 2 B 1/3 + 6kB~ 1/3 } e 2<T , (3.9) 
= [2L- 2 B- 2/3 - kB~ A/3 } e 2CT , (3.10) 



a ,yy ° ; 



B^ yy + B TT = 2a jT B T + 2a t yB t y, (3-H) 
B >Ty = a y B jT + a T B tV . (3.12) 



Switching to the light-cone coordinates 



u=~(r-y), v=~(r + y), (3.13) 



we find 



-B tU0 = [12L~ 2 B 1/3 + 6kB~ 1/3 } e 2 \ (3.14) 

-a uv = [2L~ 2 B~ 2/3 - kB~ 4/3 } e 2CT , (3.15) 

B >uu = 2a jU B iU , (3.16) 

B iVV = 2a tV B tV . (3-17) 

The latter two equations can then be directly integrated to give 

e 2 ° = V'(v)B >u = U'(u)B tV , (3.18) 

where U'{u) and V'(v) are arbitrary nonzero functions, and we will use primes to 

denote ordinary differentiation wherever the argument of a function is unique. It 
follows that B and a take the form 

B = B(U{u) + V(v)), e 2a = B'U'V, (3.19) 

which reduces the sole remaining partial differential equation (3.14) to the ordi- 
nary differential equation 

B" + {\2lT 2 B x l 3 + 6kB~ 1/3 )B' = (3.20) 

=> B' + 9L~ 2 B i/3 + 9kB 2/3 = 9/i, (3.21) 

where \i is a constant of integration. 
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The metric then takes the form 

ds 2 = B'U'V'B- 2 ^(-Adudv) + B 2/3 dfi 2 (3.22) 
= -AB- 2/3 B'dUdV + B 2/3 dn 2 3 . (3.23) 

Upon changing coordinates to 

r = B 1/3 , t = 3(V-U), (3.24) 

we recover the AdS-Schwarzschild metric 

ds 2 = -/(r)dt 2 + /" 1 (r)dr 2 + r 2 dfi 2 , (3.25) 

where, from (3.21), 

3.3.2 The modified Friedmann equations 

The trajectory of a brane embedded in AdS-Schwarzschild can be specified in 
parametric form as t — T{t) and r = R(t), where r is the proper time. The 
functions T(r) and R(t) are then constrained by 

g ab u a u b = -ft 2 + f~ l R 2 = -1, (3.27) 

where u a = (T, R, 0) is the brane velocity, dots denote differentiation with respect 
to r, and / = f(r) is as defined in (3.26). Thus 

T = f- 1 (f + R 2 ) 1 / 2 , (3.28) 

and the components of the unit normal vector (defined such that n a u a = and 
n a n a = 1) are 

nt = ±(R,-f- 1 (f + R 2 ) 1/ \0), (3.29) 

where the choice of sign corresponds to our choice of which side of the bulk we 
keep prior to imposing the Z2 symmetry. If we keep the side for which r < R{r) 
(leading to the creation of a positive-tension brane), then the normal points in 
the direction of decreasing r and we must take the positive sign in the expression 
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above. If, instead, we choose to retain the r > R(t) side of the bulk (resulting 
in the formation of a negative-tension brane), then the normal points in the 
direction of increasing r and we must take the negative sign. Either way, the 
four-dimensional induced metric on the brane is 

ds 2 ± = -dr 2 + R 2 (r)dnl (3.30) 

and hence the cosmological scale factor on the brane can be identified with the 
radial coordinate R(t). The motion of the brane through the bulk therefore 
translates into an apparent cosmological evolution on the brane. 

To calculate the brane trajectory, we must solve the Israel matching conditions 
[56, 57]. There are two nontrivial relations; one deriving from the rr component, 
and the other from the 3-spatial components of the extrinsic curvature. Beginning 
with the latter, writing dfi 2 = 7^ dx l dx : > for i, j = 1, 2, 3, the ij components of 
the extrinsic curvature tensor are 



1 ( f 1 f?2\l/2 

K% = V t nf = - n ±a d a9i3 = T U R 1 9ij, (3.31) 



where the 3-metric gfy = r 2/ ~fij, and the superscript ± indicates whether the 
brane is of positive or negative tension. Assuming a perfect fluid matter content 
in addition to the background tension cr ± , the stress-energy tensor on the brane 
takes the form 

T% = {P+P) u a u b + (p - a ± )^, (3.32) 

where p is the energy density of the fluid, p is its pressure, and g^ b is the induced 
metric on the brane. The ij components of the matching condition (3.3) then 
yield the relation 

±mUf + R 2 ) 1 / 2 = 1 -(p + a ± )R. (3.33) 
o 

After squaring this expression, substituting for / using (3.26), replacing the brane 
tensions by their critical values a ± = ±6m>l/L, and relabelling R(t) as a(r), we 
find 

^ = ^ = ±^ + ^-^ + 4- (3-34) 



This is known as the modified Friedmann equation [58], and is the braneworld 
counterpart of the usual Friedmann equation. Provided the energy density is 
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much less than the critical tension (as is typically the case for standard matter 
or radiation at sufficiently late times), the linear term in p dominates over the 
quadratic term. Then, if we set m\L = m\, where is the four-dimensional 
Planck mass, the modified Friedmann equation is a good approximation to its 
conventional four- dimensional analogue. (At least on a positive-tension brane 
where matter gravitates with the correct sign!) The additional /ia~ 4 term is 
known as the dark radiation term, since it scales like radiation but does not 
interact with standard matter directly, on account of its gravitational origin. 

A second relation follows from the rr component of the extrinsic curvature, 
given by 

K TT = K ab u a u b = u a u h S7 a n b = -n c a c , (3.35) 

where the acceleration a c is defined by a c = u fe VftU c . Since a c u c = 0, the accelera- 
tion can also be written as a c = an c , where a = —K TT . Then, since d t is a Killing 
vector of the static background 1 , we have a = (du t /dr). 

The rr component of the Israel matching condition (3.3) therefore reads 

ml K± =±^i-(/ + #)V* = -I (2p + 3p - ^). (3.36) 
R dr 6 

Combining this with our first relation (3.33) then yields the usual cosmological 
conservation equation 

p + 3H(p + p) = 0, (3.37) 
where H = R/R — a/a is the Hubble parameter as above. 

3.4 A big crunch/big bang spacetime 

Having understood the dynamics of a single brane embedded in AdS-Schwarzschild, 
we turn now to the behaviour of a pair of positive- and negative-tension branes 
when embedded in this same bulk spacetime. 

A useful alternative coordinate system for the bulk may be found by setting 
/ _1 (r)dr 2 = dY 2 . Restricting our attention to the spatially flat case in which 

1 For a Killing vector £ c , an^ = dn c f = a c £ c = ^ c u fc VfcU c = u fc Vfc(u c ^ c ) — u c u b Vf,£ c , where 
the last term vanishes by Killing's equation, V(£,£ c ) = 0, hence an^ = u b dbu^ = (du^/dr). 



11 



3. Braneworld gravity 



k = 0, the bulk line element (3.6) then takes the form 

ds 2 = dY 2 - N 2 (Y)dT 2 + b 2 (Y)dx 2 , (3.38) 
where for pure AdS (fi = 0), 

b 2 (Y) = N 2 (Y) = e 2Y/L - (3.39) 
for AdS-Schwarzschild (/i > 0) with horizon at Y = 0, 

6 2 (F) = cosh2F/L, N 2 (Y) = ^Z'ln ' ^ 3 - 40 ) 

cosh(2r /L) 

and for AdS with a naked singularity at y = (// < 0), 

(In the latter two cases, the value of \x has been absorbed into the definition of 
T and the spatial x l ). 

For a brane with trajectory Y = Y ± (T), the induced metric is 

ds 2 ± = -{N 2 ± - Y± 2 )dT 2 + b 2 ± dx 2 , (3.42) 

where N± = NlY 1 ^) and similarly b± = biY^). The proper time r on the brane 
then satisfies 



dr = JiVi - FJ 2 dT. (3.43) 

In these new coordinates, the spatial components of the Israel matching condi- 
tions, (3.33), can be put in the form 

N± A , p±L\ 



^-V^SJV 1 -^ 2 /^, (3.44) 

where to obtain this expression we have performed the coordinate transformation 
explicitly on a case by case basis, starting from (3.33) with the critical tension, 
then using (3.28) (not confusing the T referring to the parameterisation t = T(r) 
with our present coordinate T) and (3.43). 
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3.4.1 Trajectories of empty branes 

In the special case where the branes are empty (p± = 0) as well as flat, the modi- 
fied Friedmann equations (3.34) imply that the bulk must be AdS-Schwarzschild 
with /i > in order to have moving branes (and hence a nontrivial cosmological 
evolution). The relevant bulk solution is then (3.40), and the Israel matching 
condition (3.44) simplifies to 

= ±smh(2Y ± /L) cosh~ 3/2 (2F ± /L), (3.45) 

where we have identified the positive root with the positive-tension brane and vice 
versa. Equivalently, in terms of the dimensionless brane scale factor b± defined 
in (3.40), we have 

Lb ± T = ±(l-bl A ). (3.46) 
This is easily integrated to yield the brane trajectories: 

±j = b± - \ tan- 1 b ± - i In (^|) - C±, (3-47) 

where the C± are constants of integration. From the form of these expressions, we 
see that the branes must inevitably collide at some instant in time. Translating 
our time coordinate so that this collision occurs at T = 0, the constants of 
integration are then identical; 

C± = b - i tan- 1 b - i In , (3-48) 

where bo is the joint scale factor of both branes at the collision. (If desired, 
we could of course set bo = 1 by a re-scaling of the brane spatial coordinates). 
The brane trajectories are illustrated in Figure 3.1. After emerging from the 
collision, the positive-tension brane escapes to the boundary of AdS, while the 
negative-tension brane asymptotes to the event horizon of the bulk black hole. 

3.4.2 Brane-static coordinates 

Our analysis of braneworld cosmological dynamics has centred so far on coordi- 
nate systems in which the bulk is static and the branes are moving. Yet it is 
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Figure 3.1: The background brane scale factors b± plotted as a function of the 
Birkhoff-frame time T, where b± have been normalised to unity at T = 0. In these 
coordinates the bulk is AdS-Schwarzschild. The brane trajectories are then deter- 
mined by integrating the Israel matching conditions, yielding the result (3.47). In 
the limit as T — > oo, the negative-tension brane asymptotes to the event horizon 
of the black hole, while the positive-tension brane asymptotes to the boundary 
of AdS. 
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also possible, however, to find alternative coordinates in which the branes are 
static, and the bulk evolves dynamically with time [46]. This has the advantage 
of simplifying the form of the Israel matching conditions, which will be especially 
useful when we consider cosmological perturbations. The disadvantage is that the 
explicit form of the bulk solution in these coordinates is not known. Neverthe- 
less, it is easy to see that such a transformation from bulk-static to brane-static 
coordinates must exist, as we now show. 

Starting from the Birkhoff-frame metric (3.38), with b and N given by (3.40), 
we change coordinates from Y to Z, where dZ = dY/N and Z is chosen to be 
zero at the collision event. Then we have 

ds 2 = N 2 (-dT 2 + dZ 2 ) + b 2 dx 2 , (3.49) 

where N and b are now functions of Z. Passing to lightcone coordinates, T ± = 
T ± Z, we have 

ds 2 = N 2 (-dT + dT~) + b 2 dx 2 . (3.50) 
Then, under the lightcone coordinate transformation T ± = f±(T ± ) —T ± Z, 

N 2 

ds 2 = — — (-df 2 + dZ 2 ) + b 2 dx 2 . (3.51) 
J+J — 

Setting y = Z and t = ±exp(±T) (to describe post- and pre-collision spacetimes 
respectively), and defining n 2 (t,y)t 2 = N 2 /f' + f'_, the metric takes the general 
form 

ds 2 = -n 2 (t, y)(-dt 2 + t 2 dy 2 ) + b 2 (t, y)dx 2 . (3.52) 

Through a suitable choice of the functions f±, we can always make the branes 
static in the new coordinates. To see this, observe that the new spatial coordinate 
y satisfies 

y= l -[f + {T + Z)-f^T-Z)l (3.53) 

and hence is a solution of the two-dimensional wave equation. From the general 
theory of the wave equation, we can always find a solution y(T, Z) for arbitrarily 
chosen boundary conditions on the branes (which are themselves described by 
timelike curves Z = Z±(T)). In particular, we are free to choose y = y + on the 
positive-tension brane and y = y_ on the negative-tension brane, for constant y±. 
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(Even after this choice there is additional coordinate freedom, since to completely 
specify y(T, Z) we need additional Cauchy data, e.g. on a constant-time slice). 
For the case of empty flat branes discussed in the previous subsection, 



where b 2 = cosh(2Y/L). The brane trajectories b = b±(T) are given implicitly by 
(3.47), from which, in principle at least, one could determine the trajectories in 
the form Z = Z±(T). The two equations 



then determine the necessary coordinate transformation. Equivalently, since the 
y± are assumed to be constant, differentiating with respect to T gives 



f + (T + Z ± (T))(l + V±{T)) — /_(T — Z±(T))(1 - V±(T)), (3.56) 



where V±(T) = Z T are the brane velocities. The solution of these equations is 
not, however, immediately apparent. Later, when we return to this coordinate 
system in Chapter §5, we will simply construct the metric functions n and b in 
(3.52) from scratch, by solving the bulk Einstein equations. 

3.5 Four-dimensional effective theory 

In this section, we will use the moduli space approximation to derive the four- 
dimensional effective theory describing the Randall- Sundrum model at low en- 
ergies. Before we begin, however, it is interesting to consider how the nature 
of the dimensional reduction employed in braneworld theories differs from its 
counterpart in Kaluza-Klein theory. 

3.5.1 Exact versus inexact truncations 

In a dimensional reduction of the Kaluza-Klein type, one expands the higher- 
dimensional fields in terms of a complete set of harmonics on the compact internal 
space, before truncating to the massless sector of the resulting four-dimensional 




(3.54) 



y ± = - [f + (T + Z±(T)) -f-(T- Z±(T))] 



(3.55) 
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theory. Crucially, this truncation is consistent, in the sense that all solutions of the 
lower-dimensional truncated theory are also solutions of the higher-dimensional 
theory. In practice, this means that when one considers the equations of motion 
for the lower-dimensional massive fields prior to truncation, there are no source 
terms constructed purely from the massless fields that are to be retained. Thus, 
if one starts the system off purely in the massless zero mode, it will remain in 
this mode for all time. The subsequent dynamics are then exactly described by 
the four-dimensional effective theory. 

Braneworld gravity, in contrast, does not in general possess a consistent trun- 
cation down to four dimensions: the ansatz used in the dimensional reduction, 
rather than being an exact solution, is typically only an approximate solution 
of the higher- dimensional field equations 1 . In effect, the warping of the bulk 
ensures that the massive higher Kaluza-Klein modes are sourced by the zero 
mode. If the branes are moving at nonzero speed, then generically these higher 
Kaluza-Klein modes will be continuously produced. The regime of validity of 
the four-dimensional effective theory is therefore limited to asymptotic regions in 
which the branes are moving very slowly, or else are very close together (in which 
case the warp factor and tension on the branes become negligible and the theory 
reduces to Kaluza-Klein gravity). 

3.5.2 The moduli space approximation 

The moduli space approximation applies to any field theory whose equations 
of motion admit a continuous family of static solutions with degenerate action. 
This family of static solutions is parameterised by the moduli, which correspond 
to 'flat' directions in configuration space, along which slow dynamical evolution 
is possible. During this evolution, the excitation along other directions is consis- 
tently small, provided these other directions are stable and are characterised by 
large oscillatory frequencies. 

The action on moduli space may be obtained from the full action by inserting 
as an ansatz the functional form of the static solutions, but with the moduli 
promoted from constants to slowly varying functions of spacetime. Variation 

■"■For an interesting counterexample, sec however [59]. 
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with respect to the moduli then yields the equations of motion governing the low 
energy trajectories of the system along moduli space 1 . 

As we have already discussed, the two brane Randall-Sundrum model does 
indeed possess a one-parameter family of degenerate static solutions (3.4), param- 
eterised by a single modulus which is the interbrane separation. The spectrum of 
low energy degrees of freedom therefore consists of a single four-dimensional mass- 
less scalar field corresponding to this modulus, as well as the four-dimensional 
graviton zero mode (captured by promoting rj^ in (3.4) to some generic Ricci-flat 
g^(x), as in (3.5)). 

We will now proceed to calculate the moduli space action for this theory. 
Two approaches are possible: the first is to allow the interbrane separation to 
slowly fluctuate whilst preserving the form of the bulk metric given in (3.5). In 
this approach [62, 63], the kinetic terms in the moduli space action derive from 
the Gibbons-Hawking boundary terms in (3.1). The second method [2, 63] is to 
transform to an alternative coordinate system in which the branes are held at 
fixed locations, with the relevant modulus being instead encoded into the bulk 
geometry. The kinetic terms in the moduli space action then derive purely from 
the bulk Einstein-Hilbert term in (3.1). For ease of computation we will follow 
this second method here, using a brane-static coordinate system that permits 
the bulk Ricci scalar to be evaluated by a standard conformal transformation 
formula. 

We start with the bulk solution (3.5), re-expressed as 

ds 2 = |J (dZ 2 + (3.57) 

where Z = Lexp(—Y/L), and the branes are located at constant Z = Z , with 
Z + < Z~. The coordinate transformation Z = {Z~ — Z + )(z/L) + Z + then maps 
the positive- and negative-tension branes to z = and z = L respectively. The 
bulk line element now takes the form 

ds 2 = [z/L + c}- 2 (dz 2 + giiv {x)dx»dx v ), (3.58) 

^dn fact, the moduli space approximation underpins much of our knowledge of the classical 
and quantum behaviour of solitons. such as magnetic monopoles and vortices [60, 61]. 
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where the dimensionless modulus c = Z + j{Z" — Z + ), and we have re-absorbed 
a factor of L 2 /(Z~ — Z + ) 2 into g^ix). 

To apply the moduli space approximation, we simply promote the modulus c 
to a spacetime function c(x), yielding the variational ansatz 

ds 2 = [z/L + c(x)}~ 2 (dz 2 + g^(x)dx^dx u ). (3.59) 

The five- dimensional Ricci scalar may now be evaluated using the standard con- 
formal transformation formula [45] 

R = n~ 2 [R-2(n- l)g ab V a V b \nn- (n - 2)(n- l)g ab {V a In fi)V fe In Q\ , 

(3.60) 

where R and R are the Ricci scalars of the two metrics g a b and g a b respectively, 
which are themselves related by the conformal transformation g^ = £l 2 g a b, for 
some smooth strictly positive function Q(x a ). In this formula, n denotes the total 
number of spacetime dimensions, and it is understood that covariant derivatives 
are to be evaluated with respect to g ab . In the present case, n = 5, and we will 
take Q = [z/L + c(x)]" 1 . The five-dimensional Ricci scalar of the metric g a t, 
(with line element dz 2 + g fJiU (x)dx lJ- dx u ) then reduces to the Ricci scalar of the 
four-dimensional metric g^, henceforth denoted by R. Thus 

m l I V 73 ^ 5 ^ = ™>l J V^gtt 3 (R+ 12£ a6 (V a lnfi)V fe lnfi) 

-8m 3 5 J ^~g [n 3 V, In Vt] ^ , (3.61) 

where we have integrated by parts yielding a boundary term on the branes. (The 
boundary terms at spatial infinity, x^ — > oo, are assumed to vanish). Integrating 
over the ^-coordinate, we obtain 

m\L jf Q (Q 2 + - Q 2 _)R + (fi* - Qi) (3(dc) 2 - 5L~ 2 )^ , (3.62) 

where 

tt + (x) = c~ l (x), Q_(x) = (1 + c(x))-\ (3.63) 
are the dimensionless scale factors on the (±) branes respectively. 
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Dispensing with the tildes, the remaining terms in the action (3.1) are now 
easily computed: 

|v°fl(-2A) = ML' 1 f y/=g (3.64) 
~J2j ± y^ a± = SmlL- 1 (3.65) 
-Y.J ± ^-^ 2m l K± = ML- 1 J^(nX-{li), (3.66) 

where A = — 6m|/L 2 , and we are assuming the critical tuning <j ± = ±6m|/L. 
The total moduli space action is thus 

S mod = m\ jf Q R{n\ - fi 2 „) + 3(dc) 2 (nl - fi 4 _)) , (3.67) 

where we have set m\L = m 2 to recover the four- dimensional Planck mass m^. 
(Note that the terms proportional to m\jL all sum to zero). The Einstein frame 
action is then obtained by a further conformal transformation (this time setting 
n = 4 in (3.60)), yielding 

S mod = m lJ 4 ^9 (R ~ 6(<9c) 2 (l + 2c)~ 2 ) . (3.68) 

Setting (f) = —a/3/2 In (1 + 2c), so that — oo < (j) < 0, the moduli space action 
takes the simple form 

Smod = \rniJ^(R- («90) 2 ) , (3.69) 

corresponding to gravity and a minimally coupled massless scalar field. The 
induced metrics on the branes (to which matter couples) are given by 

g+ v = cosh 2 (0/v / 6)^, g^ = S3nh 2 (<f>/y/6)g lu/ , (3.70) 

and the proper separation of the branes is 

d= [ dziz/L + c)- 1 = Llncoth(-0/ v / 6), (3.71) 
Jo 
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so that the branes collide as (j) — > — oo. In this limit, g^ u ~ exp(— a/2/3 0) </ Mi ,, 
whereupon the four- dimensional effective theory reduces to standard Kaluza- 
Klein theory. 

The moduli space approximation is expected to be valid whenever the motion 
of the system is sufficiently slow that non-moduli degrees of freedom are only 
slightly excited. In practice, this means the typical scale of curvature on the 
branes must be much less than the AdS curvature scale (i.e., four-dimensional 
derivatives must be small compared to 1 /L) . In the case of cosmology, this implies 
that the four-dimensional Hubble parameter H must be small compared to 1/L. 
Equivalently, the matter density on the branes must be much less than the brane 
tensions, since p ~ m\LH 2 <C m\jL ~ a. 

Finally, it should be borne in mind that the above derivation of the four- 
dimensional effective theory via the moduli space approximation, while possessing 
the virtue of simplicity, is far from rigorous. In particular, performing the varia- 
tion of the full five- dimensional action with a restricted ansatz is not guaranteed 
to converge on correct solutions of the field equations 1 . 

Nevertheless, the validity of the four- dimensional effective action (3.69) up 
to quadratic order in derivatives can be justified by more rigorous (yet more 
laborious) techniques in which the bulk field equations are solved in a controlled 
expansion, before backsubstituting into the action and integrating out the extra 
dimension [65-68]. Alternatively, one could start out with a more general metric 
ansatz for the bulk, then attempt to restrict the form of the metric functions 
via the bulk field equations and the requirement that the dimensionally reduced 
action be in Einstein frame [69]. 

In the next chapter we will take a fresh approach, identifying a hidden sym- 
metry of the four-dimensional effective action that completely constrains its form 
up to quadratic order in derivatives. 



1 In fact, the moduli space ansatz, which may equivalently be written ds 2 = T(x) 2 dY 2 + 
cxp(2T(x)Y/L) gfj,v(x) dx^dx", is not even sufficiently general to describe linearised perturba- 
tions about a flat background [64] . 
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Chapter 4 



Conformal symmetry of 
braneworld effective actions 



In my work, I have always tried to unite the true with the beautiful; but 
when I had to choose one or the other, I usually chose the beautiful. 



Hermann Weyl 



In this chapter, we show how the brane construction automatically implies 
conformal invariance of the four- dimensional effective theory. This explains the 
detailed form of the low energy effective action previously found using other meth- 
ods. The AdS/CFT correspondence may then be used to improve the effective 
description. We show how this works in detail for a positive- and negative-tension 
brane pair. 

4.1 Derivation of the effective action 

For the general, non-static solution to the Randall- Sundrum model (3.1) it is 
convenient to choose coordinates in which the bulk metric takes the form 



ds 2 = dY 2 + g^ v (x, Y)dx^dx v . 



(4.1) 



The brane loci are now Y ± {x) and the metric induced on each brane is 



9%{x) = d IM Y ± (x)d I/ Y ± (x) + g^{x, Y ± 



(*)). 



(4.2) 
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At low energies we expect the configuration to be completely determined by the 
metric on one brane and the normal distance to the other brane, Y + — Y~. That 
is, we are looking for a four- dimensional effective theory consisting of gravity plus 
one physical scalar degree of freedom. What we will now show is that this theory 
may be determined on symmetry grounds alone. (See also [70] for related ideas). 

The full five- dimensional theory is diffeomorphism invariant. This invariance 
includes the special set of transformations 

Y' = Y + i\x), x'» = x" + ^{x,Y), (4.3) 

with Y) satisfying 

d Y e(x, Y) = -sT{x, Y)d u e(x), (4.4) 

which preserve the form (4.1) of the metric. The transformation (4.3) displaces 
the Y ± (x) coordinates of the branes, 

Y ± (x) -> Y ± (x) + £ 5 - FdrY* (x) , (4.5) 

and alters g^ y (x,Y) via the usual Lie derivative. Using (4.4), one finds that the 
combined effect on each brane metric (4.2) is the four-dimensional diffeomorphism 

x'^ = x^ + e{x,Y ± {x)). (4.6) 

In fact, by departing from the gauge (4.1) away from the branes, we can construct 
a five- dimensional diffeomorphism for which £ M vanishes on the branes. To see 
this, we can set 

e(x, Y) = (x, Y) - f(Y)$(x, r + ), (4.7) 

where £q(x, Y) is the solution to (4.4) which vanishes on the negative-tension 
brane and f(Y) is a function chosen to satisfy f(Y~) = 0, f{Y + ) = 1, and 
f'(Y-) = f'(Y+) = for all x. 

We conclude that the four- dimensional theory, in which Y ± {x) are represented 
as scalar fields, must possess a local symmetry £ 5 (x) acting nontrivially on these 
fields. The dimensionless exponentials ^(x) = exp(y ± (x)/L) transform as con- 
formal scalars: ^(x) — > exp(£ 5 /L) ?/; ± (x), while the induced brane metrics g^ v 
remain invariant. The only local, polynomial, two-derivative action possessing 
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such a symmetry involves gravity with two conformally coupled scalar fields. Af- 
ter diagonalising and re-scaling the fields, this may be expressed as 



m 2 



where A = □ — c± = ±1, and m is a constant with dimensions of mass. It 
should be stressed that the metric g^ v appearing in this expression is that of the 
effective theory, which is in general different to g^ u , the induced metric on the 
branes. Potential terms are excluded by the fact that flat branes, with arbitrary 
constant V ,± ; are solutions of the five-dimensional theory, i.e. the ip^ are moduli. 
By construction, the theory possesses local conformal invariance under 

^± ^ n(z)-V± 9lu , _> V- ( 4 - 9 ) 

For c + = — c_, without loss of generality we can set c + = —1. Provided 
(ip + ) 2 — {ip~) 2 > 0, we obtain the usual sign for the Einstein term, so there 
are no ghosts in the gravitational sector. We can then set tp + = y4cosh0/\/6 
and ip~ = — A sinh (j)/y/E. The field A has the wrong sign kinetic term, but it 
can be set equal to a constant by a choice of conformal gauge. Therefore, in this 
case there are no physical propagating ghost fields. In contrast, a similar analysis 
reveals that when c + = c_ the theory possesses physical ghosts, either in the 
gravitational wave sector (wrong sign of R) or in the scalar sector, no matter how 
the conformal gauge is fixed. We conclude that the low energy effective action 
must be 

m 2 J d A Xy/=g (-^ + A^ + + ^~A^/r) . (4.10) 

We know from the above argument that the brane metrics are conformally 
invariant: from this, and from general covariance, they must equal g^ u times 
homogeneous functions of order two in ip + and ip~ . Yet in the model under 
consideration, we have static solutions gz, = exp(2Y ± / L) r]^ u for all Y + > Y~ . 
The only choice consistent with this, and with (ip + ) 2 — {ip~) 2 > 0, is 



9%= K -^~9^ (4.11) 



which is a conformally invariant equation. (The numerical factor has been intro- 
duced for later convenience). 
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4.1.1 Conformal gauges 

It is instructive to fix the conformal gauge in several ways. First, set ip + = \/6, 
so that g yw = g+ u and the metric appearing in (4.10) is actually the metric on 
the positive-tension brane. The action (4.10) then consists of Einstein gravity 
(with Planck mass m) plus a conformally invariant scalar field ip~ which has to 
be smaller than 

m 2 J d 4 x (r ) 2 ) R + ~ W) a ) • (4-12) 

Changing variables to x = 1 — (V' _ ) 2 /6 produces the alternative form [66, 67] 

m 2 J d 4 xv^ (xR + - ^T^y(^) 2 ) • ( 4 - 13 ) 

Conversely, if we set ip~ = \/6, then g^ u is the metric on the negative-tension 
brane and which has to be larger than \/6, is a conformally coupled scalar 
field. (The relative sign between the gravitational and kinetic terms in the action 
is now wrong, however, and so this gauge possesses ghosts). If we add matter 
coupling to the metric on the positive- and negative-tension branes, we find that 
matter on the negative-tension brane couples in a conformally invariant manner 
to the positive-tension brane metric and the field ip~ , and conversely for matter on 
the positive-tension brane. (Note that we are not implying conformal invariance 
of the matter itself: it is simply that matter coupled to the brane metrics will 
be trivially invariant under the transformation (4.9), as the brane metrics are 
themselves invariant). 

A third conformal gauge maps the theory to Einstein gravity with a minimally 
coupled scalar field <fr, taking values in the range — oo < <ft < 0. Starting from 
(4.10), we can set = A cosh (j)/\/6 and ip~ — ~ A sinh (j)/y/E, as noted earlier, 
yielding the action 

AAA + ^-idcf)) 2 ^ . (4.14) 
= "\/6, we find 

tfOR+W), (4.15) 



— m 2 J d 4 x v / — ~g~ 
Now, choosing the conformal gauge A 

m 2 / d 4 XA 
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i.e., gravity plus a minimally coupled massless scalar. In this gauge equations 
(4.11) read: 

g+, = cosh 2 (0/\/6) fl>, g^ u =mnh 2 (<l>/V6)g IJttf , (4.16) 

in agreement with explicit calculations in the moduli space approach [2]. 

The present treatment also goes some way towards explaining the moduli 
space results. For example, the fact that the moduli space metric is flat is seen to 
be a consequence of conformal invariance. Specifically, for solutions with cosmo- 
logical symmetry, one can pick a conformal gauge in which the metric is static. 
The scale factors on the two branes are determined by ip . From (4.10), the 
moduli space metric is just two-dimensional Minkowski space. 

A couple of results for conformal gravity follow from the above discussion. 
Firstly, in the ^ + = y/E gauge, we have = — v^tanh^/v^- Any solution for 
a minimally coupled scalar 0, with metric g^ v , thus yields a corresponding solution 
for a conformally coupled scalar ip~ , with \ip~\ < \/6 and metric g+ v as in (4.16), 
and vice versa. Secondly, in the ip~ = a/6 gauge, we have tp + = — \/6 coth0/\/6 
hence we may also obtain a solution for a conformally coupled scalar with 
\ip + \ > and metric g~ u given in (4.16). Solutions to conformal scalar gravity 
therefore come in pairs: if g^ v and i/j are a solution, then (ijj 2 /6) g^ and ip = 
is another solution. In terms of branes, this merely states that if g+ v and ip~ are 
known in the gauge ip + = \/6, then it is possible to reconstruct g~ and ip + in 
the gauge = \/6- 

4.1.2 Generalisation to other models 

The argument given above establishing the conformal symmetry of the effective 
action is of a very general nature: the only step at which we specialised to the 
Randall-Sundrum model was in the identification of the brane metrics in terms 
of the effective theory variables (4.11). This required merely the knowledge of a 
one-parameter family of solutions. 

To derive the effective theory for other brane models, it is only necessary to 
generalise this last step. For example, in the case of tensionless branes compact- 
ified on an S 1 /^, the bulk warp is absent, and so we know that a family of 
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static solutions is given by the ground state of Kaluza-Klein theory (in which all 
fields are independent of the extra dimension, hence the additional Z 2 orbifolding 
present in the case of tensionless branes is irrelevant). Ignoring the gauge fields, 
the Kaluza-Klein ansatz for the five- dimensional metric is 

ds 2 = e 2 V^^)dy 2 + e-V^<K*) g ^( x )dx»dx\ (4.17) 

where <p and g^ v extremise an action identical to (4.15). For branes located at 
constant y, the induced metrics are exp(— \/2/3 0) g^, independent of y. 

Using the effective action in the form (4.14), conformal invariance of the in- 
duced brane metrics dictates that 

g% = A 2 f ± {ct>)g^ (4.18) 

for some unknown functions f ± . Upon fixing the conformal gauge to A = \^6, 
one recovers the action (4.15), which is just the standard Kaluza-Klein low energy 
effective action. The functions f^ 1 are thus both equal to (1/6) exp(— a/2/3 0), 
and we have 

9% = e-v 7 ^ V = \^ + + r) V- (4-19) 

Note that this is consistent with the — > — oo limit of the Randall-Sundrum 
theory (4.16): as the brane separation goes to zero, the warping of the bulk 
becomes negligible and the Randall-Sundrum theory tends to the Kaluza-Klein 
limit [46]. 

Conceivably, arguments similar to the above might also be used to derive the 
low energy four-dimensional effective action for the Hofava-Witten model (see 
Chapter §7), although we will not pursue this connection further here. 

4.1.3 Cosmological solutions 

We now turn to a discussion of the general cosmological solutions representing 
colliding branes. We choose a conformal gauge in which the metric is static, and 
all the dynamics are contained in ip . For flat, open and closed spacetimes, the 
spatial Ricci scalar is given by R = 6k, where k — 0, — 1 and +1 respectively. 
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The action (4.10) yields the equations of motion 

$± = -kip*, (4.20) 
(^+) a -(^-) 2 = -fc((^ + ) 2 -(^) 2 ). (4.21) 

For /c = 0, we have the solutions 

ip + = -At + B, ip~ = At + B, t < (4.22) 

representing colliding flat branes. It is natural to match ip + to ip~ across the 
collision, and vice versa, to obtain ip* = ±At + B for t > 0. This solution then 
describes two branes which collide and pass through each other, with the positive- 
tension brane continuing to a negative-tension brane, and vice versa [3, 46]. 
For k — — 1, we have the three solutions 

?/> (1) = Asinht; A cosh t; Ae\ , , 

^ = Asinh(t-t ); Acosh(t-to); Ae*"* , ^ } 

where we set equal to the greater, and equal to the lesser, of ip^ and ip^ . 
For k = +1, we find the bouncing solutions ip^ = Asmt, ip^ = Asm(t — to). 
In the absence of matter on the negative-tension brane, the sin and sinh solutions 
are singular when the negative-tension brane scale factor a_ vanishes. Matter on 
the negative-tension brane scaling faster than gC 4 , however, (e.g. scalar kinetic 
matter) causes the solution for ip~ to bounce smoothly at positive a_, because ip~ 
has a positive kinetic term. This bounce is perfectly regular. In contrast, the big 
crunch/big bang singularity, occurring when the positive- and negative-tension 
branes collide, is unavoidable. 

The above example illustrates a general feature of the brane pair effective 
action. If the positive- and negative-tension brane solutions are continued through 
the collision without re-labelling (this means that the orientation of the warp must 
flip), then the four-dimensional effective action changes sign. The re-labelling 
restores the conventional sign. The same phenomenon is seen in string theories 
obtained by dimensionally reducing eleven-dimensional supergravity, when the 
eleventh dimension collapses and reappears. For a discussion of brane world 
black hole solutions with intersecting branes, see [71]. 
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4.2 AdS/CFT 

Recently, it has been shown that the AdS/CFT correspondence [72, 73] pro- 
vides a powerful approach to the understanding of braneworlds. For a single 
positive-tension brane, the four-dimensional effective description comprises sim- 
ply Einstein gravity plus two copies of the dual CFT [74] (as the Z 2 symmetry 
implies there are two copies of the bulk). Notable successes of this program in- 
clude reproducing the 0(l/r 3 ) corrections to Newton's law on the brane [75], and 
reproducing the modified Friedmann equation induced on the brane [76, 77]. 

Consider, for simplicity, a single positive-tension brane containing only radi- 
ation. Taking the trace of the effective Einstein equations, we find 

-R = 2(87rG 4 ) < T CFT >, (4.24) 

as the stress tensor of the radiation is traceless. The trace anomaly of the dual 
N = 4 SU (N) super- Yang Mills theory must then be evaluated. With the help 
of the AdS/CFT dictionary, this quantity may be calculated for the case of cos- 
mological symmetry as shown in [78], giving 

~ R =Y ( R ^ - \&) • ( 425 ) 

Here, the usual R 2 counterterm has been added to the action in order to eliminate 
the OR term in the trace, thus furnishing second order equations of motion. 
For a cosmological metric, with scale factor a, this becomes 

2(aa + ka 2 ) = L 2 (k + h 2 )h, (4.26) 

where h = a/a and the dot denotes differentiation with respect to conformal 
time. Re-expressing the left-hand side as h~ 1 dt(a 2 + ka 2 ), we can then integrate 
to obtain 

h 2 + k = 1 (B - \ k 2 L 2 ) + -(h 2 + k) 2 L (4.27) 
a z 4 4 a 1 

where B is an integration constant. Now, we can expect to recover Einstein 
gravity on the brane in the limit when L — > 0, with other physical quantities 
held fixed. Expanding all terms in powers of L, at leading order we must obtain 
four-dimensional Einstein gravity, for which 87rG4 = 87TG5/L. We therefore set 
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B ~ (87rG 5 po/3£) + C, where p = po/a 4 is the energy density of conventional 
radiation, and C is a constant independent of L as L — > 0. From (4.27), we then 
obtain the first correction to h 2 + k, namely 

3Lcr or ooa D 
which, thanks to the CFT contribution, now includes the well-known dark energy 
and p 2 corrections [58]. 

It should come as no surprise that the AdS/CFT correspondence only ap- 
proximates the Randall- Sundrum setup up to first nontrivial order in an expan- 
sion in L. The AdS/CFT scenario involves string theory on AdS*, x £5. Since 
a' ~ £ 2 ~ L 2 at fixed 't Hooft coupling, and the masses squared of the Kaluza- 
Klein modes on the S5 are of order 1/L 2 , we expect nontrivial corrections at sec- 
ond order in an expansion in L. Furthermore, one can show from the AdS/CFT 
dictionary that, in order for the p 2 term to dominate in the modified Friedmann 
equation, the temperature of the conventional radiation must be greater than 
the Hagedorn temperature of the string. Clearly, the AdS/CFT correspondence 
cannot describe this situation. 

Let us now consider the extension of the AdS / CFT approach to the case of a 
pair of positive- and negative-tension branes, using the ideas developed in earlier 
in this chapter. The effective action for a single positive-tension brane is 

— J^- J d 4 x^/^R + + 2W CFT [g + ] + S m [g + }, (4.29) 

where g+ is the induced metric on the brane, S m is the brane matter action, and 
Wcft is the CFT effective action (including the appropriate R 2 counterterms) . 
Substituting now for g+ using (4.11), the Einstein-Hilbert term \Z—g + R + be- 
comes — y/—g ip + Aip + . A negative tension brane may then be incorporated as 
follows: 



16nG 



l — J d 4 xv^ (-^ + AV> + + tj;-Ajp-) + 2W C FT[g + ] 

-2W C FT[g-} + S m [g + ] + S m [g-}. (4.30) 



The action for the positive- and negative-tension brane pair must take this form 
in order to correctly reproduce the Friedmann equation for each brane. To see 
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this, consider again the conformal gauge in which the effective theory metric is 
static and all the dynamics are contained in ip^, which play the role of the brane 
scale factors. Variation with respect to the ^ yields the scalar field equations 

(^ ± )" 3 A^ ± = 2(87rG 4 ) < T± FT >, (4.31) 

where the trace anomaly must be evaluated on the induced brane metric g^ u , 
but A is evaluated on the effective metric g^ v . The left-hand side evaluates to 
— (tp ± )~ 3 ('ip ± + k(i/j ± ) 2 ). After identifying ^/y/E with a± according to (4.11), 
we recover equation (4.26), upon dropping the plus or minus label. From the 
necessity of recovering the Friedmann equation on each brane, we may also deduce 
that cross-terms in the action between and ip~ are forbidden. 

The signs associated with the gravitational parts of the action are required 
to achieve consistency with (4.10). Consequently, the relative sign between the 
gravity plus CFT part of the action, and that of the matter, is reversed for the 
negative-tension brane, consistent with the modified Friedmann equations [58] 



± 3L 36 a 2 a 4 ' 



where plus and minus label the positive- and negative-tension branes, and C is 
again a constant representing the dark radiation. 

In summary, we have elucidated the origin of conformal symmetry in brane 
world effective actions, and shown how this determines the effective action to low- 
est order. When combined with the the AdS/CFT correspondence, our approach 
also recovers the first corrections to the brane Friedmann equations. 
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Chapter 5 

Solution of a braneworld big 
crunch/big bang cosmology 

We can lick gravity, but sometimes the paperwork is overwhelming. 

Wernher von Braun 

In this chapter we solve for the cosmological perturbations in a five-dimensional 
background consisting of two separating or colliding boundary branes, as an ex- 
pansion in the collision speed V divided by the speed of light c. Our solution 
permits a detailed check of the validity of four- dimensional effective theory in 
the vicinity of the event corresponding to the big crunch/big bang singular- 
ity. We show that the four- dimensional description fails at the first nontrivial 
order in (V/c) 2 . At this order, there is nontrivial mixing of the two relevant 
four- dimensional perturbation modes (the growing and decaying modes) as the 
boundary branes move from the narrowly-separated limit described by Kaluza- 
Klein theory to the well-separated limit where gravity is confined to the positive- 
tension brane. We comment on the cosmological significance of the result and 
compute other quantities of interest in five-dimensional cosmological scenarios. 

5.1 Introduction 

Two limiting regimes can be distinguished in which the dynamics of the Randall- 
Sundrum model simplify: the first is the limit in which the interbrane separation 
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is much greater than the AdS curvature radius, and the second is the limit in 
which the interbrane separation is far smaller. When the two boundary branes 
are very close to one another, the warping of the five-dimensional bulk and the 
tension of the branes become irrelevant. In this situation, the low energy modes 
of the system are well-described by a simple Kaluza-Klein reduction from five to 
four dimensions, i.e., gravity plus a scalar field (the Z2 projections eliminate the 
gauge field zero mode). When the two branes are widely separated, however, the 
physics is quite different. In this regime, the warping of the bulk plays a key 
role, causing the low energy gravitational modes to be localised on the positive- 
tension brane [8, 79, 80]. The four- dimensional effective theory describing this 
new situation is nevertheless identical, consisting of Einstein gravity and a scalar 
field, the radion, describing the separation of the two branes. 

In this chapter, we study the transition between these two regimes - from 
the naive Kaluza-Klein reduction to localised Randall-Sundrum gravity - at fi- 
nite brane speed. In the two asymptotic regimes, the narrowly-separated brane 
limit and the widely-separated limit, the cosmological perturbation modes show 
precisely the behaviour predicted by the four-dimensional effective theory. There 
are two massless scalar perturbation modes; in longitudinal gauge, and in the 
long wavelength (k — > 0) limit, one mode is constant and the other decays as 
tj 2 , where £4 is the conformal time. In the four- dimensional description, these 
two perturbation modes are entirely distinct: one is the curvature perturbation 
mode; the other is a local time delay to the big bang. Nonetheless, we shall show 
that in the five-dimensional theory, at first nontrivial order in the speed of the 
brane collision, the two modes mix. If, for example, one starts out in the time 
delay mode at small £4, one ends up in a mixture of the time delay and curvature 
perturbation modes as £4 — > 00. Thus the two cosmological perturbation modes 
- the growing and decaying adiabatic modes - mix in the higher- dimensional 
braneworld setup, a phenomenon which is prohibited in four dimensions. 

The mode-mixing occurs as a result of a qualitative change in the nature of 
the low energy modes of the system. At small brane separations, the low energy 
modes are nearly uniform across the extra dimension. Yet as the brane separa- 
tion becomes larger than the bulk warping scale, the low energy modes become 
exponentially localised on the positive-tension brane. If the branes separate at 
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finite speed, the localisation process fails to keep pace with the brane separation 
and the low energy modes do not evolve adiabatically. Instead, they evolve into a 
mixture involving higher Kaluza-Klein modes, and the four-dimensional effective 
description fails. 

The mixing we see between the two scalar perturbation modes would be pro- 
hibited in any local four-dimensional effective theory consisting of Einstein grav- 
ity and matter fields, no matter what the matter fields were. The mixing is 
therefore a truly five- dimensional phenomenon, which cannot be modelled with 
a local four-dimensional effective theory. There is, moreover, an independent ar- 
gument against the existence of any local four- dimensional description of these 
phenomena. In standard Kaluza-Klein theory, it is well known that the entire 
spectrum of massive modes is actually spin two [81]. Yet, despite many attempts, 
no satisfactory Lagrangian description of massive, purely spin two fields has ever 
been found [82, 83]. Again, this suggests that one should not expect to describe 
the excitation of the higher Kaluza-Klein modes in terms of an improved, local, 
four-dimensional effective theory. 

The system we study consists of two branes emerging from a collision. In 
this situation, there are important simplifications which allow us to specify ini- 
tial data rather precisely. When the brane separation is small, the fluctuation 
modes neatly separate into light Kaluza-Klein zero modes, which are constant 
along the extra dimension, and massive modes with nontrivial extra-dimensional 
dependence. Furthermore, the brane tensions and the bulk cosmological con- 
stant become irrelevant at short distances. It is thus natural to specify initial 
data which map precisely onto four-dimensional fields in the naive dimensionally- 
reduced theory describing the limit of narrowly-separated branes. With initial 
data specified this way, there are no ambiguities in the system. The two branes 
provide boundary conditions for all time and the five-dimensional Einstein equa- 
tions yield a unique solution, for arbitrary four-dimensional initial data. 

Our main motivation is the study of cosmologies in which the big bang was a 
brane collision, such as the cyclic model [1] . Here, a period of dark energy dom- 
ination, followed by slow contraction of the fifth dimension, renders the branes 
locally flat and parallel at the collision. During the slow contraction phase, grow- 
ing, adiabatic, scale-invariant perturbations are imprinted on the branes prior to 
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the collision. Yet if the system is accurately described by four-dimensional effec- 
tive theory throughout, then, as a number of authors have noted [22, 47, 84-87], 
there is an apparent roadblock to the passage of the scale-invariant perturba- 
tions across the bounce. Namely, it is hard to see how the growing mode in the 
contracting phase, usually described as a local time delay, could match onto the 
growing mode in the expanding phase, usually described as a curvature pertur- 
bation. In this chapter, we show that the four-dimensional effective theory fails 
at order (V/c) 2 , where V is the collision speed and c is the speed of light. The 
four-dimensional description works well when the branes are close together, or 
far apart. As the branes move from one regime to the other, however, the two 
four-dimensional modes mix in a nontrivial manner. 

The mixing we find demonstrates that the approach of two boundary branes 
along a fifth dimension produces physical effects that cannot properly be modelled 
by a local four- dimensional effective theory. Here, we deal with the simplest 
case involving two empty boundary branes separated by a bulk with a negative 
cosmological constant. For the cyclic model, the details are more complicated. In 
particular, there is an additional bulk stress AT|, associated with the interbrane 
force, that plays a vital role in converting a growing mode corresponding to a 
pure time delay perturbation into a mixture of time delay and curvature modes 
on the brane. We will present some preliminary considerations of the effects 
of such a bulk stress in the following chapter. For now, our main conclusion 
with regard to the cyclic model is that, to compute properly the evolution of 
perturbations before and after a brane collision, one must go beyond the four- 
dimensional effective theory to consider the full five- dimensional theory. 

The outline of this chapter is as follows. In Section §5.2, we provide an 
overview of our three solution methods. In § 5.3, we solve for the background and 
cosmological perturbations using a series expansion in time about the collision. 
In §5.4, we present an improved method in which the dependence on the fifth 
dimension is approximated using a set of higher-order Dirichlet or Neumann poly- 
nomials. In § 5.5, we develop an expansion about the small-fV/c) scaling solution, 
before comparing our results with those of the four-dimensional effective theory 
in § 5.6. We conclude with a discussion of mode-mixing in § 5.7. Detailed explicit 
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solutions may be found in Appendix C, and the Mathematica code implementing 
our calculations is available online [88]. 

5.2 Three solution methods 

In this section, we review the three solution methods employed, noting their 
comparative merits. For the model considered here, with no dynamical bulk fields, 
as we saw in §3.3.1 there is a Birkhoff-like theorem guaranteeing the existence 
of coordinates in which the bulk is static. It is easy to solve for the background 
in these coordinates. The motion of the branes complicates the Israel matching 
conditions, however, rendering the treatment of perturbations difficult. For this 
reason, it is preferable to choose a coordinate system in which the branes are 
located at fixed spatial coordinates y = ±yo, and the bulk evolves with time. 

We shall employ a coordinate system in which the five- dimensional line ele- 
ment for the background takes the form 

ds 2 = n 2 (t, y)(-dt 2 + t 2 dy 2 ) + b 2 (t, y)dx 2 , (5.1) 

where y parameterises the fifth dimension and x l (for % = 1,2,3), the three 
noncompact dimensions. Cosmological isotropy excludes dtdx 1 or dydx 1 terms, 
and homogeneity ensures n and b are independent of x. The t, y part of the 
background metric may then be taken to be conformally flat, and one may further 
choose to write the metric for this two-dimensional Minkowski spacetime in Milne 
form. Since we are interested in scenarios with colliding branes in which the bulk 
geometry about the collision is Milne, we will assume the branes to be located at 
y = ±?/o, with the collision occurring at t = 0. By expressing the metric in locally 
Minkowski coordinates, T = tcoshy and Y = tsinhy, one sees that the collision 
speed is (V/c) = tanh2yo and the relative rapidity of the collision is 2?/o- As long 
as the bulk metric is regular at the brane collision and possesses cosmological 
symmetry, the line element may always be put into the form (5.1). Furthermore, 
by suitably re-scaling coordinates one can choose 6(0, y) — n(0, y) — 1. 

In order to describe perturbations about this background, one needs to specify 
an appropriate gauge choice. Five- dimensional longitudinal gauge is particularly 
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convenient [89]: firstly, it is completely gauge-fixed; secondly, the brane trajecto- 
ries are unperturbed in this gauge [46] , so that the Israel matching conditions are 
relatively simple; and finally, in the absence of anisotropic stresses, the traceless 
part of the Einstein (spatial) equation yields a constraint among the pertur- 
bation variables, reducing them from four to three. In light of these advantages, 
we will work in five-dimensional longitudinal gauge throughout. 
Our three solution methods are as follows: 

• Series expansion in t 

The simplest solution method for the background is to solve for the metric 
functions n(t, y) and b(t, y) as a series in powers of t about t = 0. At each order, 
the bulk Einstein equations yield a set of ordinary differential equations in y, 
with the boundary conditions provided by the Israel matching conditions. These 
are straightforwardly solved. A similar series approach, involving powers of t and 
powers of t times In t suffices for the perturbations. 

The series approach is useful at small times (t/L) <C 1 since it provides the 
precise solution for the background plus generic perturbations, close to the brane 
collision, for all y and for any collision rapidity yo- It allows one to uniquely 
specify four-dimensional asymptotic data as t tends to zero. Nonetheless, the 
series thus obtained fails to converge at quite modest times. Following the system 
to long times requires a more sophisticated method. Instead of taking (t/L) as 
our expansion parameter, we want to use the dimensionless rapidity of the brane 
collision y , and solve at each order in y . 

• Expansion in Dirichlet/Neumann polynomials in y 

In this approach we represent the spacetime metric in terms of variables obey- 
ing either Dirichlet or Neumann boundary conditions on the branes. We then 
express these variables as series of Dirichlet or Neumann polynomials in y and 
yo, bounded at each subsequent order by an increasing power of the collision ra- 
pidity yo- (Recall that the range of the y coordinate is bounded by \y\ < yo). 
The coefficients in these expansions are undetermined functions of t. By solving 
the five-dimensional Einstein equations perturbatively in y Q , we obtain a series 
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of ordinary differential equations in t, which can then be solved exactly. In this 
Dirichlet/Neumann polynomial expansion, the Israel boundary conditions on the 
branes are satisfied automatically at every order in yo, while the initial data at 
small t are provided by the previous series solution method. 

The Dirichlet/Neumann polynomial expansion method yields simple, explicit 
solutions for the background and perturbations as long as (t/L) is smaller than 
1/yo- Since yo <C 1, this considerably improves upon the naive series expansion 
in t. For (t/L) of order 1/yo, however, the expansion fails because the growth 
in the coefficients overwhelms the extra powers of yo at successive orders. Since 
(t/L) ~ 1/yo corresponds to brane separations of order the AdS radius, the 
Dirichlet/Neumann polynomial expansion method fails to describe the late-time 
behaviour of the system, and a third method is needed. 

• Expansion about the scaling solution 

The idea of our third method is to start by identifying a scaling solution, whose 
form is independent of yo for all yo <C 1. This scaling solution is well-behaved 
for all times and therefore a perturbation expansion in yo about this solution is 
similarly well-behaved, even at very late times. To find the scaling solution, we 
first change variables from t and y to an equivalent set of dimensionless variables. 
The characteristic velocity of the system is the brane speed at the collision, V = 
ctanh2yo ~ 2q/o for small yo, where we have temporarily restored the speed of 
light c. Thus we have the dimensionless time parameter x = yod/L ~ Vt/L, of 
order the time for the branes to separate by one AdS radius. We also re-scale the 
y-coordinate by defining to = y/yo, whose range is — 1 < oj < 1, independent of 
the characteristic velocity. 

As we shall show, when re-expressed in these variables, for small yo, the bulk 
Einstein equations become perturbatively ultralocal: at each order in yo one only 
has to solve an ordinary differential equation in u, with a source term determined 
by time derivatives of lower order terms. The original partial differential equations 
reduce to an infinite series of ordinary differential equations in ui which are then 
easily solved order by order in yo- 

This method, an expansion in yo about the scaling solution, is the most power- 
ful and may be extended to arbitrarily long times t and for all brane separations. 
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In light of the generalised Birkhoff theorem, the bulk in between the two branes 
is just a slice of five- dimensional AdS-Schwarzschild spacetime, within which the 
two branes move [27, 50, 51]. (The bulk black hole is itself merely virtual, how- 
ever, as it lies beyond the negative-tension brane and hence is excluded from the 
physical region). As time proceeds, the negative-tension brane becomes closer 
and closer to the horizon of the virtual AdS-Schwarzschild black hole. Even 
though its location in the Birkhoff-frame (static) coordinates freezes (see Figure 
3.1), its proper speed grows and the y expansion fails. Nonetheless, by analytic 
continuation of our solution in uj and x, we are able to circumvent this temporary 
breakdown of the yo expansion and follow the positive-tension brane, and the 
perturbations localised near it, as they run off to the boundary of anti-de Sitter 
spacetime. 

Our expansion about the scaling solution is closely related to derivative- 
expansion techniques developed earlier by a number of authors [66, 90, 91]. In 
these works, an expansion in terms of brane curvature over bulk curvature was 
used. For cosmological solutions, this is equivalent to an expansion in L"K + , where 
'K + is the Hubble constant on the positive-tension brane. In the present instance, 
however, we specifically want to study the time-dependence of the perturbations 
for all times, from the narrowly-separated to the well-separated brane limit. For 
this purpose, it is better to use a time- independent expansion parameter (y ), and 
to include all the appropriate time-dependence order by order in the expansion. 

Moreover, in these earlier works, the goal was to find the four-dimensional 
effective description more generally, without specifying that the branes emerged 
from a collision with perturbations in the lowest Kaluza-Klein modes. Conse- 
quently, the solutions obtained contained a number of undetermined functions. 
In the present context, however, the initial conditions along the extra dimension 
are completely specified close to the brane collision by the requirement that only 
the lowest Kaluza-Klein mode be excited. The solutions we obtain here are fully 
determined, with no arbitrary functions entering our results. 

Returning to the theme of the four- dimensional effective theory, we expect on 
general grounds that this should be valid in two particular limits: firstly, as we 
have already discussed, a Kaluza-Klein description will apply at early times near 
to the collision, when the separation of the branes is much less than L. Here, the 
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warping of the bulk geometry and the brane tensions can be neglected. Secondly, 
when the branes are separated by many AdS lengths, one expects gravity to 
become localised on the positive-tension brane, which moves ever more slowly as 
time proceeds, so the four- dimensional effective theory should become more and 
more accurate. 

Equipped with our five- dimensional solution for the background and perturba- 
tions obtained by expanding about the scaling solution, we find ourselves able to 
test the four- dimensional effective theory explicitly. We will show that the four- 
dimensional effective theory accurately captures the five-dimensional dynamics 
to leading order in the 2/0-expansion, but fails at the first nontrivial order. Our 
calculations reveal that the four- dimensional perturbation modes undergo a mix- 
ing in the transition between the Kaluza-Klein effective theory at early times 
and the brane-localised gravity at late times. This effect is a consequence of the 
momentary breakdown of the effective theory when the brane separation is of 
the order of an AdS length, and cannot be seen from four-dimensional effective 
theory calculations alone. 

5.3 Series expansion in time 

As described above, we find it simplest to work in coordinates in which the 
brane locations are fixed but the bulk evolves. The bulk metric is therefore 
given by (5.1), with the brane locations fixed at y = ±yo for all time t. The 
five-dimensional solution then has to satisfy both the Einstein equations and the 
Israel matching conditions on the branes [52]. 

The bulk Einstein equations read G b a = — AS b , where the bulk cosmological 
constant is A = —6/L 2 (we work in units in which the four-dimensional grav- 
itational coupling 87rG4 = 8irG^/L = 1). Evaluating the linear combinations 
Gq + G\ and G® + G\ — {1/2)G\ (where denotes time, 5 labels the y direction, 
and i runs over the noncompact directions), we find: 
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(5.3) 



70 



5.3 Series expansion in time 



where (t/L) = e T , (3 = 31n& and v = In (nt/L). The Israel matching conditions 
on the branes read [46, 89] 

b n L ' 1 J 

where all quantities are to be evaluated at the brane locations, y = ±y - 

We will begin our assault on the bulk geometry by constructing a series ex- 
pansion in t about the collision, implementing the Israel matching conditions on 
the branes at each order in t. This series expansion in t is then exact in both y 
and the collision rapidity y$. It chief purpose will be to provide initial data for the 
more powerful solution methods that we will develop in the following sections. 

The Taylor series solution in t for the background was first presented in [46]. 
Expanded up to terms of 0(t 3 /L 3 ), 

t 1 t 2 
n = 1 + (sechyo sinh?/) — + - sech 2 y (— 3 + cosh2y + 2 cosh2y) — , (5.5) 

L 4 L z 

t 1 t 2 
b = 1 + (sech y sinh y) — H — sech 2 y (— cosh 2y + cosh 2y) — . (5.6) 

L 2 L 

(Note that in the limit as t — ► we correctly recover compactified Milne space- 
time). 

Here, however, we will need the perturbations as well. Working in five- 
dimensional longitudinal gauge for the reasons given in the previous section, the 
perturbed bulk metric takes the form (see Appendix B) 

ds 2 = n 2 (-(1 + 2$ L ) dt 2 - 2W L dtdy + t 2 (l- 2T L ) dy 2 ) +b 2 (l- 2* L ) dx 2 , 

(5.7) 

with = $l — ^ l being imposed by the five-dimensional traceless equation. 
The Israel matching conditions at y — ±y$ then read 

Tit Tit 

*L,„ = r i — , ® L ,y = -r L -, w L = o. (5.8) 

Performing a series expansion, we find 

$ L = -g + B SeCh y ° Sinh V + (A - I - ^- + V In \kt\ 



t 2 t V 8 4 6 

H Scosh2y ( — 1 + 6yocoth2yo) sech 2 y B sech 2 y sinh 2y ) , (5.9) 

16 8 / 

B sech Un sinh w / , B Bk 2 , B , , 9 \ , 

^ L = ^ -+ (2A-j + — ln|fct| + jcosh2ysech 2 ?/oj, (5.10) 

3 

IVl = — -5 sech 2 y (y cosh 2y — y cosh2y sinh2?/) t, (5.11) 
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where the first two equations are accurate to 0(t) and the third is accurate to 
0(t 2 ). We have moreover set L = 1 for clarity; except for a few specific instances, 
we will now adopt this convention throughout the rest of this chapter. (To restore 
L, simply replace t — > t/L and k — > kL). The two arbitrary constants A and B 
(which may themselves be arbitrary functions of k) have been chosen so that, on 
the positive-tension brane, to leading order in y , $^ goes as 



5.4 Expansion in Dirichlet and Neumann poly- 
nomials 

5.4.1 Background 

Having solved the relevant five-dimensional Einstein equations as a series expan- 
sion in the time t before or after the collision event, we now have an accurate 
description of the behaviour of the bulk at small t for arbitrary collision rapidi- 
ties. In order to match onto the incoming and outgoing states, however, we really 
want to study the long-time behaviour of the solutions, as the branes become 
widely separated. Ultimately, this will enable us to successfully map the system 
onto an appropriate four- dimensional effective description. Instead of expanding 
in powers of the time, we approximate the five- dimensional solution as a power 
series in the rapidity of the collision, and determine each metric coefficient for all 
time at each order in the rapidity. 

Our main idea is to express the metric as a series of Dirichlet or Neumann 
polynomials in y and y, bounded at order n by a constant times yfi, such that 
the series satisfies the Israel matching conditions exactly at every order in y . To 
implement this, we first change variables from b and n to those obeying Neumann 
boundary conditions. From (5.4), b/n is Neumann. Likewise, if we define N(t,y) 



then one can easily check that N(t, y) is also Neumann on the branes. Notice 
that if N and b/n are constant, the metric (5.1) is just that for anti-de Sitter 




(5.12) 



by 



1 



(5.13) 



nt = 



N(t,y)-y 1 
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spacetime. For fixed yo, N describes the the proper separation of the two branes, 
and b is an additional modulus describing the three-dimensional scale-factor of 
the branes. 

Since N and b/n obey Neumann boundary conditions on the branes, we can 
expand both in a power series 

oo oo 

N = N (t) + Y,N n (t)P n (y), b/n = q (t) + $^<Z„(i)P„(y), (5.14) 

n=3 n=3 

where P n {y) are polynomials 

Tl 

Pn(y) = y n ^y n ~ 2 yl n = 3,4,... (5.15) 

n — 2 

satisfying Neumann boundary conditions, each bounded by |P n (y)| < ly^jin — 2) 
for the relevant range of y. Note that the time-dependent coefficients in this 
ansatz may also be expanded as a power series in y§. By construction, our ansatz 
satisfies the Israel matching conditions exactly at each order in the expansion. 
The bulk Einstein equations are not satisfied exactly, but as the expansion is 
continued, the error terms are bounded by increasing powers of yo- 

Substituting the series ansatze (5.14) into the background Einstein equations 
(5.2) and (5.3), we may determine the solution order by order in the rapidity 
yo- At each order in y , one generically obtains a number of linearly independent 
algebraic equations, and at most one ordinary differential equation in t. The 
solution of the latter introduces a number of arbitrary constants of integration 
into the solution. 

To fix the arbitrary constants, one first applies the remaining Einstein equa- 
tions, allowing a small number to be eliminated. The rest are then determined 
using the series expansion in t presented in the previous section: as this solution 
is exact to all orders in y , we need only to expand it out to the relevant order 
in yo, before comparing it term by term with our Dirichlet/Neumann polynomial 
expansion (which is exact in t but perturbative in y ), taken to a corresponding 
order in t. The arbitrary constants are then chosen so as to ensure the equiva- 
lence of the two expansions in the region where both t and yo are small. This 
procedure suffices to fix all the remaining arbitrary constants. 
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The first few terms of the solution are 



and 




(5.19) 



(5.16) 



(5.17) 



(5.18) 



The full solution up to O(?/q ) may be found in Appendix C.l. 
5.4.2 Perturbations 

Following the same principles used in our treatment of the background, we con- 
struct the two linear combinations 



both of which obey Neumann boundary conditions on the branes, as may be 
checked from (5.4) and (5.8). In addition, Wl already obeys simple Dirichlet 
boundary conditions. 

The two Neumann variables, 04 and £4, are then expanded in a series of 
Neumann polynomials and Wl is expanded in a series of Dirichlet polynomials, 



D n (y)=y n -y%, n = 2,4,..., D n (y) = yD^y), n = 3,5,..., (5.21) 



each bounded by \D n (y)\ < y$ for n even and y^n — l)/n n//( - n_1 ^ for n odd, over 
the relevant range of y. As in the case of the background, the time- dependent 
coefficients multiplying each of the polynomials should themselves be expanded 
in powers of y . 

To solve for the perturbations it is sufficient to use only three of the per- 
turbed Einstein equations (any solution obtained may then be verified against 
the remainder). Setting 



04 = 2^ + £ 4 = b 2 (^ L - $ L ) = b 2 T L , 



(5.20) 



W L 



(5.22) 
(5.23) 
(5.24) 
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where t = e T , (3 = 3 In 6 and v = In rat, the Gf , G® and G\ equations take the form 

W>T = 2<f), y -4e 3u/2 {4je u/2 ) t y, (5.25) 
<t>,r = ^w >y -e 3v (^e-") )T , (5.26) 

-^e- 2 ^ 3 (HP 2 r + Pl)-wP, T P, y ). (5.27) 

Using our Neumann and Dirichlet ansatze for </>4, £4 and Wl, the Israel match- 
ing conditions are automatically satisfied and it remains only to solve (5.25), 
(5.26) and (5.27) order by order in the rapidity. The time-dependent coefficients 
for 04, £4 and Wl are then found to obey simple ordinary differential equations, 
with solutions comprising Bessel functions in kt, given in Appendix C.2. Note 
that it is not necessary for the set of Neumann or Dirichlet polynomials we have 
used to be orthogonal to each other: linear independence is perfectly sufficient to 
determine all the time-dependent coefficients order by order in yo. 

As in the case of the background, the arbitrary constants of integration re- 
maining in the solution after the application of the remaining Einstein equations 
are fixed by performing a series expansion of the solution in t. This expansion 
can be compared term by term with the series expansion in t given previously, 
after this latter series has itself been expanded in y . The arbitrary constants are 
then chosen so that the two expansions coincide in the region where both t and 
yo are small. The results of these calculations, at long wavelengths, are 



<D L = A-B[^-^ln\kt\) + (At + j)y + ... (5.28) 



k 2 f B\ 

* L = 2A + B—\n\kt\ - f At + jj y+ ... (5.29) 

W L = 6At 2 (y 2 -y 2 ) + ... (5.30) 

where the constants A and B can be arbitrary functions of k. The solutions for 
all k, to fifth order in yo, are given in Appendix C.2. 
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5.5 Expansion about the scaling solution 

It is illuminating to recast the results of the preceding sections in terms of a set of 
dimensionless variables. Using the relative velocity of the branes at the moment 
of collision, V = 2ctanhyo — 2ct/o (where we have temporarily re-introduced 
the speed of light c), we may construct the dimensionless time parameter x = 
y ct/L ~ Vt/L and the dimensionless ^-coordinate u = y/yo ~ y(c/V). 

Starting from the full Dirichlet/ Neumann polynomial expansion for the back- 
ground given in Appendix C.I, restoring c to unity and setting t = xL/y Q and 
y = uy , we find that 

n~ l = N(x) -ux + 0{y$), (5.31) 
q{x) + 0{yl), (5.32) 



n 



where 



~ . . 1 x 2 3x 4 25x 6 343x 8 2187 x 10 12 , 

N(x) = 1 hO x 12 , (5.33) 

v ; 2 8 48 384 1280 v v ; 

3x 2 7x 4 55 x 6 245 x 8 4617x 10 n . 12 . . 

q(x) = 1 hO r . (5.34) 

qK ' 2 8 48 128 1280 v ; v ) 

The single term in (5.31) linear in u is necessary in order that n" 1 satisfies 
the correct boundary conditions. Apart from this one term, however, we see 
that to lowest order in yo the metric functions above turn out to be completely 
independent of uj. Similar results are additionally found for the perturbations. 

Later, we will see how this behaviour leads to the emergence of a four- 
dimensional effective theory. For now, the key point to notice is that this series 
expansion converges only for x <C 1, corresponding to times t L/yo- In order 
to study the behaviour of the theory for all times therefore, we require a means 
of effectively resumming the above perturbation expansion to all orders in x. Re- 
markably, we will be able to accomplish just this. The remainder of this section, 
divided into five parts, details our method and results: first, we explain how to 
find and expand about the scaling solution, considering only the background for 
simplicity. We then analyse various aspects of the background scaling solution, 
namely, the brane geometry and the analytic continuation required to go to late 
times, before moving on to discuss higher-order terms in the expansion. Finally, 
we extend our treatment to cover the perturbations. 
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5.5.1 Scaling solution for the background 

The key to the our method is the observation that the approximation of small 
collision rapidity (y <C 1) leads to a set of equations that are perturbatively 
ultralocal: transforming to the dimensionless coordinates x and u, the Einstein 
equations for the background (5.2) and (5.3) become 

P, U0J + P%-l2e 2S> = y*(x(xP, x ) tX + x 2 pl), (5.35) 
v,^-\pl + 2e 2 » = y 2 (x(xv x ), x -^x 2 (3 2 x ), (5.36) 

where we have introduced v = v + lny - Strikingly, all the terms involving x- 
derivatives are now suppressed by a factor of y 2 relative to the remaining terms. 
This segregation of x- and cu-derivatives has profound consequences: when solv- 
ing perturbatively in y , the Einstein equations (5.35) and (5.36) reduce to a 
series of ordinary differential equations in u>, as opposed to the partial differential 
equations we started off with. 

To see this, consider expanding out both the Einstein equations (5.35) and 
(5.36) as well as the metric functions j3 and i/asa series in positive powers of yo. 
At zeroth order in y , the right-hand sides of (5.35) and (5.36) vanish, and the 
left-hand sides can be integrated with respect to u> to yield anti-de Sitter space. 
(This was our reason for using v = v + In y rather than v: the former serves to 
pull the necessary exponential term deriving from the cosmological constant down 
to zeroth order in y , yielding anti-de Sitter space as a solution at leading order. 
As we are merely adding a constant, the derivatives of v and v are identical.) 
The Israel matching conditions on the branes (5.4), which in these coordinates 
read 

2 = = e u , (5.37) 

are not, however, sufficient to fix all the arbitrary functions of x arising in the 
integration with respect to uj. In fact, two arbitrary functions of x remain in 
the solution, which may be regarded as time-dependent moduli describing the 
three-dimensional scale factor of the branes and their proper separation. These 
moduli may be determined with the help of the Einstein equation as we will 
demonstrate shortly. 
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Returning to (5.35) and (5.36), at y$ order now, the left-hand sides amount 
to ordinary differential equations in u for the y% corrections to (3 and v. The 
right-hand sides can no longer be neglected, but, because of the overall factor 
of ?/q, only the time-derivatives of (3 and v at zeroth order in y are involved. 
Since (3 and v have already been determined to this order, the right-hand sides 
therefore act merely as known source terms. Solving these ordinary differential 
equations then introduces two further arbitrary functions of x\ these serve as y\ 
corrections to the time- dependent moduli and may be fixed in the same manner 
as previously. 

Our integration scheme therefore proceeds at each order in y via a two-step 
process: first, we integrate the Einstein equations (5.35) and (5.36) to determine 
the a>-dependence of the bulk geometry, and then secondly, we fix the x-dependent 
moduli pertaining to the brane geometry using the G\ equation. This latter step 
works as follows: evaluating the G\ equation on the branes, we can use the Israel 
matching conditions (5.37) to replace the single ^-derivatives that appear in this 
equation, yielding an ordinary differential equation in time for the geometry on 
each brane. Explicitly, we find 



where five- dimensional considerations (see Section § 5.6) further allow us to fix 
the constants of integration on the (±) brane as 



where the brane conformal time t± is defined on the branes via ndt = bdt±. When 
augmented with the initial conditions that n and b both tend to unity as x tends 
to zero (so that we recover compactified Milne spacetime near the collision), these 
two equations are fully sufficient to determine the two x-dependent moduli to all 
orders in y . 

Putting the above into practice, for convenience we will work with the Neu- 
mann variables N and q, generalising (5.31) and (5.32) to 




(5.38) 




(5.39) 



-i 



N(x, uj) — ujx, 



b 



q(x,u). 



(5.40) 



n 



n 
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Seeking an expansion of the form 

N(x,u) = Nofa^ + ytNfa^ + Oiy*), (5.41) 
q(x,u) = q (x,uj)+ylq 1 (x,uj) + O(yQ), (5.42) 

the Einstein equations (5.35) and (5.36) when expanded to zeroth order in yo 
immediately restrict Nq and go to be functions of x alone. The bulk geometry to 
this order is then simply anti-de Sitter space with time- varying moduli, consistent 
with (5.31) and (5.32). The moduli Nq(x) and qo(x) may be found by integrating 
the brane equations (5.39), also expanded to lowest order in y . In terms of the 
Lambert W- function [92], W(x), defined implicitly by 

W(x)e w{x) = x, (5.43) 

the solution is 

jVo(x) = e* w (-*\ q Q (x) = (l + W(-x 2 )) e \ w ^ 2 \ (5.44) 

Thus we have found the scaling solution for the background, whose form is 
independent of y , holding for any y <C 1. Using the series expansion for the 
Lambert W-function about x = W(x) = 0, namely 1 

w(x) = S^ xm ' (5 - 45) 

m=l 

we can immediately check that the expansion of our solution is in exact agreement 
with (5.33) and (5.34). At leading order in y Q then, we have succeeded in resum- 
ming the Dirichlet/Neumann polynomial expansion results for the background to 
all orders in x. 

Later, we will return to evaluate the y\ corrections in our expansion about 
the scaling solution. In the next two subsections, however, we will first examine 
the scaling solution in greater detail. 

1 Note that the radius of convergence of the series (5.45) for W{x) is 1/e, and thus it 
converges for arguments in the range — 1/e < x < as required. 
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5.5 Expansion about the scaling solution 



5.5.2 Evolution of the brane scale factors 

Using the scaling solution (5.44) to evaluate the scale factors on both branes, we 
find to 0{yl) 

b ± = 1 ± xe~^ w{ - x2) = 1 ± ^-W(-x 2 ). (5.46) 

To follow the evolution of the brane scale factors, it is helpful to first understand 
the behaviour of the Lambert W-function, the real values of which are displayed 
in Figure 5.2. For positive arguments the Lambert W-function is single-valued; 
yet for the negative arguments of interest here, we see that there are in fact 
two different real solution branches. The first branch, denoted Wq(x), satisfies 
Wo 0*0 > — 1 and is usually referred to as the principal branch, while the second 
branch, W-±(x), is defined in the range W_i(x) < — 1. The two solution branches 
join smoothly at x = — 1/e, where W — — 1. 

Starting at the brane collision where x = 0, the brane scale factors are chosen 
to satisfy b± = 1, and so we must begin on the principal branch of the Lambert 
W-function for which Wo(0) = 0. Thereafter, as illustrated in Figure 5.3, b + 
increases and 6_ decreases monotonically until at the critical time x = x c , when 
Wq(—x 2 c ) = — 1 and 6„ shrinks to zero. From (5.43), the critical time is therefore 
x c = = 0.606..., and corresponds physically to the time at which the negative- 
tension brane encounters the bulk black hole 1 . 

At this moment, the scale factor on the positive-tension brane has only at- 
tained a value of two. From the Birkhoff-frame solution, in which the bulk is 
AdS-Schwarzschild and the branes are moving, we know that the positive-tension 
brane is unaffected by the disappearance of the negative-tension brane and sim- 
ply continues its journey out to the boundary of AdS. To reconcile this behaviour 
with our solution in brane-static coordinates, it is helpful to pass to t + , the con- 
formal time on the positive-tension brane. Working to zeroth order in y , this 
may be converted into the dimensionless form 

x 4 = ^± = ^/ ?dt=/4^- = *e-^- 2 ) = y^W{^). (5.47) 
L L J b J q (x) 

1 From the exact solution in bulk-static coordinates, the scale factor on the negative-tension 
brane at the horizon obeys b 2 _ — scch2Y"o/^ = tanhyo, and so 6_ ~ j/q ' . 
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Figure 5.2: The real values of the Lambert W-function. The solid line indicates 
the principal solution branch, Wq(x), while the dashed line depicts the W-i(x) 
branch. The two branches join smoothly at x = — 1/e where W attains its 
negative maximum of —1. 
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Inverting this expression, we find that the bulk time parameter x = i 4 e"^. 
The bulk time x is thus double-valued when expressed as a function of £4, the 
conformal time on the positive-tension brane: to continue forward in £ 4 beyond 
X4 = 1 (where x = x c ), the bulk time x must reverse direction and decrease 
towards zero. The metric functions, expressed in terms of x, must then continue 
back along the other branch of the Lambert W-function, namely the branch. 
In this manner we see that the solution for the scale factor on the positive-tension 
brane, when continued on to the W-\ branch, tends to infinity as the bulk time 
x is reduced back towards zero (see dotted line in Figure 5.3), corresponding to 
the positive-tension brane approaching the boundary of AdS as £4 — > 00. 

For simplicity, in the remainder of this chapter we will work directly with the 
brane conformal time £4 itself. With this choice, the brane scale factors to zeroth 
order in yo are simply b± — 1 ± £4. 



5.5.3 Analytic continuation of the bulk geometry 

In terms of £4, the metric functions n and b are given by 

n = -^^ + 0{yll b = ±—^ + 0{yl). (5.48) 
1 — ujx^ 1 — UJX^ 

At £4 = 1, the three-dimensional scale factor b shrinks to zero at all values of uj 
except uj = 1 (i.e. the positive-tension brane). Since b is a coordinate scalar under 
transformations of £4 and uj, one might be concerned that that the scaling solution 
becomes singular at this point. When we compute the y$ corrections, however, 
we will find that these corrections become large close to £4 = 1, precipitating a 
breakdown of the small-?/o expansion. Since it will later turn out that the scaling 
solution maps directly on to the four- dimensional effective theory, and that this, 
like the metric on the positive-tension brane, is completely regular at £4 = 1, 
we are encouraged to simply analytically continue the scaling solution to times 
£4 > 1. 

When implementing this analytic continuation careful attention must be paid 
to the range of the coordinate uj. Thus far, for times £4 < 1, we have regarded uj as 
a coordinate spanning the fifth dimension, taking values in the range — 1 < uj < 1. 
The two metric functions n and b were then expressed in terms of the coordinates 
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Figure 5.3: The scale factors b± on the positive-tension brane (rising curve) and 
negative-tension brane (falling curve) as a function of the bulk time parameter x, 
to zeroth order in y . The continuation of the positive-tension brane scale factor 
on to the W_\ branch of the Lambert W-function is indicated by the dashed line. 
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£4 and uj. Strictly speaking, however, this parameterisation is redundant: we 
could have chosen to eliminate uj by promoting the three-dimensional scale factor 
b from a metric function to an independent coordinate parameterising the fifth 
dimension. Thus we would have only one metric function n, expressed in terms 
of the coordinates £4 and b. 

While this latter parameterisation is more succinct, its disadvantage is that 
the locations of the branes are no longer explicit, since the value of the scale 
factor b on the branes is time-dependent. In fact, to track the location of the 
branes we must re-introduce the function uj{x^b) = (b + x\ — l)/bx4 (inverting 
(5.48) at lowest order in y ). The trajectories of the branes are then given by the 
contours uj — ±1. 

The contours of constant uj as a function of £4 and b are plotted in Figure 
5.4. The analytic continuation to times £4 > 1 has been implemented, and the 
extent of the bulk is indicated by the shaded region. From the figure, we see 
that, if we were to revert to our original parameterisation of the bulk in terms of 
£4 and uj, the range of uj required depends on the time coordinate £4: for early 
times £4 < 1, we require only values of uj in the range — 1 < uj < 1, whereas 
for late times £4 > 1, we require values in the range uj > 1. Thus, while the 
positive-tension brane remains fixed at uj — 1 throughout, at early times £4 < 1 
the value of uj decreases as we head away from the positive-tension brane along 
the fifth dimension, whereas at late times £4 > 1, the value of uj increases away 
from the positive-tension brane. 

While this behaviour initially appears paradoxical if uj is regarded as a coordi- 
nate along the fifth dimension, we stress that the only variables with meaningful 
physical content are the brane conformal time £4 and the three-dimensional scale 
factor b. These physical variables behave sensibly under analytic continuation. 
In contrast, uj is simply a convenient parameterisation introduced to follow the 
brane trajectories, with the awkward feature that its range alters under the ana- 
lytic continuation at £ 4 = 1. 

For the rest of this chapter, we will find it easiest to continue parameterising 
the bulk in terms of £4 and uj, adjusting the range of the uj where required. Figure 
5.5 illustrates this approach: at early times £4 < 1 the three-dimensional scale 
factor b is plotted for values of uj in the range — 1 < uj < 1. At late times £4 > 1, 
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0.5 1 1.5 2 2.5 

Figure 5.4: The contours of constant uj in the (b, x±) plane. Working to zeroth 
order in yo, these are given by x^ = \{buj ± \Jb 2 uj 2 — 4(6 — 1)) , where we have 
plotted the positive root using a solid line and the negative root using a dashed 
line. The negative-tension brane is located at uj — — 1 for times x^ < 1, and 
the trajectory of the positive-tension brane is given (for all time) by the positive 
root solution for uj = 1. The region delimited by the trajectories of the branes 
(shaded) then corresponds to the bulk. From the plot we see that, for < x± < 1, 
the bulk is parameterised by values of uj in the range — 1 < u < 1. In contrast, 
for £4 > 1, the bulk is parameterised by values of uj in the range uj > 1. 
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Figure 5.5: The three-dimensional scale factor 6, plotted to zeroth order in y Q as a 
function of £4 and u, for x 4 < 1 (top) and £4 > 1 (bottom). The positive-tension 
brane is fixed at u — 1 for all time (note the evolution of its scale factor is smooth 
and continuous), and for x 4 < 1, the negative-tension brane is located at uj — — 1. 
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we must however plot b for values of uj in the range uj > 1. In this fashion, the 
three-dimensional scale factor b always decreases along the fifth dimension away 
from the brane. 

We have argued that the scaling solution for the background, obtained at 
lowest order in yo, may be analytically continued across X4 — 1. There is a coor- 
dinate singularity in the x^, uj coordinates but this does not affect the metric on 
the positive-tension brane which remains regular throughout. The same features 
will be true when we solve for the cosmological perturbations. The fact that the 
continuation is regular on the positive-tension brane, and precisely agrees with 
the predictions of the four-dimensional effective theory, provides strong evidence 
for its correctness. Once the form of the the background and the perturbations 
have been determined to lowest order in y , the higher-order corrections are ob- 
tained from differential equations in y with source terms depending only on the 
lowest order solutions. It is straightforward to obtain these corrections for x± < 1. 
If we analytically continue them to x± > 1 as described, we automatically solve 
the bulk Einstein equations and the Israel matching conditions on the positive 
tension brane for all x±. The continued solution is well behaved in the vicin- 
ity of the positive-tension brane, out to large distances where the yo expansion 
eventually fails. 

5.5.4 Higher-order corrections 

In this section we explicitly compute the y^ corrections. The size of these cor- 
rections indicates the validity of the expansion about the scaling solution, which 
perforce is only valid when the y$ corrections are small. 

Following the procedure outlined previously, we first evaluate the Einstein 
equations (5.35) and (5.36) to 0{%/q) using the ansatze (5.41) and (5.42), along 
with the solutions for Nq(x) and qo(x) given in (5.44). The result is two second- 
order ordinary differential equations in uj, which may straightforwardly be inte- 
grated yielding Ni(x,u) and qi(x,u>) up to two arbitrary functions of 24. These 
time-dependent moduli are then fixed using the brane equations (5.39), evaluated 
at O(yo) higher than previously. 
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Xa 



To 0(i/q), we obtain the result: 

I z a I x 2 2 

n(X4 ' W) = I^T + on/ ; 7 ^ X 7T 2 x4 (^(5^ (~3 + O; 2 ) 

1 — a>X4 30(— l+c;x 4 ) (— I+X4) 

-5x 4 + 40c; (-3 + c; 2 ) x 2 - 5 (-14 + 9c; 2 (-2 + c; 2 )) x^ 

+3c; 3 (-5 + 3c; 2 ) x^ - \§x\ + 5x1) - 5(-l + xl) 3 ln(l - , 

(5.49) 

1 ^ I-WX4 30(-l+u;x4) 2 (-l+x 2 ) 3 1 1 ' 
-20x 4 + 5c; (-7 + 4c; 2 ) xl - 10 (l - 12c; 2 + 3c; 4 ) 

+3c; (-20 - 5c; 2 + 2c; 4 ) x\ - 12x^ + 31c;x^ - 5c;x^ 

-5(-l + xl) 2 (u - 2x 4 + c;x^) ln(l -xl)). (5.50) 

In Figure 5.6, we have plotted the ratio of the corrections to the corre- 
sponding terms at leading order: where this ratio becomes of order unity the 
expansion about the scaling solution breaks down. Inspection shows there are 
two such regions: the first is for times close to X4 = 1, for all u>, and the second 
occurs at late times X4 > 1, far away from the positive-tension brane. In neither 
case does the failure of the yo expansion indicate a singularity of the background 
metric: from the bulk-static coordinate system we know the exact solution for 
the background metric is simply AdS-Schwarzschild, which is regular everywhere. 
The exact bulk-static solution in Birkhoff frame tells us that the proper speed of 
the negative-tension brane, relative to the static bulk, approaches the speed of 
light as it reaches the event horizon of the bulk black hole. It therefore seems 
plausible that a small-yo expansion based on slowly moving branes must break 
down at this moment, when X4 = 1 in our chosen coordinate system. 

Analytically continuing our solution in both X4 and u around X4 = 1, the 
logarithmic terms in the y$ corrections now acquire imaginary pieces for times 
X4 > 1. Since these imaginary terms are all suppressed by a factor of y 2 , however, 
they can only enter the Einstein-brane equations (expanded as a series to y$ 
order) in a linear fashion. Hence the real and imaginary parts of the metric 
necessarily constitute independent solutions, permitting us to simply throw away 
the imaginary part and work with the real part alone. As a confirmation of 
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Figure 5.6: The ratio of the corrections to the leading term in the small-yo 
expansion for b, plotted for x 4 < 1 (top) and £4 > 1 (bottom), for the case where 
yo = 0.1. Where this ratio becomes of order unity the expansion about the scaling 
solution breaks down. The analogous plots for n display similar behaviour. 
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this, it can be checked explicitly that replacing the In (1 — x 4 ) terms in (5.49) 
and (5.50) with In 1 1 — x±\ still provides a valid solution to 0{y^) of the complete 
Einstein-brane equations and boundary conditions. 

Finally, at late times £4 > 1, note that the extent to which we know the bulk 
geometry away from the positive-tension brane is limited by the y$ corrections, 
which become large at an increasingly large value of uj, away from the positive- 
tension brane (see Figure 5.6). The expansion about the scaling solution thus 
breaks down before we reach the horizon of the bulk black hole, which is located 
at uj — > 00 for X4 > 1. 

5.5.5 Treatment of the perturbations 

Having determined the background geometry to 0(jJq) in the preceding subsec- 
tions, we now turn our attention to the perturbations. In this subsection we 
show how to evaluate the perturbations to 0(yfy by expanding about the scaling 
solution. The results will enable us to perform stringent checks of the four- 
dimensional effective theory and moreover to evaluate the mode-mixing between 
early and late times. 

In addition to the dimensionless variables x = yod/L and uj = y/yo, when 
we consider the metric perturbations we must further introduce the dimension- 
less perturbation amplitude B = By^c 2 / L 2 ~ BV 2 /L 2 and the dimensionless 
wavevector k = kL/y ~ ckL/V. In this fashion, to lowest order in y Q and k, we 
then find = A — B/t 2 = A — B/x 2 and similarly kct = kx. (Note that the 
perturbation amplitude A is already dimensionless however). 

Following the treatment of the perturbations in the Dirichlet /Neumann poly- 
nomial expansion, we will again express the metric perturbations in terms of Wl 
(obeying Dirichlet boundary conditions), and the Neumann variables 04 and £4, 
defined in (5.20). Hence we seek an expansion of the form 



4 (x 4 ,u;) 
£4(2:4, u) 

W L (X4,UJ) 



04o(^4, uj) + yl 04i (x 4 , w) + 0(j/q), 
£40(2:4,^) +Vq €*l(x4,u) +O(yo), 
HV^, ^) + Vo W L i(x 4 , u) + 0(yi) 



(5.51) 
(5.52) 
(5.53) 
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As in the case of the background, we will use the Gf equation evaluated on the 
brane to fix the arbitrary functions of £4 arising from integration of the Einstein 
equations with respect to uj. By substituting the Israel matching conditions into 
the G5 equation, along with the boundary conditions for the perturbations, it is 
possible to remove the single ^-derivatives that appear. We arrive at the following 
second-order ordinary differential equation, valid on both branes, 

= 2n[x\- l)(2b 2 <p A + ^)b 2 + b 2 (nx 4 (x^ - 3) - {x\- l)n) (26 2 4 - £ 4 ) 

+bb Unx 4 (x 2 4 - 3)(6 2 04 + &) - (4 ~ 1)(4(& 2 04 + &)n - n(lO6 2 4 + £4))) 
+bn(x 2 - 1) (4(& 2 4 + U)b + 26 3 4 - &| 4 ) , (5-54) 

where dots indicate differentiation with respect to £4, and where, in the interests 
of clarity, we have omitted terms of 0(k 2 ). 

Beginning our computation, the Gf and G\ Einstein equations when evaluated 
to lowest order in y immediately restrict 4O and £40 to be functions of x 4 only. 
Integrating the G° equation with respect to to then gives Wlo in terms of 04 O 
and £ 40 , up to an arbitrary function of x 4 . Requiring that Wlo vanishes on both 
branes allows us to both fix this arbitrary function, and also to solve for £40 in 
terms of 040 alone. Finally, evaluating (5.54) on both branes to lowest order in 
Ho and solving simultaneously yields a second-order ordinary differential equation 
for 040, with solution 

/3A 9B\ B 

where the two arbitrary constants have been chosen to match the small-i series 
expansion given in Section §5.3. With this choice, dropping terms of 0(k 2 ) and 
higher, 

U = -A + lf-^ + (A- 3 f) X l (5.56) 
W L0 = (1 1^2)1 (3^ 2 (-2 + ux 4 ) + Bx 4 {^x 4 + u{l - jjx 2 ))) e~K (5.57) 

The resulting behaviour for the perturbation to the three-dimensional scale factor, 
b 2 ^L, is plotted in Figure 5.7. 
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Figure 5.7: The perturbation to the three-dimensional scale factor, b 2 ^> l, plotted 
on long wavelengths to zeroth order in y for early times (top) and late times 
(bottom). Only the B mode is displayed {i.e. A = and B = 1). Note how the 
perturbations are localised on the positive-tension brane (located at u — 1), and 
decay away from the brane. 
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In terms of the original Newtonian gauge variables, an identical calculation 
(working now to all orders in k but dropping terms of 0{y 2 )) yields, 

$l = 2 k S 2 ~ (A Jo(kx 4 ) + B Y (kx 4 )) 

O \X 4 1 J 

+ {l ~^ X f ) (A Ji(kx 4 ) + BoY^kx,)), (5.58) 

^ J- *^4 

m L = 1 (2fc(i - ux 4 ) 2 (A j (kx 4 ) + ^ y (fcx 4 )) 

3(l-x|)V 

-3 (x4 + ^(-2 + ^X4))(AoJi(fcr 4 ) + B Yi(fcc 4 ))), (5.59) 
W L = 2x\e-? x * (fS x 2)2 (k(l-uxJ(AoJ (kx4)+B Y (kx A )) 

+uj {AMkxi) + B Q Y y {kx 4 ))\ (5.60) 
where the constants A and 5 are given by 

A ° = T " + o 5fc dn2 - 7 ) + 0(y 2 ), B = — + 0(y 2 ). (5.61) 
k 8k * 4 

To evaluate the y 2 corrections, we repeat the same sequence of steps: inte- 
grating the G\ and G\ Einstein equations (at y\ higher order) gives us 4i and 
£41 up to two arbitrary functions of 24, and integrating the G° equation then 
gives us Wlx in terms of these two arbitrary functions plus one more. Two of the 
three arbitrary functions are then determined by imposing the Dirichlet bound- 
ary conditions on Wli, and the third is found to satisfy a second-order ordinary 
differential equation after making use of (5.54) on both branes. Solving this dif- 
ferential equation, the constants of integration appearing in the solution are again 
chosen so as to match the small-t series expansion of Section §5.3. 

Converting back to the original longitudinal gauge variables, the results to 
O(yo) an d t° 0(k 2 ) take the schematic form 

$l = fo+Vo(fi +f! In (l+x 4 ) + / 3 *ln(l- x 4 ) +ft In (l-u>x 4 ), (5.62) 

*l = fo+Vo(fi +fi In (l + x 4 ) + / 3 * In (1- x 4 ) + ff In (1-cjx 4 ), (5.63) 
W L = e-te(f™ + y 2 (/r + fY ln(l + x 4 ) 

+f™ In (1 - x 4 ) + /f In (1 - ux 4 ))), (5.64) 
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where the / are rational functions of X4 and oj which, due to their length, have 
been listed separated separately in Appendix C.3. (If desired, more detailed 
results including the 0{k 2 ) corrections are available [88]). 

It is easy to check that the results obtained by expanding about the scaling 
solution are consistent with those obtained using our previous method based 
upon Dirichlet/ Neumann polynomials. Taking the results from the polynomial 
expansion given in Appendix C.2, substituting t = (x4/yo)e~2 x * and y = uy , 
retaining only terms of 0(yfy or less, one finds agreement with the results listed 
in Appendix C.3 after these have been re-expressed as a series in X4. This has 
been checked explicitly, both for the background and the perturbations. 

Just as in the case of the background, the small- yo expansion breaks down for 
times close to £4 = 1, when the y$ corrections to the perturbations become larger 
than the corresponding zeroth order terms. Again, we will simply analytically 
continue the solution in £4 and uo around this point. In support of this, the 
induced metric on the positive-tension brane is, to zeroth order in y , completely 
regular across x 4 = 1, even including the perturbations as can be seen from (5.58) 
and (5.59). 

As in the case of the background, any imaginary pieces acquired from ana- 
lytically continuing logarithmic terms are all suppressed by order y$. Thus they 
may only enter the Einstein-brane equations (when these are expanded to order 
2/o) in a linear fashion, and hence the real and imaginary parts of the metric con- 
stitute independent solutions. We can therefore simply drop the imaginary parts, 
or equivalently replace the In (1 — £4) and In (1 — uox^) terms with In |1 — x 4 | and 
In 1 1 — ux^\ respectively. We have checked explicitly that this still satisfies the 
Einstein-brane equations and boundary conditions. 

5.6 Comparison with the four-dimensional ef- 
fective theory 

We have now arrived at a vantage point from which we may scrutinise the predic- 
tions of the four-dimensional effective theory using our expansion of the bulk ge- 
ometry about the scaling solution. We will find that the four-dimensional effective 
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theory is in exact agreement with the scaling solution. Beyond this, the y$ correc- 
tions lead to effects that cannot be described within a four-dimensional effective 
framework. Nonetheless, the higher-order corrections are automatically small at 
very early and very late times, restoring the accuracy of the four-dimensional 
effective theory in these limits. 

In the near-static limit, the mapping from four to five dimensions may be cal- 
culated from the moduli space approach [2, 93, 94]: putting the four-dimensional 
effective theory metric into Einstein frame, the mapping reads 

g+ v = cosh 2 (</>/V%)L g- v = smh 2 (4>/V6)g*„, (5.65) 

where and g~ are the metrics on the positive- and negative-tension branes 
respectively, and <fi is the radion. As we showed in Chapter § 4, on symmetry 
grounds this is the unique local mapping involving no derivatives [31], and that 
to leading order, the action for g and is that for Einstein gravity with a 
minimally coupled scalar field. 

Solving the four- dimensional effective theory is trivial: the background is 
conformally flat, = fef (£4) r)^, and the Einstein-scalar equations yield the 
following solution, unique up to a sign choice for <f> , of the form 

b\ = C 4 U, ev 7 !^ = A 4 t 4 , (5.66) 

with 0o the background scalar field, and A 4 and C 4 arbitrary constants. (Through- 
out this chapter we adopt units where 87rG*4 = 1). 

According to the map (5.65), the brane scale factors are then predicted to be 

b± = -he~$s ± ev 7 !^ j = 1 ± A 4 U, (5.67) 

where we have chosen C4 = AA^, so that the brane scale factors are unity at 
the brane collision. As emphasised in [46], the result (5.67) is actually exact for 
the induced brane metrics, when t 4 is identified with the conformal time on the 
branes. From this correspondence, one can read off the five-dimensional meaning 
of the parameter A 4 . it equals L _1 tanhyo (our definition of yo differs from that 
of [46] by a factor of 2). 
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With regard to the perturbations, in longitudinal gauge (see e.g. [34, 41]) the 
perturbed line element of the four- dimensional effective theory reads 

dsl = bl(t 4 ) [-(1 + 2$ 4 )d^ + (1 - 2$ 4 )df 2 ] , (5.68) 

and the general solution to the perturbation equations is [10, 46] 

$4 = Y 4 (AiMktJ+BoY^ktS), (5-69) 
-| = | k (A J (kU) + B Y (ktS) (AoJi(kt 4 ) + BoY^ktS), (5.70) 

with A and B being the amplitudes of the two linearly independent perturbation 
modes. 

5.6.1 Background 

In the case of the background, we require only the result that the scale factors 
on the positive- and negative-tension branes are given by 

b± = l±A 4 t i , (5.71) 

where the constant A^ = L~ x tanh yo and £4 denotes conformal time in the four- 
dimensional effective theory. (Note this solution has been normalised so as to 
set the brane scale factors at the collision to unity). Consequently, the four- 
dimensional effective theory restricts b + + 6_ = 2. In comparison, our results 
from the expansion about the scaling solution (5.50) give 



2 x\ [x 4 3) i/q 



b + + 6_ = 2 + *-±L£± + 0{y ,y (5 . 72) 



Thus, the four-dimensional effective theory captures the behaviour of the full 
theory only in the limit in which the yfi corrections are small, i.e., when the 
scaling solution is an accurate description of the higher-dimensional dynamics. 
At small times such that 14 C 1, the corrections will additionally be suppressed 
by 0(2:4), and so the effective theory becomes increasingly accurate in the Kaluza- 
Klein limit near to the collision. Close to 14 = 1, the small- yo expansion fails, 
hence our results for the bulk geometry are no longer reliable. For late times 
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£4 > 1, the negative-tension brane no longer exists and the above expression is 
not defined. 

We can also ask what the physical counterpart of £4, conformal time in the 
four- dimensional effective theory, is: from (5.50), to 0(i/q), we find 

2A 4 y 30(1 - x\) V / 

(5.73) 

In comparison, the physical conformal times on the positive- and negative-tension 
branes, defined via bdt± = ndt = (n/yo)(l — x 2 ) e~ x ^ 2 dx^, are, to 0(uq), 

t , = — + y° ( 10 - 30x1 - x 3 A5 - Uxl + hx\) + 5x 4 (l - xlf ln(l - x 2 S) 

Vo 30(l-x2) 3 V 4 V 

(5.74) 

and 



( 10 - 30x1 + xl(5- Uxl + 5x 5 4 ) -5x^(1- xlf \n(l-xl) J, 



yo 30(1 - x 4 

(5.75) 

where we have used (5.49) and (5.50). 

Remarkably, to lowest order in y , the two brane conformal times are in agree- 
ment not only with each other, but also with the four-dimensional effective theory 
conformal time. Hence, in the limit in which y 2 corrections are negligible, there 
exists a universal four-dimensional time. In this limit, t± = Xa/vq and the brane 
scale factors are simply given by b± = 1 ± ^4^4 = 1 ± x 4 . The four-dimensional 
effective scale factor, 6 4 , is given by 

[b A f = b\-b 2 _= 4A 4 t 4 = 4x 4 h = 2x] /2 . (5.76) 

In order to describe the full five-dimensional geometry, one must specify the 
distance between the branes d as well as the metrics induced upon them. The 
distance between the branes is of particular interest in the cyclic scenario, where 
an interbrane force depending on the interbrane distance d is postulated. In the 
lowest approximation, where the branes are static, the four-dimensional effective 
theory predicts that 

b+ \ , . ( 1 + A 4 t 4 



d = Llncoth|^j=Ll„^j=L.„^ T - 1 ^|. (5.77) 
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Substituting our scaling solution and evaluating to leading order in y , we find 



d = L\n^-^)+O(y 2 ). (5.78) 

(Again, this quantity is ill-defined for x± > 1). 

In the full five-dimensional setup, a number of different measures of the inter- 
brane distance are conceivable, and the interbrane force could depend upon each 
of these, according to the precise higher-dimensional physics. One option would 
be to take the metric distance along the extra dimension 

ryo ryo rl 1 2 

d m — L / y/g~y~y~({y = L ntdy — L / nx 4 e~^ Xi duj. (5.79) 
J -yo J —yo J— l 

Using (5.49), we obtain 

d m = L [ —^duj + 0{yl) = Lin (\±^±\ + O (y 2 ), (5.80) 

in agreement with (5.78). 

An alternative measure of the interbrane distance is provided by considering 
affinely parameterised spacelike geodesies running from one brane to the other at 
constant Birkhoff- frame time and noncompact coordinates x l . The background 
interbrane distance is just the affine parameter distance along the geodesic, and 
the fluctuation in distance is obtained by integrating the metric fluctuations along 
the geodesic, as discussed in Appendix C.4. One finds that, to leading order in 
yo only, the geodesic trajectories lie purely in the ^-direction. Hence the affine 
distance d a is trivially equal to the metric distance d m at leading order, since 

d a = L J \J g ab x a x b d\ = L J ntydX = d m , (5.81) 

where the dots denote differentiation with respect to the affine parameter A. 
Both measures of the interbrane distance therefore coincide and are moreover 
in agreement with the four- dimensional effective theory prediction, but only at 
leading order in y . 
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5.6.2 Perturbations 

Since the four-dimensional Newtonian potential $4 represents the anticonformal 
part of the perturbed four-dimensional effective metric (see (5.68)), it is unaffected 
by the conformal factors in (5.65) relating the four-dimensional effective metric 
to the induced brane metrics. Hence we can directly compare the anticonformal 
part of the perturbations of the induced metric on the branes, as calculated in 
five dimensions, with 2$4 in the four-dimensional effective theory. The induced 
metric on the branes is given by 



ds 2 = b 2 (-(l + 2$ L )d4 + (l-2* L )df 2 ) 

= b 2 (1 + $ L - tf £ ) (-(1 + $ L + * L ) dt 2 ± + (l-y L - $ L ) dx 2 ) , (5.82) 



where the background brane conformal time, t±, is related to the bulk time via 
bdt± = ndt. The anticonformal part of the metric perturbation is thus simply 
$£, + It is this quantity, evaluated on the branes to leading order in y , that 
we expect to correspond to 2$4 in the four- dimensional effective theory. 

Using our results (5.51) and (5.55) from expanding about the scaling solution, 
we have to 0{uq), 



with Aq and B$ as given in (5.61). On the other hand, the Newtonian potential 
of the four- dimensional effective theory is given by (5.69). Since t±, the conformal 
time in the four-dimensional effective theory, is related to the physical dimension- 
less brane conformal time £4 by £4 = 24/7/0 (to lowest order in 7/0) > an d moreover 
k = k/yo, we have kx± = H4. Hence, the four-dimensional effective theory predic- 
tion for the Newtonian potential is in exact agreement with the scaling solution 
holding at leading order in y , upon identifying A with A Q /y and B with B /y . 
The behaviour of the Newtonian potential is illustrated in Figure 5.8. 

Turning our attention now to the radion perturbation, 5<f), we know from our 
earlier considerations that this quantity is related to the perturbation Sd in the 
interbrane separation. Specifically, from varying (5.77), we find 



-($ L + ^ L )+ = -($ L + ^ L )_ = <p i0 (x 4 ) 




1 



A Ji(fcx 4 ) +B Y 1 (kx 4 )) , (5.83) 




(5.84) 
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Figure 5.8: The four-dimensional Newtonian potential $ 4 on the positive-tension 
brane, plotted to zeroth order in yo as a function of the time £4 for wavelength 
k = 1. The upper plot illustrates the mode with A = 1 and B = 0, while the 
lower plot has A = and B = 1. 
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Inserting the four- dimensional effective theory predictions for and 5<f>, we obtain 
Sd ( 4A 4 t 4 \ /2 . , ~ , . s ~ ^ 1 



-J ^- fc (A Jo(H 4 ) + B Y (kU)) - - (A Ji(kt 4 ) + B Y x {kt 4 ) 



(5.85) 

where to lowest order A 4 = y + 0(?/q). 

In comparison, the perturbation in the metric distance between the branes is 

5d m f yo f 1 nx 4 £ 4 _i x 2 

—— = J n tT L dy = j — — e ^du, (5.86) 

where we have used (5.20). Evaluating the integral using (5.56), to an accuracy 
of 0(]Jq) we obtain 

5d m f 1 nx 4 £ 40 (x 4 ) _i x 2 2x 4 ^ 40 (x 4 ) 

1 '_ p 2 4Hm = . 

1 & 2 (i-^!) 2 

1 / 8 

kx 4 (A J (kx 4 ) + B Y (kx 4 )) 



(x 2 4 -l)\3 



-4(A Ji{kx 4 ) + B Y 1 (kx 4 ))^j , (5.87) 



which is in agreement with (5.85) when we set A 4 t 4 ~ y t 4 = x 4 , along with 
A = A /y , B = -Bo/2/0 and k = kyo- The calculations in Appendix C.4 show 
moreover that the perturbation in the affine distance between the branes, 5d a , is 
identical to the perturbation in the metric distance Sd m , to lowest order in y . 

The four-dimensional effective theory thus correctly predicts the Newtonian 
potential $4 and the radion perturbation 5(j), but only in the limit in which the 
2/q corrections are negligible and the bulk geometry is described by the scaling 
solution. While these corrections are automatically small at very early or very 
late times, at intermediate times they cannot be ignored and introduce effects 
that cannot be described by four- dimensional effective theory. The only five- 
dimensional longitudinal gauge metric perturbation we have not used in any of 
the above is Wl' this component is effectively invisible to the four-dimensional ef- 
fective theory, since it vanishes on both branes and has no effect on the interbrane 
separation. 
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5.7 Mixing of growing and decaying modes 

Regardless of the rapidity of the brane collision y , one expects a four- dimensional 
effective description to hold both near to the collision, when the brane separation 
is much less than an AdS length, and also when the branes are widely separated 
over many AdS lengths. In the former case, the warping of the bulk geometry is 
negligible and a Kaluza-Klein type reduction is feasible, and in the latter case, 
one expects to recover brane-localised gravity. At the transition between these 
two regions, when the brane separation is of order one AdS length, one might 
anticipate a breakdown of the four-dimensional effective description. 

When the brane separation is of order a few AdS lengths, however, the 
negative-tension brane reaches the horizon of the bulk black hole and the small- 
yo expansion fails. This failure hampers any efforts to probe the breakdown of 
the four- dimensional effective theory at X4 = 1 directly; instead, we will look 
for evidence of mixing between the four- dimensional perturbation modes in the 
transition from Kaluza-Klein to brane-localised gravity. 

To see this in action we have to compare the behaviour of the perturbations 
at very small times with that at very late times: in both of these limits a four- 
dimensional effective description should apply, regardless of the collision rapidity 
yo, in which the four-dimensional Newtonian potential $4 satisfies 



$ 4 = — ( A Ji(kx 4 ) + B Y 1 (~kx i )) • (5.* 



■i' 1 

Expanding this out on long wavelengths k <C 1, taking in addition foe 4 C 1, we 
find 

$ 4 = - % + Al + - B% 2 In kx 4 - - A°~k 2 x 2 d + C^fo 4 ), (5.89) 
x% 2 8 

where the dimensionless constants A\ and B\ are given in terms of the five- 
dimensional perturbation amplitudes A and B by 

A °, = -A-—B--Bk 2 , B° 4 = -B, (5.90) 
4 2 16 8 4 2 ' v ; 

where we have used (5.61) and (5.83), recalling that $ 4 = </> 40 at leading order in 



103 



5. Solution of a braneworld big crunch/big bang cosmology 



In comparison, using our results from expanding about the scaling solution, 
we find the Newtonian potential on the positive-tension brane at small times 
x 4 <C 1 is given by 

" A + I A - ^ B - \ Bk 2 + \ Bk 2 In kx A - A Ak 2 x\ + fiP^ 
2x4 2 16 8 4 16 128 

+l/o ( — — -B — 3 Ar 4 — — Bx A — — Bk 2 x/± In fcx 4 + 6AX4 

+ Bx\ - ^ Sfc 2 ^ + Bk 2 x\ In /cx 4 ] f O(xl) + O(y 4 ). 
105 960 

Examining this, we see that to zeroth order in yo the result is in exact agreement 
with (5.89) and (5.90). At y 2 order, however, extra terms appear that are not 
present in (5.89). Nonetheless, at sufficiently small times the effective theory is 
still valid as these 'extra' terms are subleading in x 4 : in this limit we find 

B E 1 

$ 4 = - -§- + Af + - Bf k 2 In ~kx A + 0(x 4 ) + 0(Ai 4 ) (5.92) 
X4 2 

(the superscript E indicating early times), in accordance with the four-dimensional 
effective theory, where 

Af = Al + ^Byl Bf = Bl (5.93) 

At late times such that x± 3> 1 (but still on sufficiently long wavelengths that 
kxi <C 1), we find on the positive-tension brane 

- -A + I A - — B - - Bk 2 + \B~k 2 In k Xi - — A~k 2 x\ + — Bk 2 x\ 
2x1 2 16 8 4 16 4 128 4 

A B Ak 2 llZBk 2 Bk 2 \nk Bk 2 In x± 
+Vo I ~~ 3xf ~ 24xf ~ ~8xj + 960x1 18^1 12x1 

~B + 2 -Bk 2 + l -Bk 2 \nx 4 + ^Bk 2 xij + O Qj) +O(y 4 ). (5.94) 

To zeroth order in y Q , the results again coincide with the effective theory pre- 
diction (5.89) and (5.90). At y 2 order, however, extra terms not present in the 
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four-dimensional effective description once more appear. In spite of this, at suf- 
ficiently late times the effective description still holds as these 'extra' terms are 
suppressed by inverse powers of x^ relative to the leading terms, which are 

$ 4 = A\ - % - - A\ k 2 x\ + OCk 2 In kxA (5.95) 
x\ 8 

(where the superscript L indicates late times), in agreement with the four-dimen- 
sional effective theory. (Since x& 3> 1, we find k <C kx^ <C 1, and we have chosen 
to retain terms of 0{k 2 x 2 i ) but to drop terms of 0(k 2 ). The term of 0(x7 2 ) is 
much larger than 0(k 2 ) and so is similarly retained). Fitting this to (5.94), we 
find 

A\ = Al - | Byl Bl = Bl + (| + ^)y 2 . (5.96) 

Comparing the amplitudes of the two four-dimensional modes at early times, 
Af and with their counterparts A\ and B\ at late times, we see clearly that 
the amplitudes differ at y\ order. Using (5.90), we find: 



A\\ = ( 1 -Wo 

KhZ 1 + d + fiH 2 




(5.97) 



Hence the four-dimensional perturbation modes (as defined at very early or very 
late times) undergo mixing. 



5.8 Summary 

In this chapter we have developed a set of powerful analytical methods which, we 
believe, render braneworld cosmological perturbation theory solvable. 

Considering the simplest possible cosmological scenario, consisting of slowly 
moving, flat, empty branes emerging from a collision, we have found a striking 
example of how the four- dimensional effective theory breaks down at first non- 
trivial order in the brane speed. As the branes separate, a qualitative change 
in the nature of the low energy modes occurs, from being nearly uniform across 
the extra dimension when the brane separation is small, to being exponentially 
localised on the positive-tension brane when the branes are widely separated. If 
the branes separate at finite speed, the localisation process fails to keep up with 
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the brane separation and the low energy modes do not evolve adiabatically. In- 
stead, a given Kaluza-Klein zero mode at early times will generically evolve into 
a mixture of both brane-localised zero modes and excited modes in the late-time 
theory. From the perspective of the four-dimensional theory, this is manifested 
in the mixing of the four- dimensional effective perturbation modes between early 
and late times, as we have calculated explicitly. Such a mixing would be impossi- 
ble were a local four-dimensional effective theory to remain valid throughout cos- 
mic history: mode-mixing is literally a signature of higher- dimensional physics, 
writ large across the sky. 

The strength of our expansion about the scaling solution lies in its ability to 
interpolate between very early and very late time behaviours, spanning the gap in 
which the effective theory fails. Not only can we solve for the full five-dimensional 
background and perturbations of a colliding braneworld, but our solution takes us 
beyond the four-dimensional effective theory and into the domain of intrinsically 
higher- dimensional physics. 
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Chapter 6 

Generating brane curvature 



If you can look into the seeds of time, 

And say which grain will grow, and which will not... 

MacBeth, Act I. 

6.1 Introduction 

A key quantity of cosmological interest is (, the curvature perturbation on co- 
moving slices. Well away from the collision, £4, the curvature perturbation in the 
four-dimensional effective theory, agrees well with the same quantity ( defined on 
the branes. Yet whereas the brane curvature perturbation is precisely conserved 
(in the absence of additional bulk stresses), the curvature perturbation in the 
four-dimensional effective theory changes as the branes approach at order (V/c) 2 , 
due to the mixing of four- dimensional perturbation modes. 

In light of this failure of the four- dimensional effective theory, it is interesting 
to revisit the generation of curvature on the branes from a five- dimensional per- 
spective. In the ekpyrotic and cyclic models, this takes place through the action 
of an additional bulk stress, AT|, on top of the background negative cosmologi- 
cal constant [1, 2, 9]. In this chapter, we show how ( on the branes is generated 
as the two branes approach by an effect proportional to AT| times an 'entropy' 
perturbation. The latter measures the relative time delay between the interbrane 
distance and a quantity measuring the time delay on the branes. The entropy 
perturbation is nonzero even in the model with no AT|, for which we have solved 
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the five-dimensional Einstein equations. We are thus able to give an expression 
for the final brane curvature which is accurate to leading order in AT|. The 
iV/c) 2 violations of the four-dimensional effective theory cause the entropy per- 
turbation to acquire a scale-invariant spectrum, which is then converted, under 
the influence of AT|, into a scale-invariant brane curvature perturbation. 



6.1.1 C m the four-dimensional effective theory 

Working in longitudinal gauge, from (2.24) the comoving curvature perturbation 
in the four-dimensional effective theory is given by 

C4 = $4 ~ W-W =3 Aa > ((U) 

where primes denote differentiation with respect to the four-dimensional con- 
formal time £4. (Recall from (5.76) that the comoving Hubble parameter "K = 
h'Jb^ = l/2t 4 , and that $4 = A 4 — B±/t\ on long wavelengths). From the four- 
dimensional effective theory, we expect that £4 is constant on long wavelengths. 
As a consequence of the mode-mixing in (5.97), however, the value of A4 (and 
hence £4) differs at 0{uq) between early and late times. 



6.1.2 ( on the branes 

To find a quantity that is exactly conserved on long wavelengths, let us instead 
define the comoving curvature perturbation on the branes directly in five dimen- 
sions. The following linear combination of the five- dimensional isotropic spatial 
metric perturbation \1/ and the lapse perturbation $ is gauge-invariant, when 
evaluated on either brane: 

C = tt- + HQ), (6.2) 

where T is now the proper time on the brane and H = h^jb. In four-dimensional 
cosmological perturbation theory, the three-curvature of spatial slices is propor- 
tional to k 2 ^, and the second term in (6.2) is proportional to the matter velocity 
with respect to those slices [41]. The latter may be gauged away in a suitable 
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comoving time-slicing. Thus, k 2 ( may be interpreted as the curvature pertur- 
bation of spatial slices which are comoving with respect to the matter. In our 
case, no matter is present on the branes. Nevertheless, ( is still a well-defined 
gauge-invariant quantity. For empty branes and an empty bulk (apart from the 
background negative cosmological constant), as we shall see below, £ is a con- 
served quantity at long wavelengths. We can therefore calculate £ at small t from 
the results (5.6) and (5.11) from the previous chapter. Evaluating (6.2) on each 
brane, we find that, to order y^, 

C + = C- = 2A- \Byl -\ B yt = \A\ = \a e a - g^f, (6.3) 

where the latter two equalities follow from (5.96), (5.90) and (5.97). Thus, in the 
limit of late times, (± and £4 coincide. At early times, however, £4 differs from 
C± at O(y 2 ). 

6.2 Inducing brane curvature from a five-dimen- 
sional perspective 

In the ekpyrotic and cyclic models, the primordial ripples on the branes are 
generated by a force between the branes [2, 10]. In the four- dimensional effective 
theory, this is described by an assumed potential V(4>) for the radion field <fi. In 
the five-dimensional description, the force between the branes is associated with 
an additional bulk stress AT|, on top of the background negative cosmological 
constant. This extra stress enters the background evolution equations as follows. 
The background G\ Einstein equation, when evaluated on each brane with the 
appropriate boundary conditions [89], yields an equation for the Hubble constant 
H on each brane [58]: 

H T + 2H 2 = --^ AT* (6.4) 

where AT| is the five-dimensional stress evaluated on the brane, and the five- 
dimensional Planck mass m 5 is defined so that coefficient of the Ricci scalar in the 
five-dimensional Einstein-Hilbert action is ml/2. Likewise for the linearised per- 
turbations, the G5 equation, when supplemented with the appropriate boundary 
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conditions, yields [89] 

k 2 1 

* TT + H ($ T + 4* T ) + 2$ (H T + 2H 2 ) = — ($ - 2*) + — - <5Tf, (6.5) 

where 5T| is the linearised perturbation of the bulk stress [89]. 

From equations (6.2), (6.4) and (6.5), one can derive the following evolution 
equation for £: 

Ct = -^3-^ [(AT|), T S-A ; 2 r 2 i/($-2M/)], (6.6) 

where we define a gauge-invariant 'entropy' perturbation 

H ST 5 
§ = + + # — -f— , (6.7) 

which measures the difference between the dimensionless 'time delay' in 5Tf 
and a quantity transforming in an identical manner constructed from the four- 
dimensional metric perturbations. Therefore § can be thought of as a measure 
of the extent to which JTf is synchronised with the metric perturbations on the 
brane. For adiabatic perturbations in four dimensions, § is zero on long wave- 
lengths for all scalar variables such as ATf. 

In a fully five- dimensional description of the attractive interbrane force, (6.6) 
and (6.7) could be used to study how £ is generated on the branes. Such a 
calculation is beyond the scope of the present work. Instead, we will content 
ourselves with an estimate of the induced ( on the branes, under the assumption 
that AT| is a function of the interbrane separation d. Since § does not depend on 
the magnitude of ATf (nor indeed on its functional dependence upon d), as long 
as ATf is small we can still use our solution from the preceding chapter, in which 
the effects of ATf have been ignored, to reliably calculate S to lowest order in 
ATf. Nonetheless, some ambiguity remains in our result because it is not obvious 
a priori precisely which measure of the interbrane distance d should be used. In 
fact, this will depend on the detailed microphysics generating the interbrane force, 
which will be model-dependent. We will investigate two possible cases here and 
use the difference between the answers as a measure of the uncertainty in the 
final result. 
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As a first choice, let us assume that the bulk stress AT| is a function of the 
geodesic separation between the two branes, along spacelike geodesies chosen to 
be orthogonal to the four translational Killing vectors of the static background 
(corresponding to shifts in x and Birkhoff- frame time). The relevant distance d 
and its perturbation are calculated in the final section of Appendix D. We find 
that, as the collision approaches, S tends to a constant on both branes: 

S = A + #$) + H ^ « CsB, (V 2 /L 2 ) {V/cf + 0(t 2 ), (6.8) 
£i t a t 

where V is the relative speed of the branes as they collide, is the growing 
mode amplitude in the four-dimensional effective theory 1 , and C,s = 7 for this 
choice of d. We have also calculated § under the assumption that the relevant d 
is just the naive metric separation between the two branes, J dy^fg^, in which 
case we find Cs — 1. 

The absence of any fc-dependence in this formula means that in the cyclic 
model, where the coefficient B± has a scale-invariant spectrum in the incoming 
state, the entropy perturbation is scale-invariant as the two branes collide. We 
emphasise that S is a fully gauge-invariant quantity, measuring the relative time 
delay between the perturbations on the branes, and the time to the collision as 
measured by the distance d, upon which the interbrane force depends. 

We can see from (6.6) how a bulk stress AT|, acting on the entropy pertur- 
bation S, leads directly to a scale-invariant £ on the branes. If we work pertur- 
batively in AT|, then using (6.4) and (6.6), to leading order in k we find 

fjkh^*- (6 - 9) 

When the branes are close, § is approximately constant. Since the integrand is 
not a total derivative, and therefore does not cancel, it follows that ( acquires 
the same ^-dependence as S. For example, taking AT 5 5 to be a narrow negative 
square well of strength A and width r, we find 



1 From (5.90) and (5.93), to leading order in (V/c) 2 we have = B/2 as the collision 
approaches, where B is the five-dimensional growing mode amplitude. 
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to leading order in r, where So is the constant value in (6.8). 

In the cyclic model, however, AT| is not small. Its presence, as discussed in [1], 
qualitatively alters the evolution of the background solution so that the negative 
tension brane, initially contracting after the collision, turns around and starts 
expanding so that the interbrane distance never becomes large. An accurate 
calculation of the amplitude of the final perturbation spectrum will therefore 
require the effect of AT| on the background, and on the perturbations, to be 
included. 

It is nevertheless not hard to estimate the parametric dependence of the result. 
In the scaling solution of the four- dimensional effective theory, one has H ~ T _1 

and m^ 3 AT| T~ 2 from (6.4). Equation (6.6) then reads T(d(/dT) ~ §. Well 

before the collision, when the branes are moving slowly with respect both to each 
other and the bulk, one expects the four- dimensional effective description to be 
accurate and hence § should be zero. The brane curvature is then directly related 
to the four- dimensional effective theory curvature perturbation £ 4 , which is zero 
on long wavelengths (hence A4 = 0). As the branes approach, their speed picks 
up and § rises to a value of order ~ B^(V 2 / L 2 )(V/c) 2 . Finally, as the branes 
become very close, AT| turns off and £ is, from (6.6), conserved once more. 

On long wavelengths, therefore, we obtain the following equation for the brane 
curvature as the collision approaches: up to a factor of order unity, we have 

C+((r) « C-(O-) ~ B 4 (V 2 /L 2 )(V/c) 2 . (6.11) 

In the ekpyrotic and cyclic models B 4 is scale-invariant in the incoming state, 
and so the final-state spectrum of curvature perturbations on the brane is also 
scale- invariant. 
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Colliding branes in heterotic 
M-theory 

What is it that breathes fire into the equations 
and makes a universe for them to describe? 

Hawking, A Brief History of Time. 

Our goal in this chapter is to derive a cosmological solution of the Hofava- 
Witten model with colliding branes, in which the five-dimensional geometry about 
the collision is that of a compactified Milne spacetime, and the Calabi-Yau volume 
at the collision is finite and nonzero. We construct this solution as a perturbation 
expansion in the rapidity of the brane collision, applying the formalism developed 
in Chapter § 5. As this work is currently still in progress, we will present here 
only the leading terms in this expansion, corresponding to the scaling solution 
for the background geometry. 

7.1 Heterotic M-theory 

One of the most successful models relating M-theory to phenomenology is that of 
Hofava and Witten [95], in which an eleven- dimensional spacetime is the product 
of a compact Calabi-Yau manifold with a five- dimensional spacetime consisting 
of two parallel 3-branes or domain walls, one with negative and the other with 
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positive tension. As in the case of the Randall- Sundrum model, the dimension 
normal to the 3-branes is compactified on an S l /Z2 orbifold. 

Performing a generalised dimensional reduction of the eleven- dimensional the- 
ory, such that the 4-form has non- vanishing flux on the Calabi-Yau internal space 
[96], one obtains a five-dimensional supergravity theory as shown in [97] 1 . The 
five-dimensional action for the bulk metric and the scalar field <f), representing 
the volume of the internal Calabi-Yau manifold, then takes the form 

S = jj=g[R- i(<90) 2 - 6aV^] - £ £ Jg± o±e~\ (7.1) 

where the sum runs over the two branes, with induced metrics g^ u , and we have 
dropped the Gibbons-Hawking term. The constants a^ 1 are given by a ± = ±12a, 
where a is an arbitrary constant parameterising the amount of 4-form flux thread- 
ing the Calabi-Yau. The exponential potential in the bulk action is the remnant 
of the 4-form field strength, and plays a role analogous to that of the bulk cos- 
mological constant in the Randall- Sundrum model. In the above, and for the 
remainder of this chapter, we have adopted units in which the five-dimensional 
Planck mass is set to unity. 

A static solution of this model [97], comprising two domain walls located at 
y = ±y , is given by 

ds 2 = Hiy^dx^dx" + H\y)dy\ (7.2) 
= H 3 (y), (7.3) 
H(y) = l + 2a\y + y \, (7.4) 

where, to create the second domain wall, we must additionally impose a Z 2 reflec- 
tion symmetry about y = +y - Recently, exact solutions of the Horava-Witten 
model with colliding branes were found in [98]. An analysis of these solutions 
from a four- dimensional effective perspective was given in [59]. In these solu- 
tions, however, the volume of the internal Calabi-Yau manifold shrinks to zero 
at the collision. For the purposes of constructing an M-theory model of a big 
crunch/big bang transition, it seems preferable to seek solutions in which the 

1 Notc in particular that the truncation from eleven to five dimensions is consistent, i.e., a 
solution of the five-dimensional theory is an exact solution of the full eleven-dimensional theory. 
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volume of the Calabi-Yau manifold at the collision is nonzero. This way, the 
only singular behaviour at the collision is that of the five-dimensional spacetime 
geometry, which we will take to be Milne about the collision. 

7.2 A cosmological solution with colliding branes 
7.2.1 Equations of motion 

As in the case of the Randall-Sundrum model 1 , cosmological symmetry on the 
branes permits the bulk line element to be written in the form 

ds 2 = n 2 (t, y)(-dt 2 + t 2 dy 2 ) + b 2 (t, y)dx 2 , (7.5) 

where the x % coordinates span the flat, three-dimensional spatial world volume of 
the branes. Moreover, in this coordinate system, the branes are held fixed at 
y = ±y . (Recall also that, physically, y represents the rapidity of the brane 
collision occurring at t — 0). As in Chapter §5, it is useful to introduce the 
variables uj = y/yo and x = y t, along with the metric functions j3(x, uj) = 3\nb 
and u(x, uj) = ln(nx). 

The junction conditions, evaluated at uj — ±1, are then 

u i = ae u -*, f3' = 3ae u -*, 0' = 6ae^, (7.6) 

where primes denote differentiation with respect to uj. Evaluating the 0, Gq + G^, 
G^, Gq + G\ — (1/2)G*, and G 05 equations in the bulk, we obtain the following 
set of equations: 

y 2 (4> + P<P), (7.7) 
y 2 (P + P 2 ), (7.8) 
y 2 o0+lp 2 -Pv + \<p 2 ), (7.9) 

yl{i>- l -P 2 + \tf\ (7.10) 

0, (7.11) 



<P" + P'(f>' + 12a 2 e 2v ~ 2 * = 
{3" + {3' 2 + 6a 2 e 2v - 2 * = 
l -(3' 2 + l3'v' -^ + 3a 2 e 2 ^ = 

3 1 A Y 
iSee Section §5.2. 



115 



7. Colliding branes in heterotic M-theory 



where, throughout this chapter, we will use a dot as a shorthand notation for 
xd x . 

In the above, both the equation (7.9) and the G05 equation (7.11) involve 
only single derivatives with respect to u. Applying the junction conditions, we 
find that both left-hand sides vanish when evaluated on the branes. The G05 
equation is then trivially satisfied, while the equation yields the relation 

h \ _ b,x n x 1 2 
b + b 2 bn + 12^ U ' {< ] 

valid on both branes. Introducing the brane conformal time r, defined on either 
brane via bdr = ndx, this relation can be re-expressed as 

T + ^ = ' (7 ' 13 > 

7.2.2 Initial conditions 

Through a suitable rescaling of coordinates, we can always arrange for the metric 
functions n and b to tend to unity as t — > 0, so that the geometry about the 
collision is that of a Milne spacetime. Furthermore, we are seeking a solution in 
which the Calabi-Yau volume tends to a finite, nonzero value at the collision. 
Since a constant shift in is a symmetry of the system (given that a is arbitrary), 
without loss of generality we can set — > 1 as t — > 0. 

To check that such a solution is possible, and to help ourselves later with fixing 
the initial conditions, it is useful to solve for the bulk geometry about the collision 
as a series expansion in t. Up to terms of order t 3 , the solution corresponding to 
the Kaluza-Klein zero mode is: 

a 2 

n = 1 + a (secht/o sinhy) t + — sech 2 ?/ (9 — cosh2y — 8 cosh 2y) t 2 , (7.14) 

8 

o? 

b = 1 + a (sechj/ sinhy) t + — sech 2 y (3 + cosh 2y — 4 cosh 2y) t 2 , (7.15) 

= 1 + 6a (sechy sinh?/) t — sech 2 y (2 — cosh2y — cosh2y) t 2 , (7.16) 

where we have used the junction conditions to fix the arbitrary constants arising 
in the integration of the bulk equations with respect to y. 
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The brane conformal times are then given by 



T± = 




dx = x + 0(x 3 ), 



(7.17) 



in terms of which the brane scale factors b± are 



b± = 1 ± a tanh y ( — 



) 



do i 2 ( T ± 

- a tanh y — 
2 \y 



) 




3 



(7.18) 



7.2.3 A conserved quantity 

In the case of the Randall- Sundrum model, where there were only two time- 
dependent moduli, the G\ equation evaluated on both branes provided sufficient 
information to determine both moduli. In the present setup, however, we have 
an extra time-dependent modulus due to the presence of the scalar field 0. We 
therefore require a further equation to fix this. 

Introducing the variable x — 4> ~ 2/3, upon subtracting twice (7.8) from (7.7) 
we find 



which is simply the two-dimensional massless wave equation Ox — 0. Since the 
junction conditions imply x' = on the branes, the left-hand side vanishes upon 
integrating over uj. A second integration over x then yields 



for some constant 7, which we can set to zero since our initial conditions are such 
that xe 13 -> — > Let us now consider solving (7.19) for x, as a perturbation 
expansion in y^. Setting x = Xo + UoXi + 0(Vo)i an d similarly for j3, at zeroth 
order we have (x'o e ^ )' = 0. Integrating with respect to to introduces an arbitrary 
function of x which we can immediately set to zero using the boundary condition 
on the branes, which, when evaluated to this order, read x'o = 0- This tells us 
that x'o = throughout the bulk; xo is then a function of x only, and can be 
taken outside the integral in (7.20). Since 7 = 0, yet the integral of e 13 across the 
bulk cannot vanish, it follows that xo must be a constant. 

At order y^, the right-hand side of (7.19) evaluates to y\ (xo^°)\ which van- 
ishes. Evaluating the left-hand side, we have (Xi e/3 °)' — 0, and hence, by a 




(7.19) 




(7.20) 
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sequence of steps analogous to those above, we find that Xi must also be con- 
stant. It is easy to see that this behaviour continues to all orders in y\. We 
therefore deduce that \ = <fi — 2/3 is exactly constant. Since both and (3 tend 
to zero as x — > 0, this constant must be zero, and so we find <p = 2(3. 

The essence of this result is that the scaling solution, and consequently any 
expansion about it in powers of y$, exists only when x is i n the Kaluza-Klein zero 
mode. This may be seen from (7.19), which reduces, in the limit where x — > and 
[3 — > 0, to x" = UoX- F° r a perturbative expansion in y\ to exist, the right-hand 
side of this equation must vanish at leading order. This is only the case, however, 
for the Kaluza-Klein zero mode: all the higher modes have a rapid oscillatory 
time dependence such that the right-hand side does contribute at leading order. 
This means that for the higher Kaluza-Klein modes a gradient expansion does 
not exist. 

Setting <fi = 2/3 from now on, returning to (7.13) and recalling that f3 — 3 In 6, 
we immediately obtain the equivalent of the Friedmann equation on branes; 

K " +4f%) 2 = 0. (7.21) 



b J,T V b 

Integrating, the brane scale factors are given by 

b± = A ± (T ± -c ± y/\ (7.22) 

To fix the arbitrary constants A± and c±, we need only to expand the above in 
powers of t± and compare with (7.18). We find 

b ± = {l±Ar ± fl\ (7.23) 

where A = (4a /yo) tanhyo- This equation determines the brane scale factors to 
all orders in y^ in terms of the conformal time on each brane, and is one of our 
main results. (As a straightforward check, it is easy to confirm that the 0(t±) 
terms in (7.18) are correctly reproduced). 

7.2.4 The scaling solution 

Having set up the necessary prerequisites, we are now in a position to solve the 
bulk equations of motion perturbatively in y^. Setting 

P = Po + VoPi + O(y 4 ), u = u + ylv x + O(y 4 ), (7.24) 
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the scaling solution corresponds to the leading terms in this expansion, i.e., to 
Po and vq. To determine the tu-dependence of these functions, we must solve the 
bulk equations of motion at zeroth order in y^. Evaluating the linear combination 
(7.8) — 3 [(7.9) + (7.10)], noting that at this order the right-hand sides all vanish 
automatically, we find 

m - Hy°y = o. (7.25) 

Using the boundary condition j3' — 3u' = on the branes, we then obtain /3q = 
3u + f(x), for some arbitrary function f(x). Backsubstituting into (7.9) and 
taking the square root then yields 

v' Q = ae- 5u °- 2f , (7.26) 

where consistency with the junction conditions (7.6) forced us to take the positive 
root. Integrating a second time, and re-writing = B~ 5 , we find 



e 



v ° = B\x)h l '\ = B(x)h 3/5 , (7.27) 



where h = 5au + C(x) and C(x) is arbitrary. 

In the special case where B and C are both constant, we recover the exact 
static domain wall solution (7.2), up to a coordinate transformation. In general, 
however, these two moduli will be time dependent. Inverting the relation b\ = 
B 5 / 3 (±5a + C) to re-express B and C in terms of the brane scale factors b±(x), 
we find 

B V3.J_ = 5a(£±y. (7.28) 

Furthermore, at zeroth order in t/q, the conformal times on both branes are equal, 
since 

dr = ? dx = e v °-M 3 — = B^ 3 {x)— (7.29) 
b xx 

is independent of to. To this order then, (7.23) reduces to 

b± = (l±4ar) 1/4 , (7.30) 

allowing us to express the moduli in terms of r as 

5 5 / 3 (r) = ^[(l+4ar) 5/4 -(l-4ar) 5 / 4 ], (7.31) 
10a 



C(r) = 5a 



(1 +4ar) 5 / 4 + (1 -4ar) 5 / 4 
(1 + 4ar) 5 /4 _ (1 -4aT) 5 /4 



(7.32) 
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The brane conformal time r and the Milne time x are then related by 



As we do not know how to evaluate this integral analytically, it seems that the 
best approach is simply to adopt r as our time coordinate 1 in place of x. The 
complete scaling solution is then given by (7.27), (7.31) and (7.32). 

7.3 Discussion 

In this chapter we have taken some first steps towards a cosmological solution of 
the Hofava-Witten model with colliding branes, in which the Calabi-Yau volume 
at the collision is finite and nonzero, and the five-dimensional spacetime geometry 
about the collision is Milne. Employing the techniques developed in Chapter §5, 
we have shown how to construct the scaling solution for the background geometry. 
All the relevant tools are now in place to compute the corrections at 0(uq). Only 
two steps are required: firstly, integrating the bulk equations of motion at 0(yfy 
to find the ^-dependence of the solution; and secondly, fixing the time dependence 
using our explicit knowledge of the scale factors on the branes, according to (7.23). 
With this accomplished, it is then feasible to consider solving for cosmological 
perturbations about this background, along the lines of the work presented in 
Chapter § 5. 

An interesting feature of our solution is the vanishing of the scale factor (and 
also the volume of the Calabi-Yau manifold) on the negative-tension brane at 
r = (l/4a). In more general models incorporating matter on the negative-tension 
brane, four- dimensional effective theory arguments suggest that under similar cir- 
cumstances the scale factor on the negative-tension brane, rather than shrinking 
to zero, would instead undergo a bounce at some small nonzero value. Since this 
behaviour persists even in the limit of only infinitesimal amounts of matter on the 
negative-tension brane, one might consider implementing a bounce in the present 
model by replacing the factors of (1 — 4ar) in (7.31) and (7.32) with |1 — 4ar\. 

^n support of this, note that the relation between x and r is monotonic. At small times 




(7.33) 



x = t + {a 2 /4)t 3 + 0(r 5 ). 
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The scale factor on the positive-tension brane, along with the bulk time x, would 
then continue smoothly across the bounce and out to late times as r — > oo. The 
implications of this continuation remain to be explored. 
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Conclusions 

All of this danced up and down, like a company of gnats, 
each separate, but all marvellously controlled in an invisible 
elastic net - danced up and down in Lily's mind... 

Virginia Woolf, To the Lighthouse. 

In this thesis we have explored the impact of higher-dimensional physics in the 
vicinity of the cosmic singularity Taking a simple Randall- Sundrum braneworld 
with colliding branes as our model of a big crunch/big bang universe, we have 
explicitly calculated the full five- dimensional gravitational dynamics, both for the 
background solution, and for cosmological perturbations. Our solution method, 
a perturbative expansion in (V/c) 2 (or equivalently, in the collision rapidity), 
provides a powerful analytical tool enabling us to assess the validity of the low 
energy four-dimensional effective theory. 

While this four-dimensional effective theory holds good at leading order in our 
expansion of the bulk geometry, at order (V/c) 2 new effects appear that cannot be 
accounted for within any local four- dimensional effective theory. Our principal 
example is the mixing of four- dimensional effective perturbation modes in the 
transition from early to late times (and vice versa). 

Despite the small size of this effect in the case of a highly non-relativistic 
collision, mode-mixing is nonetheless an important piece of the puzzle concerning 
the origin of the primordial density perturbations. If these primordial fluctuations 
were indeed generated via the ekpyrotic mechanism in a prior contracting phase of 
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the universe, then mode-mixing suggests a resolution to the problem of matching 
the scale-invariant spectrum of growing mode time-delay perturbations in the 
collapsing phase, to the growing mode curvature perturbations post-bang. In 
such a scenario, higher-dimensional physics, rather than being insignificant, is 
instead responsible for all the structure we see in the universe today. 

The failure of the four- dimensional effective theory, however, reminds us that 
in order to address this issue satisfactorily, a fully five-dimensional description of 
the generation of curvature perturbations on the brane is required. Initiating this 
programme, we have shown how brane curvature is generated through the action 
of an additional five- dimensional bulk stress pulling the branes together. Our 
understanding is sufficient to provide a quantitative estimate of the final brane 
curvature prior to the collision, but to progress further it will be necessary to 
construct a specific five-dimensional model of the additional bulk stress. Work 
along these lines is already underway. 

Another direction of active research lies in the adaptation of our solution 
methods to other braneworld scenarios. Of particular interest are the cosmological 
colliding-brane solutions of the Hofava-Witten model. Here, we have constructed 
the leading terms of a background solution in which the five-dimensional geometry 
about the collision is Milne, and the Calabi-Yau volume tends to unity at the 
collision. Evaluating the corrections to the background at subleading order, as 
well as solving for cosmological perturbations, are both the subject of current 
research. Our methods should also extend to scenarios in which matter is present 
on the brane. A more challenging extension would be to probe the evolution of 
a black string [53] in an expanding cosmological background consisting of two 
separating branes. As the branes separate, the gradual stretching of the black 
string will eventually trigger the onset of the Gregory-Laflamme instability [99]. 

We conclude in the only fashion possible - with a call-to-arms. Higher- 
dimensional physics, beyond the four-dimensional effective theory, necessarily 
holds the key to the decisive experimental signature over which braneworld cos- 
mologies will one day stand or fall. 
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Appendix A 



Israel matching conditions 



In this appendix we derive the Israel matching conditions [52] from scratch, based 
on the treatment presented in [100]. 

We begin with a timelike (or spacelike) hypersurface E, pierced by a congru- 
ence of geodesies which intersect it orthogonally. Denoting the proper distance 
(or proper time) along these geodesies by £, we can always adjust our parameteri- 
sation to set £ = on £. (One side of the hypersurface is therefore parameterised 
by values of I < 0, and the other by values i > 0). Introducing n a , the unit 
normal to £ in the direction of increasing £, the displacement away from the 
hypersurface along one of the geodesies is given by dx a = n a d£. From this it 
follows that n a n a = e, where e — +1 for spacelike E, and e = —1 for timelike E, 
and that 

n a = e d a £. (A.l) 

It is also useful to introduce a continuous coordinate system x a spanning both 
sides of the hypersurface, along with a second set of coordinates y a installed on 
the hypersurface itself. (Henceforth Latin indices will be used for hypersurface 
coordinates and Greek indices for coordinates in the embedding spacetime). The 
hypersurface may then be parameterised as x a = x a (y a ), and the vectors (in a 
Greek sense) 
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are tangent to curves contained in £ (hence e^n a = 0). For displacements within 

ds| = g af3 dx a dxP (A.3) 

= ^dy'dy", (A.5) 

and so the induced metric h ab on £ is given by 

h ab = g af3 e^. (A.6) 

To denote the jump in an arbitrary tensorial quantity A across the hypersur- 
face, we will use the notation 

[A] = lim (A) - lim (A). (A.7) 

The continuity of x a and £ across £ immediately imply [n Q ] = [e°] = 0. 

To decompose the bulk metric in terms of the different metrics on either side 
of £, we will require the services of the Heaviside distribution Q(£), equal to +1 
if £ > 0; if £ < 0; and indeterminate if £ = 0. Note in particular the properties: 

2 (£) = 0(£), 0(£)0(-£) = 0, 1 0(i) = (A.8) 

where S(£) is the Dirac distribution. 

We can now express the metric g a p in terms of the coordinates x a as a 
distribution-valued tensor: 

g^ = e(£)g^ + e(-£)g-^ (A.9) 

where g^p denotes the metric on the £ > side of £, and g~^ the metric on the 
£ < side. Differentiating, we find 

gap„ = ew^ + e^j^+e^^K (a.io) 

The last term is singular; moreover, this term creates problems when we compute 
the Christoffel symbols by generating the product Q(£)5(£), which is not defined 
as a distribution. In order for the connection to exist as a distribution, we are 
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A. Israel matching conditions 



forced to impose the continuity of the metric across the hypersurface: \g a p] = 0. 
This statement can be reformulated in terms of hypersurface coordinates alone 
as = [ga,p\^a e b = [9af3 e a e b\ = i^ab], the induced metric h a b must be the 
same on both sides of E. This condition is often referred to as the 'first' junction 
condition. 

To derive the 'second' junction condition (the Israel matching condition), 
we must calculate the distribution-valued Riemann tensor. Beginning with the 
Christoffel symbols, we obtain 

r; = e(i)r}; + eH)r- (A.n) 

where are the Christoffel symbols constructed from g^. Thus 

= 0(£)r+« s + e(-£)Tfif + e5(£)[r^]n s , (A.12) 
and the Riemann tensor is 

Ufa = e(£)R + ^s + ©M)^ + m^ 5 , (A.i3) 

where 

A a ^ 5 = e ([r^K - [T%]n s ) . (A.14) 

The quantities Ap s transform as a tensor since they are the difference of two 
sets of Christoffel symbols. We will now try to find an explicit expression for this 
tensor. 

Observe that the continuity of the metric across £ in the coordinates x a 
implies that its tangential derivatives must also be continuous. Thus, if g a p n is 
to be discontinuous, the discontinuity must be directed along the normal vector 
n a . We can therefore write 

[g a p,y] = K Q/3 n 7 , (A. 15) 

for some tensor K a p (given explicitly by K a p = e[g al 3 ) . y ]n' y ). We then find 

Ry] = \ (*>7 + «> " Wl a ), (A.16) 

and hence 

Ap-yS = 7, ( K 5 n P n y - K y n p n s - np S n a n^ + n^n a n 5 ). (A. 17) 
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A few lines of calculation show us that the 5-function part of the Einstein 
tensor is 

S a p = | {K^n^np + Kftfi-n^ria - Kn a rip - EK a p - (K^n^n" - en)g a p) . (A. 18) 

On the other hand, the total stress-energy tensor is of the form 

= 6(£)T+ + + S(£)T a/3 , (A.19) 

where and T~„ represent the bulk stress-energy in the regions where £ > 
and £ < respectively, while T a/3 denotes the stress-energy localised on the 
hypersurface S itself. From the Einstein equations, we find T Q/ g = (8irG)~ 1 S a/ 3. 
It then follows that T a 0nP = 0, implying T a p is tangent to S. This allows us to 
decompose T a/3 as 

T Q/3 = T ah e a a el , (A.20) 

where T ab = T a pe^e b re-expresses the hypersurface stress-energy tensor in terms 
of coordinates intrinsic to S. From (A. 18), we have 

\QnGT ab = -K af3 e a a el - e(n^ u n^n v - en)h ab 

= -Kafi^A ~ - h mn e^ n )h ab + Kh ab 

= -^e a A + h mn ^e^ n h ab . (A.21) 
Finally, we can relate T ab to the jump in extrinsic curvature across E. From 
[V a np] = -[r^]n 7 = - {sK aP - K^npn 1 - K^n^n 7 ), (A.22) 
we deduce that 

[K ab ] = [V a n p ]eA = \ **f*%4- (A.23) 
This leads us to our goal; 

8nGT ab = -e ([K ab ] - [K]h ab ) , (A.24) 

which is the Israel matching condition (or 'second' junction condition). 



127 



Appendix B 



Five-dimensional longitudinal 
gauge 



Starting with the background metric in the form (5.1), the most general scalar 
metric perturbation can be written as [89] 

ds 2 = n 2 (-(1 + 2$)dt 2 -2Wdtdy + t 2 {\ -2Y)dy 2 ) -2V i adx i dt 

+2t 2 V i (3dydx i + b 2 ((1 - 2$) - 2V i V i x) ^ dx j . (B.l) 

Under a gauge transformation x A — > x A + £ A , these variables transform as 

$_>$_£*_£*™_£vl^ 
n n 

t n n 

W -> W - £ n + t 2 i y , 

b 2 ■ 

n z 

b b' 

X^X + r, (B.2) 
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where dots and primes indicate differentiation with respect to t and y respectively. 
Since a five- vector £ A has three scalar degrees of freedom £*, £ y and = V^ 5 , 
only four of the seven functions ($, T, W, a, (3, \&, x) are physical. We can therefore 
construct four gauge-invariant variables, which are 



$inv 


= $ - a - 


n ~n 

a p — , 

n n 


r. 

- 1 - mv 


= r + p" + 


1 h ~n' 
- a + a — hp — 
t n n 




= W-oi 






T i> ~ 

= * + - a 
b 





where a = a — (b 2 /n 2 ) x and /3 — f3 + (b 2 /n 2 t 2 ) \' ■ 

In analogy with the four-dimensional case, we then define five-dimensional 
longitudinal gauge byx = tt = /5 = 0, giving 

^inv ^ Li ^ inv - 1 - 

Winv = ^L, *inv = ^L, (B.4) 

i. e. the gauge-invariant variables are equal to the values of the metric perturba- 
tions in longitudinal gauge. This gauge is spatially isotropic in the x l coordinates, 
although in general there will be a non-zero t, y component of the metric. 

As for the locations of the branes, this will in general be different for different 
choices of gauge. In the case where the brane matter has no anisotropic stresses, 
the location of the branes is easy to establish. Working out the Israel matching 
conditions, we find that (3 on the branes is related to the anisotropic part of 
the brane stress-energy. If we consider only perfect fluids, for which the shear 
vanishes, then the Israel matching conditions give (3{y = ±yo) = 0. 

From the gauge transformations above, we can transform into the gauge a = 
X = using only a £ s and a transformation. We may then pass to longitudinal 
gauge (a — f3 — x — 0) with the transformation £ y = (3 alone. Since (3 (and 
hence (3) vanishes on the branes, £ y must also vanish leaving the brane trajectories 
unperturbed. Hence, in longitudinal gauge the brane locations remain at their 
unperturbed values y = ±yo Transforming to a completely arbitrary gauge, we 
see that in general the brane locations are given by y = ±yo — (3. 
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Appendix C 



Detailed results 

By the pricking of my thumbs, 
Something wicked this way comes. 

MacBeth, Act IV. 

In this appendix, we list the results of detailed calculations of the background 
geometry and cosmological perturbations in a big crunch/big bang universe, ob- 
tained using a variety of methods. For clarity, the AdS radius L has been set to 
unity throughout. (To restore L, simply replace t —*t/L and k — > kL). 



C.l Polynomial expansion for background 

Using the expansion in Dirichlet/Neumann polynomials presented in Section § 5.4 
to solve for the background geometry, we find 

N = I-i^ 2 + lt(8-9f 2 )y 4 -^t(l36 + 900f 2 + 375t 4 )y 6 

+ (3968 + 354816 t 2 - 348544 1 4 - 36015 1 6 ) y* + 0(yjj°) (C.l) 

N 3 = --+( — -2t 2 )y 2 — (61 - 20880 1 2 + 19440 t 4 )yt 

6 \72 J 2160 V ; y ° 

/ 277 743t 2 677 1 4 101 1 6 \ 6 ^. 8 . n , 
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C.l Polynomial expansion for background 



iV 4 = \t*y 2 + ^ 3 (-28 + 33t 2 )y 4 

+1 1 3 ( 19g4 _ 6776t 2 + 2715t 4) y e + 0(%8) (C 3) 

iV 5 = — + — — (7 - 1800 1 2 - 540 t 4 ) y 2 

5 120 1800 v Jy ° 

— (323 - 990528 t 2 + 2207520 t 4 + 362880 t 6 ) yf, + (C.4) 

201600 V > yo KyoJ V ; 

N 6 = t 3 2 / 2 + ^t 3 (-142 + 37H 2 )y 4 + O(y 6 ) (C.5) 

^7 = -^-^ (-9 + 20384 1 2 + 23520 t 4 ) Z/0 2 + O(y 4 ) (C.6) 

Ns = ^tV + 0{yi) (C.7) 

N, = ^ + 0(, 2 ) (C.8) 

A^io = 0(y 2 ) (C.9) 

and 

3 2 2 f 2 7t4 \ 4 / — 17t 2 17t 4 55£ e 



31 1 2 9t 4 233 1 6 245 1 8 , 



105 5 90 128 
3~ 



1/5 + W) (C10) 



-2t 3 ^ ('^-8t 5 )y 4 



_^_^» +27tr j rf + 0M) (cn) 

94 = ^ 4 y 2 +(^ + ^y 4 + O(, 6 ) (C.12) 

95 = -^o 2 +( Z T^-^2/o 4 + 0( % 6 ) (C.13) 
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% = \t'yl + 0{yt) (C.14) 
1 3 

9r = -^yl + 0{yt) (C.15) 

98 = O(y 2 ) (C.16) 

99 = 0(y 2 ), (C.17) 



This solution explicit satisfies all the Einstein equations up to 0(yl 



io\ 
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C. Detailed results 



C.2 Polynomial expansion for perturbations 
C.2.1 All wavelengths 

Using the Dirichlet/Neumann polynomial expansion to solve for the perturba- 
tions, the solution may be expressed in terms of the original longitudinal gauge 
variables as 

$l = ^\y,t)F^(t) +yU(y,t)FW(t), (C.18) 

where 

F^Xt) = AJ n (kt) + BY n (kt), (C.19) 

for n — 0, 1 and 7 = 0.577. . . is the Euler-Mascheroni constant. The constants 
A and B are arbitrary functions of k. In order to be consistent with the series 
expansion in t presented in Section §5.3, we must set 

A = 12A + 25fc 2 (ln2- 7 )-^ + ^M + 0( % 6 ), (C.20) 
B = Bk 2 7i + O(y 6 ). (C.21) 

The polynomials 7^ are then given (for all k and t) by 
Of (*,!/) = -^ + ^y + ^M-2y 2 + y 2 ) + lty(lV + 3(-ll + 3t 2 )y 2 ) 
+ 2160 (~ 52%4 + 90(19 " 5t2 ^ 2 yo + (511 - 180t 2 + Ahk 2 t 4 )y 4 ) 



ty 



3 (-92 + (9 + 4k 2 )t 2 ) y 4 



2160 

+30 (92 - (219 + Ak 2 )t 2 + k 2 t 4 ) y 2 y 2 

- (6900 - (20087 + 300£; 2 )t 2 + 90 (13 + k 2 ) t 4 - 90 k 2 1 6 ) y 4 
+O(y 6 ), 

(C.22) 
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C.2 Polynomial expansion for perturbations 



and by 



5-| + ^(V + (-3 + W)rf) 



| ((3 + fcV) y 2 + 3 (-3 + (3 - A; 2 ) t 2 + 2*¥) y 2 ) 



f 2 



+ 2i6o 1 75 (6 + fc2t2) 2/4 + 90 (- 12 + 3 ( 2 - k2 ) t2 + 2fcV ) y 2 y 2 o 



(718 + A; 2 t 2 (-101 + 225t 2 )) ^ 



3(3 + 2(-9 + 31A; 2 )t 2 + 2fc 2 t 4 ) y 4 



2160 

+30 (-3 + (219 - 62k 2 ) t 2 + 16k 2 t 4 ) y 2 y 2 

+ (225 - (20104 - 4650A; 2 ) t 2 + (1215 - 1822A; 2 ) t 4 + 765A; 2 t 6 ) 

+O(y 6 ). (C.23) 

Since the are of zeroth order in y , the solution for $£, to a given order less 
than O(yo) is found simply by truncating the polynomials above. (Should they 
be needed, results up to 0(?/q 4 ) can in addition be found at [88]). 
In a similar fashion we may express the solution for as 

9 L = yP(y,t)F^(t) + ^\y,t)F^(t), (C.24) 

where F^ is defined as above and 

iffry) = 1 -- t j + t ^(y 2 + y 2 ) + f 6 (-2y 2 + 3(2-3t 2 )y 2 ) 

+ 2l60 ( 12%4 + 450 ( ~ 2 + ^ ^ + (644 + 45 ° t2 " 45A;2t4) ^ 

+ 2ifo f 3 ( ~ 2 + (9 + 4fc2) ' 2) + 30 ( 2 + 4 ( 9 - fc2 ) * 2 + fc2 * 4 ) y 2 y 2 o 

- (150 + (2863 - 300A; 2 ) t 2 + 90 (13 + A; 2 ) t 4 - 90A; 2 t 6 ) y 4 
+O(y 6 ), (C.25) 
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C. Detailed results 



^\t,y) 



i 

kt 



|-^(3 2 / 2 + (3 + fcV) 2/0 2 ) 



+ S ( (3 + kH2) y2 + 3 (- 3 + ( 3 - k ") t2 + 2kH4 ) Vo) 



2160 



15 (12 + 5k 2 t 2 ) y 4 - 90 (6 - 3 (2 - A; 2 ) t 2 - 2k 2 t A ) y 2 yl 



(752 - (540 - 101& 2 ) t 2 - 360fc 2 t 4 ) y% 



ty 



(9 + 6 (-9 + k 2 ) t 2 + 6Pt 4 ) y 4 



2160 

-30 (3 + (33 + 2A; 2 ) t 2 - 7k 2 t 4 ) y 2 y 2 

+ (225 + 2 (1288 + 75A; 2 ) t 2 + (1215 - 1012A; 2 ) t 4 + 765£; 2 t 6 ) 

+O(y 6 ). 

(C.26) 



Finally, writing 



W L = ?%\y,t)F^(t)+7^(y,t)F^\t) 



(C.27) 



we find 



^ (*>*/) = -^t 2 (2/ 2 -yo)(- 30 + 30^-25( 2 / 2 +(-5 + 3t 2 )y 2 ) 
+ty (21y 2 + (-149 + 75t 2 ) y 2 ) ) + O(y 6 ), 



(C.28) 



t 2 (y 2 - y 2 ) (30y - 5k 2 t (2y 2 + (-10 + 3t 2 ) y 2 ) 



60k 

-y ((12 + llkH 2 ) y 2 + (-38 + (60 - 69A; 2 ) t 2 + IhkH 4 ) y 2 ) 
-O(y 6 ). 

(C.29) 



134 



C.2 Polynomial expansion for perturbations 



C.2.2 Long wavelengths 

On long wavelengths, reduces to 

F^\t) = 12 A _^_M + ^M + 0(fc2) + 0{yll (C30) 

F U(t) = (6At-0 + ^-^)i?^ + O(A; 2 ) + O(, o 6 ). (C.31) 
For convenience, we list below the metric perturbations truncated at 0(k 2 ). 
*l = (A-§) + (^ + Aty) + 1 -(B(y>-4y*)-4At*(y 2 + y *)) 



^ {B(3y 2 (-A + t 2 ) + Ay 2 ) + AAt 2 (3y 2 (-19 + 3t 2 ) + 1% 2 ) 



21 

+ (^At 2 (y 4 (29 - Qt 2 ) + 3y 2 (13 - 2t 2 )y 2 - lOy* 

2 -L(5(^(56 - 45t 2 ) + 15y 2 (-16 + 5*V + 100y 4 )) 
1/ 



+ -^-(2At 2 (5^(905 - 1338t 2 + 75t 4 ) + 10y 2 (-181 + 219t 2 )y 2 + 181?/ 4 ) 
+ J B(y 4 (50 - 3509t 2 + 135t 4 ) + 5y 2 (-4 + 235t 2 )y 2 + (2 - 12t 2 )y 4 )) 
+O(y 6 ), (C.32) 
¥ L = 2A-|( J B + At 2 ) + i(2^ 2 ( 2/0 2 + y 2 ) + J B(-y 2 + 2y 2 )) 

(4At 2 (3y 2 (-l + 3t 2 ) + y 2 ) + 5(3 Z / 2 (-4 + t 2 ) + Ay 2 )) 
+ ^ (8At 2 (2y 4 (17 + 3t 2 ) + 3y 2 (-7 + 2t 2 )y 2 + y*) 
+B(y*(8 + 15t 2 ) + 3y 2 (-8 + 5t 2 )y 2 + 8y 4 )) 
-^JL (2At 2 (25y 4 (l + 42t 2 + 15t 4 ) - 10y 2 (l + 39t V + y 4 ) 

+B(y%(50 + 721i 2 + 135t 4 ) - 5y 2 (4 + 47t 2 )y 2 + (2 - 12t 2 )?/ 4 )) 
+0(yg), (C33) 
W L = QAt 2 (-y 2 + y 2 )-t(B + 3At 2 )y(-y 2 + y 2 ) 

+ \t 2 (-y 2 + y 2 )(-9By 2 + 20A(y 2 (-5 + 3t 2 ) + y 2 )) 
1 -ty{-vl + ?/ 2 )(12(L4t 2 ( % 2 (-26 + 9t 2 ) + 3y 2 ) 



120 

+B(y 2 (-152 + 105t 2 ) + 48y 2 )) + O(y e ). (C.34) 
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C. Detailed results 



C.3 Perturbations from expansion about the scal- 
ing solution 

Following the method of expanding about the scaling solution presented in Section 
§5.5, we have computed the perturbations to 0{yfy. On long wavelengths, the 
five-dimensional longitudinal gauge variables take the form 

$l = fo+Vo(fi +/ 2 *ln(l + X4) + / 3 *ln(l-x 4 ) + / 4 ln(l-a;x4), (C.35) 

*l = /o I '+y 2 (A* + / 2 *ln(l + X4) + /3 I 'ln(l-X4) + /4*ln(l-^4), (C.36) 
W L = e-tt^ + tiUY + fY In (l + XA )+fW ^(1-^) 

+/f In (l-wa*)), (C.37) 

where the / are rational functions of X4 and u. For we have 

/o* = 77; — 57 — ~: orfl65- 165wx 4 -2(8A + i?-4ik> 2 )x 2 

16x|(-l + x|)\ 

+2(-8A + 3B)o;a;|+ (8A - 3I?)(3 + cj 2 )xf) , (C.38) 

/,* = c (8Ax 5 J- 580x4 + 95x1 - 576uj 5 xt + 281x1 

960x^(-l + xl) 5 V 4V 

+96u/x| - 60x1 + 5cA 4 (40 + 167x2 - 39x 4 ) + 20w 3 (-19 - 5x\ + 48x 4 ) 
-10cj 2 x 4 (78 + 117a^ - a>4 - 2 4) + 4cj(285 + 100x^ - 25x\ - 29x| + 5x1)) 
+B( - 1920 + x 4 (480cu + 80(91 + 15cu 2 )x 4 

-160cj(5 + Auj 2 )x\ + 40(-273 - 116c; 2 + 7u 4 )a% + 20cj(649 - 127c; 2 )x 4 
+4(1231 - 3285oj 2 + 2060a/)x 4 i - 4c;(2036 - 2115cj 2 + 1152cj 4 )x^ 
+ (1107 + 8102cu 2 - 4905c; 4 + 768w 6 )a£ + 12cu(233 + 48cj 2 (-5 + 3u 2 ))xl 
-(3131 + 1582cj 2 - 585a/ + 288u 6 )x 9 4 - 772uxf + 20(85 + 29uj 2 )x^ 
+180o;xf - 120(3 + cu 2 )xf))), (C.39) 

/ 2 * = 7f J B(48 + (36 -48cj)x 4 +12(-7-6cj + 2cj 2 )x^ 

48x|(-l + xf) \ 

+4(-13 + 6u + 9uj 2 )x\ + (60 + 80cj - 12c; 2 )x 4 + 4(8 + Quo - 9uj 2 )x\ 
-3(8 + 7uj + Auo 2 )x\ + (-11 + 5uo 2 )x\ + 3ua%) + 8Axf(x 4 + u; 2 x 4 
-^(1+x 2 ))), (C.40) 
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C.3 Perturbations from expansion about the scaling solution 



ff = = (j3 (- 48 + 12(3 + 4u)x 4 + 12(7 - Quj - 2uj 2 )x\ 

48xf(-l + x^) V 

-4(13 + 6u - 9u 2 )xl - 4(15 - 20u - 3c/)x 4 + 4(8 - 6u - 9u) 2 )x\ 
+3(8 - 7u + Auj 2 )x\ + (-11 + hu 2 )xl + 3uxl) 

+8Axt(x 4 + uj 2 X4-u(l + x 2 i ))), (C.41) 



ft = ™{-l + u>xtf (C.42) 
2x 4 (-l + x 2 ) 2 V ' 

For ^>l, we find 

fo = - , \ ^(l6Bou-4(8A + B(-l + 2c/))x 4 + 2(8,4 - 35>x 2 

16x 4 (— 1 + x\) \ 

+ (-8A + 35)(-3 + c/)xf), (C.43) 



f? = ■ — g ( - 480Bu - 2405(-7 + 5c/)x 4 + 1605cj(5 + 4c/)x 2 

960x|(— 1 + xf) \ 

-40.0(143 - 104c; 2 + oj a )x\ + 20cj(5(197 - 155c/) + 8A(-3 + c/))x 4 

+20(-8A(34 - 21c/ + c/) + 5(98 - 369c/ + 101u 4 ))x 5 4 

-4w(200A(-14 + 5c/) + 5(34 - 975c/ + 288u 4 ))x 6 4 

+ (40A(55 - 234c/ + 67a/) + 5(2455 + 838c/ - 1005a/ + \92u Q ))x\ 

-4oj(35(53 + 120c/ - 36c; 4 ) + 8A(155 - 120c/ + 36c/))x!j 

-(5(3515 - 1042c/ - 225c/ + 72c/) - 8,4(89 + 170c/ - 75c/ + 24c/))x 4 ' 

+4(232A + 193B)ux\° - 20(8A(1 + c/) 

+5(-91 + 29u 2 ))x\ 1 - 20(8,4 + 9B)ujx 12 + 1205(-3 + uj 2 )xf ) , (C.44) 

f? = — ? ( 8Ax\(x 4 + c/x 4 - u(l + xl)) + 5(36 + 60x 4 - 36x1 

A8x\{-l + x\Y^ 

-8Ax\ + 24x^ + 5x1 + w 2 x 4 (24 + 36x 4 - 12x^ - 36x 3 4 - 12x 4 + 5x 5 4 ) 
+cj(-48 - 72x 4 + 24x^ + 80x^ + 24x^ - 21xt + 3x1))) > (C45) 



/ 3 * = 3 (8Axl(-x 4 - c/x 4 + w(l + x\)) + £(-36 + 60x 4 + 36x^ 

48x\(-l + x 2 4 ) v 

-8Ax\ + 24x^ - 5x^ - c/x 4 (-24 + 36x 4 + 12x^ - 36x^ + 12x^ + 5x\) 
+cj(-48 + 72x 4 + 24x 2 - 80x^ + 24x 4 + 21x 5 4 - 3x1))) > ( c - 46 ) 



137 



C. Detailed results 



9 _ -3B(-l+u X4 ) 
2x|(-l + x|) 

Finally, for H^, we have 

, w (-l + ^ 2 )x4(-24Ax 4 (-2 + a;x4) + J B(-18x4 + c«;(-8 + 9xl))) . 

8(-l + xi) 

W = ("I + W ) / 8Ar 4/ 15QQ g4 4 4 ( _ 12 + 2) 

1 480x 2 (-l+x 2 ) 7 V 4V 41 4; 

-6cj 2 (50 + lOOx^ - 427x 4 + x\) + 3aA 4 (60 + 585x 2 - 160a; 4 + 3x^) 

+c/(168x 4 ' - 36x1) - 6x 4 (-590 + 243^ + 201x 4 - Ax\ + lUx\) 
+cux 4 (-1560 - 4935x2 + 2000x 4 - 2Qhx\ + 92xf)) 
+5(1440 + x 4 (-84u; 4 x 4 i (24 + 28^ + 3x 4 ) 

+12cj 5 x 4 i (40 + 6x1 + 9x 4 ) + 6cj 2 x^(-450 + 964x2 + 863x 4 + 2>x\) 
-3cj 3 x 2 (-40 - 748x 2 - 2333x 4 + 672x° + 9x^) 
-6x 4 (2160 - 5490x 2 + 3770x 4 - 2249x^ + 597x2 - 588xf + 90xf ) 
+cu(-2160 + 18920x2 - 53216x 4 + 20629x2 - 11216z| + 5579x2° 
-22364 2 + 360x^ 4 )))), (C.49) 



ff = —il — (b(- 144 + 108(-1 + 2cj)x 4 -24(-5 + o;)(3 + 2c ( j)x2 
24(x 4 - xl) K - 



— 36(— 11 + w(16 + uj))x\ + 24cj(-29 + 7cj)x 4 
+36(-2 + w(-2 + 7u))x 5 4 + 3(1 + w)(3 + 32u)x e A + 7w(l + cj)x2 
+27(1 + w)x2 - 9w(l + cj)x2) + 24,4(1 + w)x|(-l - 3x2 + cj(x 4 + x 3 )) 

(C.50) 

( 1 + UJ ) Jb( - 144+ 108(1 + 2cu)x 4 - 24(5 + u) (-3 + 2uo)x\ 

2Ax\(-l+x\) K - 

+36(-ll + (-16 + uj)uj)x\ + 24cj(29 + 7uj)x\ - 36(-2 + uj(2 + 7uj))x\ 

+3(-l + w)(-3 + 32cj)x2 - 7(-l + u)ux\ - 27(-l + cj)x2 

+9(-l + cj)^) - 24A(-1 + wJxK-l - 3x2 + cj(x 4 + xfj)) , (C.51) 



, w 3S(-1 + oux A y(A - 10x1 + wx 4 (-l + 7x1)) 

x 4 ^ 1 + X 4 J 



Results including the corrections at 0(k 2 ) can be found at [88]. 
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C.4 Bulk geodesies 

To calculate the affine distance between the branes along a spacelike geodesic we 
must solve the geodesic equations in the bulk. Let us first consider the situation 
in Birkhoff-frame coordinates for which the bulk metric is static and the branes 
are moving. The Birkhoff-frame metric takes the form (see Chapter § 3) 

ds 2 = dY 2 - N 2 (Y) dT 2 + A 2 (Y) dx 2 , (C.53) 

where for AdS-Schwarzschild with a horizon at Y = 0, 

2 cosh(2l7L) cosh(2VL) / sinh (2Y/L) \ 2 

{ ' cosh(2F /L)' 1 ' cosh(2F/L) \mx)x{2Y Q / L) ) ' 1 ' 

At T = 0, the ^-coordinate of the branes is represented by the parameter Y . 
The subsequent brane trajectories Y±(T) can then be determined by integrating 
the Israel matching conditions, which read tanh (2Y±/L) = ± a/ 1 — V± , where 
V± = (dY±/dT)/N(Y±) are the proper speeds of the positive- and negative- 
tension branes respectively. From this, it further follows that Y"o is related to 
the rapidity y of the collision by tanhy = sech(2F /L). 

For the purpose of measuring the distance between the branes, a natural choice 
is to use spacelike geodesies that are orthogonal to the four translational Killing 
vectors of the static bulk, corresponding to shifts in x and T. Taking the x and 
T coordinates to be fixed along the geodesic, we find that Y\ is constant for an 
affine parameter A along the geodesic. 

To make the connection to our original brane-static coordinate system, recall 
that the metric function b 2 (t,y) = A 2 (Y), and thus 

y j= (^+_W = „.(-«. +<M ). (C.55) 

where we have introduced the constant 9 = tanhyo = V/c. Adopting y now as 
the affine parameter, we have 

= (b 2 t b 2 + n 2 (b 4 - 6 2 ))t 2 y + 2b t b iy b 2 t, y + (b 2 y b 2 - n 2 t 2 (b A - 9 2 )), (C.56) 

where t is to be regarded now as a function of y. 
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Cl 



We can solve this equation order by order in y using the series ansatz 

oo 

t{y) = Y,c n y\ (C.57) 

n=0 

where the constants c n are themselves series in y . Using the series solution for the 
background geometry given in Section Cl, and imposing the boundary condition 
that t(yo) = t , we obtain 

co = t + t ^-2tly* + ^ + lf )y ° -tl(l + 5tl)yS 

to 17tg 5 \ 6 tg(13 + 250t§ + 795 4) ^ 
.720 + ^ + 4 *V y ° 60 
+O(y 8 ) 1 (C.58) 

2t^ 2 +(^ + 5^ y^-Ulyl 

+ ( ii + ^ + ^r) %6 + 0(2/ ° 7) ' (c - 59) 
+ (ji+^M +0(y a ); (C60) 

c 3 = -^- iig™lM + 0(^ (C.61) 

o lo 

«« = l i + (^ + z r) ! '»- i # l+0(! '» ) ' (a62) 

c 5 = 51|M + o«), (C.63) 

c 6 = -^ + O(y 2 ), (C.64) 
c 7 = + O(y ). (C.65) 



Substituting t = ^o/l/o and ?/ = uyo, we find x(o;) = £oA/o + O(yo), i.e. to 
lowest order in y 0) the geodesies are trajectories of constant time lying solely 
along the u direction. Hence in this limit, the affine and metric separation of 
the branes (defined in (5.80)) must necessarily agree. To check this, the affine 
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distance between the branes is given by 



da 

L 



/no 
nVt 2 - t' 2 dy 
-yo 



9 , , (to + 5^g) g/g 2 4 , (tp-10tg + 159tg) j/Q 5 
^ J o Uo H g 4 1 y H — 

2 (tg + 30 *g) y% (t + 31115 tg — 5523 ^ + 12795^) y£ 



3 2520 
+O(y 8 ), (C.66) 

which to lowest order in y reduces to 

T = 2xo+5 -f +5 -^ + 8 -^ + °^ + °^ (c - 67) 

in agreement with the series expansion of (5.80). (Note however that the two 
distance measures differ nontrivially at order y%). 

To evaluate the perturbation 5d a in the affine distance between the branes, 
consider 

/l f d A 
y/g^xWdX =2j ^ gp ( 6 9^ ±U + 9^,Jx K x^x u + 2g^5x u ) 



x v 8x h 



y / g p(7 xPx a 



2 J \/gpai p x a 



dA, (C.68) 



where dots indicate differentiation with respect to the affine parameter A, and 
in going to the second line we have integrated by parts and made use of the 
background geodesic equation x a = \g fJLV ^x tJ 'x v and the constraint g iiv x li x v = 1. 
If the endpoints of the geodesies on the branes are unperturbed, this expression is 
further simplified by the vanishing of the surface term. Converting to coordinates 
where to = ^0/2/0 and y = uyo, to lowest order in y the unperturbed geodesies 
lie purely in the u direction, and so the perturbed affine distance is once again 
identical to the perturbed metric distance (5.86). 
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Explicitly, we find 

6da 2 {B + Atp y 

L in 



B (2 + 2169t 2 1 + 135* l ) + 2A* 2 , (1 + lllOtg + 375t£ 



120 1 

6 



4Atg (i+42*g) -B (4 + 57^ 
1 



+ ( ^ ] Vo 



B(A + 88885tg + 952866^ + 28875*°) 



10080*o 

+4A* 2 , (1 - 152481^ + 293517^ + 36015^) ) y 7 

+O(y 8 ), (C.69) 

which, substituting * = ^o/yo and dropping terms of O reduces to 

8d a 2B A B A , 9 ~ o 25 , , 

-T = 2Ax --x -3Axg--5x^-— Ai^ 

Li Xq 4 o 4 

-^5^-^Ax 7 + O(x«), (C.70) 

where B = By%. Once again, this expression is in accordance with the series 
expansion of (5.86). At 0(yfy, however, the perturbed affine and metric distances 
do not agree. 
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